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ABSTRACT 


A  whistler  is  a  natural  audio -frequency  electromagnetic  signal  which 
.results  from  the  dispersion  of  "broad-band  energy  radiated  by  lightning 
discharges.  Whistlers  travel  along  discrete  paths  aligned  with  the  earth's 
magnetic  field.  Such  a  path  takes  the  energy  several  thousand  kilo¬ 
meters  out  into  space,  and  back  to  earth  in  the  opposite  hemisphere. 

The  space  between  the  earth  and  the  lower  edge  of  the  ionosphere  acts  as 
a  waveguide.  The  whistler  then  propagates  in  the  earth -ionosphere  wave¬ 
guide  to  the  receiver.  Thus  it  should  be  possible  to  explain  certain 
properties  of  whistlers  in  terms  of  the  properties  of  the  waveguide.  It 
should  also  be  possible  to  locate  the  point  at  which  the  energy  enters 
the  waveguide  (the  exit  point)  by  direction-finding  (d-f)  measurements  at 
spaced  stations. 

The  first  possibility  is  investigated  by  expressing  the  theoretical 
intensity,  polarization,  and  apparent  direction  of  arrival  of  the  received 
signal  as  the  sum  of  multiply -reflected  rays  propagating  in  the  earth- 
ionosphere  waveguide.  The  boundaries  of  the  guide  are  assumed  to  be 
sharply  bounded  homogeneous  media.  The  earth  is  described  by  a  complex 
dielectric  constant.  The  ionosphere  is  described  by  the  ql  approximation 
of  the  magneto-ionic  theory,  assuming  a  vertical,  static,  magnetic  field. 
The  coupling  introduced  between  the  two  polarization  components  of  the 
wave  by  the  magnetic  field  in  the  ionosphere  is  included. 

Curves  (parametric  in  frequency,  distance,  and  waveguide  parameters) 
of  the  response  of  a  vertical  monopole  and  of  crossed  loops  are  shown. 

The  ratio  of  the  crossed-loop  voltages  and  the  null -bearing  errors  are 
determined  for  comparison  with  experimental  data. 

Qualitative  agreement  is  found  between  the  calculated  and  the  pre¬ 
viously  unexplained  experimental  data  on  crossed-loop  voltage  ratios. 
Detailed  agreement  is  not  expected  because  of  the  dependence  of  the 
Lissajous  patterns  on  various  waveguide  parameters  which  are  not  known 
with  sufficient  accuracy.  It  is  possible  that  a  detailed  comparison 
could  be  obtained  from  a  large  number  of  calculations,  but  the  conditions 
would  probably  not  be  unique . 
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From  the  theoretical  work  it  was  found  that  apparent  lateral  devia¬ 
tion  of  bearing  can  be  caused  by  the  presence  of  multiple  rays,  even 
though  all  the  rays  travel  in  the  great -circle  plane  containing  the 
source  and  the  receiver. 

Another  new  result  is  that  either  a  normal  vertical  loop  or  vertical 
monopole  may  receive  a  higher  field  strength  when  multiple  rays  are 
present.  This  factor  is  significant  in  vlf  field-strength  measurements. 

A  two -station  d-f  experiment  was  performed,  using  crossed  vertical 
loops,  a  vertical  monopole,  and  a  goniometer  device.  Exit  points  were 
located  from  average  null  bearings  for  ten  whistlers.  The  geomagnetic 
latitudes  of  these  locations  were  compared  with  the  latitudes  calculated 
independently  from  whistler -dispersion  data.  The  individual  latitude 
differences  are  as  great  as  18  degrees,  although  the  average  difference 
is  only  3-1  degrees.  The  individual  differences  can  be  accounted  for  by 
the  d-f  errors  for  ea^h  location  and  therefore  do  not  contradict  the 
theory  of  whistler  propagation. 

Large  differences  in  the  null  bearing  occur  for  each  whistler  for 
different  frequencies.  Calculations  using  ray  theory  show  that  the  null¬ 
bearing  error  in  a  loop  system  varies  widely  with  frequency,  because  of 
polarization  and  apparent  lateral-deviation  error.  An  even  wider  vari¬ 
ation  in  error  is  expected  when  distance  and  waveguide  parameters  are 
varied.  The  theoretical  and  the  experimental  results  thus  agree  quali¬ 
tatively.  The  results  show  that  the  average  bearing  tends  to  approach 
the  true  bearing,  but  that  the  accuracy  is  insufficient  to  warrant 
further  d-f  work  with  the  present  system.  A  system  which  uses  only  one 
type  of  antenna,  and  preferably  responds  to  only  one  field  component,  is 
recommended.  This  system  would  eliminate  the  polarization  and  apparent 
lateral-deviation  errors. 
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I*  IMRODUCTION 


A.  BACKGROUND 

Whistlers  are  one  of  many  types  of  natural  electromagnetic  nolee 
oceurring  in  the  audible  range  of  frequencies,  They  may  frequently  be 
detected  by  listening  to  the  output  of  a  sensitive  audio  amplifier  con¬ 
nected  to  an  antenna,  Whistlers  are  gliding  tones  which  descend  in  fre¬ 
quency  across  the  audible  range  in  about  1  see.  . 

The  gliding  tone  results  from  the  dispersion  of  lightning  energy 

that  travels  along  discrete  paths  believed  to  be  aligned  with  the  earth's 

magnetic  field.  This  path  takes  the  energy  several  thousand  kilometers 

out  into  space  and  then  back  to  earth  in  the  opposite  hemisphere.  The 

action  of  free  electrons  in  the  presence  of  the  magnetic  field  causes  the 

velocity  of  propagation  to  vary  markedly  with  frequency.  At  the  lower 

frequencies  this  variation  is  given  approximately  by  D  =  t  where  t 

is  the  time  in  seconds  after  the  originating  impulse,  f  is  the  whistler 

frequency  in  cps,  and  D  is  the  dispersion.  The  value  of  D  depends  upon 

the  magnetic  latitude  of,  and  the  ionization  density  along,  the  path, 

(Recent  papers  and  reports  are  available  which  cover  the  history  and  de- 

. l-_5* 

tails  of  whistlers.) 

The  whistler  energy  leaves  the  lower  edge  of  the  ionosphere  at  a 
point  determined  mainly  by  the  location  of  such  a  path.  The  lower  edge 
of  the  ionosphere  and  the  earth's  surface  form  a  channel  which  acts  like 
a  waveguide  at  very  low  frequencies  (vlf).  The  whistler  is  propagated 
in  the  earth-ionosphere  waveguide  from  the  exit  point  to  the  receiver. 

These  considerations  lead  directly  to  a  number  of  interesting  possibilities 
for  further  research.  Two  of  these  areas  of  study  were  selected  for  dis¬ 
cussion  in  this  paper.  The  first  of  these  is  the  explanation  of  certain 
characteristics  of  whistlers  in  terms  of  the  properties  of  the  earth- 
ionosphere  waveguide;  the  second' is  the  location  of  the  exit  point  by 
direction-of -arrival  measurements  at  two  or  more  spaced  stations.  This 
information  would  provide  a  useful  check  on  the  determination  of  whistler- 

path  location  by  other  methods.  It  is  the  purpose  of  this  report  to 

*"  *  “  — 
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present  the  results  of  theoretical  and  experimental  investigations  in 
these  two  areas. 

Ray  theory  is  used  to  calculate  the  field  intensity,  apparent 
polarization,  and  apparent  direction  of  arrival  of  the  received  signal 
for  different  distances,  frequencies,  and  waveguide  parameters.  A  two- 
station  direction -finding  (d-f)  experiment  was  also  performed.  The  cal¬ 
culated  results  are  shown  to  be  consistent  with  this  and  other  experiments. 
In  addition,  the  theoretical  investigation  showed  that  the  mere  presence 
of  multiple  rays  can  cause  apparent  lateral  deviation  at  vlf  (see  p.  145). 

A  new  factor,  which  results  in  differences  in  field-etrength  measure¬ 
ments  with  vertical  loops  and  monopoles  in  the  presence  of  multiple  rays, 
was  also  found. 

B.  PREVIOUS  WORK 

The  wealth  of  experimental  data  on  whistlers  available  in  recent 

years,  especially  the  multistation  synoptic  recordings  made  during  the 

IGY  and  the  following  periods,  and  the  magneto -ionic -duct  propagation 

7 

studies,  provide  overwhelming  evidence  of  the  existence  of  discrete 

paths  or  ducts  of  propagation.  The  magnetic  latitude  of  Buch  a  path 

and  the  ionization  density  along  the  path  may.  be  calculated,  on  the 

assumption  that  the  paths  are  aligned  with  the  earth's  dipole  field.  The 

4 

results  are  quite  reasonable,  as  shown  by  Smith,  and  provide  models  of 
ionization  density  variation  with  height  along  the  paths. 

The  earliest  known  measurements  of  direction -of -arrival  and  polar¬ 
ization  of  whistlers  were  undertaken  by  the  author,  beginning  in  late 
Jl  8 

1955-  ’  In  the  first  experiment  an  attempt  was  made  to  compare  the 
amplitudes  and  phases  of  the  voltages  induced  in  a  pair  of  crossed  vertical 
loopes.  These  voltages  were  recorded  on  a  dual-channel  tape  recorder.  The 
amplifiers  and  tape-recording,  channels  were  equalized  in  amplitude  to 
within  0.1  db  and  in  phase  to  within  one  degree  over  the  useful  frequency 
range.  Lissajous  patterns  and  waveforms  were  later  photographed  on  con¬ 
tinuously  moving  film. 

The  best  results  were  obtained  during  the  great  magnetic  disturbance  ■ 
of  4  March  1956.  The  crossed-loop  antenna  voltages  produced  by  a  whistler 
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were  found  to  vary  rapidly  with  frequency  in  relative  phase  and  amplitude. 

Jfce  variations  of  axis  ratio  and  orientation  of  the  Li ssajous -pattern 

1  9 

ellipse  were  roughly  sinusoidal  in  tine,  as  previously  reported.  ’  The 
only  conclusion  that  could  be  drawn  was  that  the  direction  of  arrival 
or  polarization,  or  both,  seeaed  to  vary  rapidly  with  the  frequency  of 
the  whistler.  Both  directions  of  rotation  of  the  ellipse  were  found  to 
be  present  for  these  whistlers.  A  previous  whistler  recorded  on  13 
September  1955  at  approximately  1204  UT  was  found  to  contain  only  one 
direction  of  rotation,  although  the  crossed-loop  Lisaajotts  pattern  varied 
from  elliptical  to  linear  in  a  similar  manner.  These  results  will  be 
discussed  on  page  133. 

Some  of  the  early  work  on  the  directions  of  arrival  proceeded  from 

the  assumption  (shown  in  this  report  to  be  untenable)  that  the  energy 

would  arrive  from  the  direction  of  the  magnetic  field.  A  thesis  by 
9* 

Delloue  reported  partial  confirmation  of  this  expected  result  by  narrow- 

band  polarization  and  direction  measurements  at  5  kc.  3#  stated  that  the 

directions  of  arrival  tended  to  "circle  around  the  magnetic  field  line'' 

in  some  cases.  He  also  stated  that  the  direction  of  arrival  appears  to 

vary  rapidly  in  most  eases  and  this  wag  attributed,  to  "many  small 

elementary  whistlers,  arriving,  apparently,  from  different  directions. " 

The  author  then  stated  that  the  polarization  varied  rapidly.  He  attempted 

to  interpret  the  observed  polarization  in  terms  of  the  ionization  density 

in  the  outer  atmosphere.  Even  when  the  effects  of  the  earth -ionosphere 

waveguide  are  neglected,  this  interpretation  may  be  criticized  on  the 

35 

basis  that  the  observed  polarization  should  be  the  "limiting"  polaris¬ 
ation  characteristic  of.  the  lower  layers  of  the  ionosphere.  The  observed 
rapid  variation  of  apparent  direction  of  arrival  is  exactly  what  would  be 
expected  when  the  effects  of  the  earth -ionosphere  waveguide  are  taken 
into  account. 

The  crossed-loop  Lissajous  patterns  for  several  whistlers  were  obser¬ 
ved  by  'the  author  on  the. night  of  27  May  1959  between  0730  and  0755  UT. 
These  patterns,  sketched  directly  from  the  oscilloscope,  were  generally 


*  / 

Undated,  but  from  references  it  is  presumed  to  have  been  written 
in. 1956-7,  .  , 
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elliptical  ip  shape  and  same  showed  fairer  oops latent  patterns,  with  tome 
rapid  changes  In  tilt  or  axis  ratio**  These  results  will  be  discussed 
further  on  page  14?. 

* 

J.  M.  Watts  was  the  fleet  to  suggest  that  goniometer  techniques 

could  he  applied  to  whistlers,  and  he  made  the  first  measurements  of  this 
4l 

type.  Further  details  are  given  on  page  101, 

In  i960,  Ellis^’^  reported  4-f  measurements  of  vlf  noise  associated 
with  geomagnetic  disturbances.  The  measurements  were  made  with  crossed 
vertical -loop  antennas  and  a  goniometer.  Bearings  from  three  stations 
in  Hobart,  Tasmania,  and  Camden  and  Adelaide,  Australia,  were  reported, 
yielding  a. location  for  the  source.  Interpretations  were  made  of  the 
"null"  ratios  (maximum/minimum  voltage  ratios)  of  the  goniometer  in  terms 
of  the  size  of  the  source.  A  horizontal  direction  of  ‘arrival  was  assumed, 
thus  eliminating  the  effect  of  polarization  on  the  azimuth  determination. 
However,  this  assumption  is  not  valid,  as  the  calculations  given  later 
in  this  report  will  demonstrate. 

28  ’ 

A  different  system  is  described  by  Cartwright  who  assumes  that  the 
polarization  is  completely  random.  This  view  is  questioned,  since  a 
characteristic  polarization  is  probably  present  in  at  least  some  vlf 
noise.  An  example  is  the  "hiss"  (a  type  of  vlf  noise^)  of  4  March  19^6 
at  Stanford,  in  which  a  characteristic  pattern  was  present  in  crossed- ' 
loop  Lissajous  patterns.  The  results  given  later  will  show  that  the 
polarization  over  a  limited  frequency  range  is  not  random. 


* 

A  goniometer  is  a  device  constructed  to  provide  coupling  between 
two  circuits  and  an  output  as  a  sinc-coaine  function  of  the  rotor  position. 
It  is  constructed  with  two  pairs  of  fixed  plates  and  a  rotor  assembly, 
shaped  to  give  sinusoidal  coupling.  The  rotor  assembly  consists  of  two 
sections  that  are  90  deg  apart  in  position  on  the  shaft.  One  section 
meshes  with  each  set  of  fixed  plates,  so  that  sine-cosine  coupling  is 
obtained . 
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C.  PRESUTT  WORK 


It  has  been  shown  that  the  effect  of  the  earth- ionosphere  waveguide 
on  whistler  and  other  vlf  signals  is  of  great  interest.  By  mean*  of  ray 
theory  and  several  simplifying  assumptions  a  theoretical  method  is  developed 
that  yields  values  for  the  magnitude  and  polarization  of  the  electric 
field  and  the  apparent  angles  of  arrival  of  the  wave  as  seen  by  the  antenna 
system.  The  antenna-induced  voltages  and  certain  ratios  of  voltage*  that 
are  important  to  the  d-f  system  are  also  calculated.  The  response  of  a 
d-f  system,  utilizing  a  goniometer,  is  also  calculated  for  positions  of 
the  goniometer  for  0  to  180°  in  10°  steps.  The  results  are  presented 
in  curves,  and  are  discussed  in  comparison  with  available  experimental 
data,  which  include  those  mentioned  in  previous  work  and  the  results  of 
a  d-f  experiment  that  utilized  goniometer  systems  at  two  spaced  stations 
to  determine  exit -point  locations  from  the  observed  null  bearings. 

A  general  discussion  of  exit-point  location  by  means  of  inverse 
loran  (see  p.  100)  and  d-f  methods  for  vlf  is  presented..  In  addition,  a 
new  sensing  method  designed  for  the  goniometer  system  is  described. 

Sense  (removal  of  180°  ambiguity)  is  obtained  by  adding  the  output  of 
a  vertical  monopole  antenna  to  the  goniometer  output.  The  new  method 
also  improves  the  theoretical  bearing  accuracy  and  null  ratio  for  ellip- 
tically  polarized  waves. 

The  response  and  errors  of  the  goniometer  system,  with  and  without 
sense,  are  calculated  for  various  conditions  of  wave  polarization  and 
elevation  angle  of  arrival.  Some  of  the  results  of  an  experimental  d-f 
program  are  also  presented  and  discussed. 

The  conclusions  drawn  from  the  theoretical  results  are  given,  theoreti - 
cal  and  experimental  data  are  compared,  and  recommendations  for  future 
work  are  made. 

D.  CONTRIBUTIONS  OF  THE  PRESENT  WORK 

It  is  the  purpose  of  this  report  to  present  the  ray -theory  (geo¬ 
metrical-optics)  method  selected  for  the  calculation  of  the  effect  of 
the  earth -ionosphere  waveguide  on  whistlers  and  other  vlf  signals,  and 
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the  results  of  the  use  of  the  Model.  This  model  is  the  first  one  to  in¬ 
clude  the  total  field  excited  by  an  ellipticaJLly  polarised  wave  and  the 
effect  of  the  coupling  between  the  horizontal  (perpendicular)  and  vertical 
(parallel)  components  of  the  field  resulting  from  the  presenae  of  the 
earth's  magnetic  field.  The  results  show  that  this  relatively  simple 
model  can  explain  most  of  the  experimental  results,  including  field- 
strength  variations,  apparent  polarization,  and  apparent  arrival  angles. 
This  agreement  gives  significant  support  to  the  duct  theory  of  propagation. 
There  is  no  theory  that  vould  explain  most  of  the  experimental  results  in 
terms  of  a  signal  arriving  directly  at  the  receiver  from  the  direction  of 
the  magnetic  field.  Furthermore,  these  results  provide  for  the  first 
time  a  basis  for  evaluating  the  performance  of  antenna  systems  used  in 
whistler  observation  and  for  estimating  the  effective  area  of  a  whistler. 
They  also  provide  values  of  transmission  loss  from  the  whistler  mode  to 
the  antenna  that  may  be  used  to  evaluate  the  various  losses  in  the  ex¬ 
citation  and  transmission  of  whistler-mode  signals. 

The  equations  of  the  total  field  and  the  antenna  response  also  show 
that  the  response  of  different  antennas  is  often  not  compatible  with  the 
assumption  that  the  total  received  field  is  a  plane  wave,  although  the 
ray  components  which  make  it  up  are  essentially  plane  waves.  It  is  also 
shown  that  a  vertical  monopole  may  "see"  a  larger  vertically  polarized 
field  component  than  a  normal  loop  antenna  under  these  conditions.  The 
difference  in  the  response  of  the  normal  and  abnormal  loops  is  also  shown 

to  cause  apparent  lateral  deviation  at  vlf.  This  effect  is  similar  to,  / 

/ 

but  distinct  from,  polarization  error. 

/ 

To  test  the  hypothesis  that  whistler  signals  generally  arrive  trim 

/ 

the  direction  of  the  magnetic  field,  a  calculation  was  performed  at  a 
distance  such  that  the  n  =  0  ray  is  parallel  to  the  earth's  magnetic 
field  at  the  geomagnetic  latitude  of  Stanford  (44  N).  Tbe  characteristics 
of  the  actual  signal  calculated  for  this  case  show  that  the  polarisation 
and  the  phase  shift  between  crossed -loop  antenna  voltages  is  relatively 
constant  as  a  function  of  frequency.  These  results  show  that  the  experi¬ 
mental  results  generally  are  inconsistent  with  the  premise  that  the 
energy  arrives  from  the  direction  of  the  magnetic  field* 
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It.  THE  CALCULATION  OF  IBB  CHARACTERISTICS  OF  WAVBQUISS  ElflNAIfl 


A.  181  AVAILABLE  METSOOS  AND  THE  CHOICE  OF  A  MODEL 

It  is  now  neaeesaiy  to  develop  a  node!  for  the  calculation  of  the 

characteristies  of  whiatler-oode  eignals.  Two  general  methods  —  "ray 

1I4.  26  45 

theory"  *  and  "mode  theory"  -  have  been  uaed  in  the  calculation  of 
the  vertical  component  of  vlf  field  strengtha  for  a  path  between  two 
points  on  the  ground. 

The  mode  theory  utilized  a  aeries  of  solutions  to  Maxwell's  equation* 
matching  the  boundary  conditions  to  represent  the  wave  traveling  in  the 
guide.  Flat-earth  geometry  is  usually  assumed.  The  same  assumption* 
about  the  boundary  conditions  are  made  in  both  methods  of  solution. 

The  ray -theory  uses  the  principles  of  geometrical  optics  in  the 
solution  of  propagation  problems.  The  use  of  geometrical  optics  is 
strictly  valid  in  solving  refraction  and,  reflection  problems  only  for 
the  case  of  A  approaching  zero.  It  can  be  applied  to  other  cases  Of  re¬ 
flection,  however,  when  the  reflecting  or  refracting  medium  may  be  eon- 
sidered  to  have  slowly  varying  properties. 

Each  method  has  some  general  advantages.  Hay  theory  generally  re¬ 
quires  few  rays  for  short  distances,  while  mode  theory  requires  a  large 
number  of  modes.  At  long  distances  the  reverse  tends  to  be  true.  Thft 
exact  number  of  modes  or  rays  necessary  to  describe  the  field  for  a  given 
distance  depends,  of  course,  on  the  electrical  characteristies  of  tha 
ionosphere  and  the  ground,  and  on  the  frequency  of  the  wave.  For  long 
distances,  mode  theory  has  the  advantage  over  ray  theory  In  the  case  of 
the  calculation  of  transmission  from  point  to  point  on  the  ground.  Be¬ 
yond  the  critical  distance  (around  1000  km)  the  first  skyvave  component 
is  in  the  shadow  region  and  a  diffracted  wave  exists  at  the  receiving 
antenna.  Special  calculations  are  required  for  an  accurate  solution  for 
the  diffracted  ray  in  this  region.  This  problem  is  discussed  on  page  23* 

One  further  factor  gives  the  ray  theory  a  definite  advantage  for 
this  particular  problem.  The  majority  of  mode-theory^2  *  treatments 
deal  only  with  the  vertical  component  of  the  electric  field  of  the  wave, 
although  same  later  work  consider*  the  horizontal  component  of  the  field. 

Some  *work  by  Wait  and  Spies,  showing  that  the  curvature,  cannot  be 
neglected  in  mode  calculation*,  has  recently  been  published. 
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To  a  slightly  lesser  extent  the  same  is  true  of  the  ray  theory.  However, 
it  is  very  difficult  in  mode  theory  to  take  account  of  the  coupling 
introduced  between  the  two  polarizations  by  the  earth's  magnetic  field 
in  the  ionosphere.  This  coupling  can  be  handled  simply  in  the  ray 
theory.  The  necessary  ray-theory  expressions  for  the  complete  field  are 
much  simpler  to  derive  and  calculate,  and  reflection  coefficients  that 
incltide  the  effects  of  both  vertical  and  horizontal  field  components  are 
already  availaele.  For  these  reasons,  the  ray  theory  was  chosen  for 
this  problem. 

B.  ASSUMPTIONS 


It  is  now  necessary  to  make  several  simplifying  assumptions,  as 
follows : 

1.  Each  individual  ray  component  may  be  considered  a  plane  wave 
in  the  antenna  area,  and  the  antenna  is  considered  to  be  in  the  far 
field.  In  the  near  field  the  ratio  of  E/H  is  no  longer  constant  and  the 
antenna -response  equations  lose  accuracy*  Thus  the  problem  is  restricted 
to  distances  beyond  the  order  of  one  half  the  wavelength.  This  restric¬ 
tion  will  be  violated  for  short  distances  at  the  lower  end  of  the  whistler 
spectrum,  since  the  wavelength  at  2  kc  is  150  km.  The  results  obtained, 
however,  can  still  be  significant  in  terms  of  gross  behavior  of  the 
whistler  signals.  For  this  reason  the  calculations  will  be  carried  in¬ 
to  this  range  for  the  purpose  of  comparison  with  the  available  experi¬ 
mental  data.  Over  most  of  the  ranges  of  frequency  and  distance  of  in¬ 
terest,  however,  the  wave  may  be  considered  plane  and  the  equations 
given  may  be  applied  with  very  little  error. 

2.  The  quasi -longitudinal  (<j,l)  approximation  is  valid  for  this 
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problem.  According  to  Budden,  the  ql  approximation  is  valid  where 


«  1 


(2.1) 


This  may  be  expressed  as 
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where  9  is  the  angle  bet wean  tbs  wave  notasl  and  the  static  magnetic 
field,  and 

yT  "  yR  6ia  6  ' 

yL  =  ya  COS  0 


f 

•  r* 

a  IH 


v 

z  =  2itf 

«  plasma  frequency  (cycles/sec) 


f  *  ti. e-  «  gyrofrequency  (cycle/ sec) 

ri  dTtffl 

N  »  electron  density  (number/m^) 
e  =  electronic  charge  (coulooibs) 
m  *  electronic  mass  (kg) 

H  =  strength  of  static  magnetic  field  (smp/aeter) 
f  =  wave  frequency  (cpa) 
v  «=  collision  frequency  (colliaions/aec) 
kq  =  peraitivity  of  free  apace  (farad/meter) 
li'=  permeability  of  free  space  (hy /meter) 

Rationalized  aka  units  will  he  used  except  where  otherwise  noted. 
This  restriction  will  he  most  critical  at  the  higher  frequencies  because 
of  the  r  in  the  denominator  of  the  second  tern.  When  f  ■  20  kc  (the 
highest  frequency  of  interest  here)  and  for  a  typical  suaner-nlgjit -model 
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ionosphere  (see  page  37),  the  result  is: 

|  tan  0  sin  ©|  «  8  (2.3) 

When  %  error  is  chosen,  the  result  is; 

sin  0  tan  9  .  =  ’O.b 

9  =  35° 

The  angle  of  refraction  in  the  ionosphere  id  complex.  The  best 
way  to  interpret  this  fact  is  to  use  the  angle  of  the  were  normal  -ty,  de¬ 
fined  on  p.,166.  The  maximum  angle  of  refraction  vs  frequency  is 

g* 

plotted  in  Fig.  F-l.  The  value  of  for  the  summer -ci ght -model  iono¬ 
sphere  is  approximately  28  deg.  The  angle  between  the  earth's  magnetic- 

field  and  the  vertical,  measured  in  the  same  sense  as  ifr.  Is  27.5  deg 

c  er 

from  the  vertical  at  the  latitude  of  Stanford  (4b  geomagnetic).  If 
propagation  is  toward  the  magnetic  pole  in  the  magnetic  meridian  the  ql 
approximation  is  valid,  since  the  wave  normal  will  make  small  angles  with 
the  earth's  field.  This  statement  is  also  essentially  true  for  directions 
on  either  side  of  the  meridian  where  the  rays  make  small  angles  with  the 
field.  The  wave  normal  would  have  to  be  almost  perpendicular  to  the  field 
for  propagation  toward  the  equator.  The  ql  approximation  would  be  very 
poor  for  this  direction.  The  region  of  validity  increases  rapidly  as 
the  frequency  is  decreased,  however,  and  the  approximation  is  valid  for 
all  directions  up  to  90°  on  either  aide  of  the  direction  of  the  magnetic 
pole,  at  the  low  end  of  the  whistler  range. 

The  exact  conditions  for  validity  are  very  sensitive  to  plasma  fre¬ 
quency  and  collisional  frequency..  If  the  plasma  frequency  were  doubled, 
which  is  not  unreasonable  for.  night-time  conditions,  the  ql  approximation 
becomes  valid  for  essentially  all  directions.  It  must  be  recalled  that 
the  parameters  used  in  the  summer-night  model  were  determined  for  this 
same  approximation,  to  match  experimental  ground-to-ground  propagation 
data  for  all  directions  of  propagation  In  the  15-  to  20 -kc  range.  There¬ 
fore  the  range  of  application  is  probably  greater  in  the  upper  frequency 
range  than  the  Inequality  suggests. 
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It  must  also  be  considered  that  the  longitudinal  component  of  the 

Qrrofrequency  y^  is  actually  complex  in  the  case  of  the  lossy  ionosphere, 

and  its  magnitude  varies  vith  the  angle  of  refraction.  Another  result 

of  the  ql  assumption  is  that  yL  is  represented  by  a  real  value  vith  an 

effective  magnitude  determined  to  fit  experimental  results  best.  When 

the  ql  assumption  is  not  make,  a  complete  solution  ie  necessary,  in  the 
17 

manner  of  Yabroff,  for  the  case  for  which  the  refractive  index  is  a 
function  of  the  wave  normal  angle.  This  requirement  would  so  complicate 
the  calculations  that  they  would  be  no  longer  practical.  The  justification 
for  this  assumption  is,  first,  that  previous  results^' of  calculations 
using  these  assumptions  for  ground-to-ground  propagation  problems  are 
consistent  with  experimental  data;  second,  the  results  obtained  in  this 
report  are  consistent  with  the  experimental  data  on  whistlers. 

3.  The  third  assumption  is  that  the  lower  edge  of  the  ionosphere 
may  be  adequately  represented  by  an  equivalent  homogeneous  model  with  a 
sharp  lower  boundary.  This  assumption  is  very  difficult  to  treat  quanti¬ 
tatively.  It  is  generally  accepted  that  a  sharp  boundary  approximation 
is  reasonable  if  there  is  a  large  change  in  the  properties  of  the  medium 
within  the  distance  of  a  fraction  of  a  wavelength.  There  is  a  sharp 
gradient  of  electron  density  at  the  lower  edge  of  the  E  layer,  so  that 
this  condition  is  fulfilled  over  most  of  the  frequency  range  of  interest. 

The  homogeneous  model  may  be  used  if  the  changes  in  the  medium  above 
the  sharp  boundary  are  slow  enough  that  there  is  not  a  significant  change- 
in  the  properties  of  the  medium  within  a  distance  of  the  order  of  a 
wavelength.  A  superposition  of  a  number  of  these  slabs  of  different' 
properties  may  be  used  to  describe  a  rapidly  varying  medium  in  the  limit 
as  the  number  of  slabs  becomes  very  large  and  their  thickness  very  small. 

In  this  report,  one  homogeneous  medium  is  used.  The  medium,  which 
approximates  the  ionosphere,  must  be  sufficiently  uniform  that  the  wave 
is  attenuated  to  the  point  of  being  negligible  before  it  encounters  a 
significant  change.  This  uniformity  would  thus  have  to  exist  over  a 
distance  of  at  least  the  order  of  a  wavelength.  The  wavelength  in  the 
medium  decreases  approximately  inversely  as  the  square  root  of  the  fre¬ 
quency.  The  actual  distance  over  which  the  medium  must  not  change  signifi¬ 
cantly  thus  depends  on  the  frequency,  the  refractive  index,  and  the  angle 
of  refraction  in  the  medium. 
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An  tottnl  quantitative  statement  of  the  criteria  for  the  validity 
of  this  assumption  has  not  been  made  and  cannot  be  made  with  the  present 
state  of  knowledge  of  propagation  in  a  conductive,  inhomogeneous  medium. 
The  only  Justification  for  tMs  assumption  is  that  the  results  are  con¬ 
sistent  with  the  available  experimental  data  on  ground-to-ground  trans- 
mission  and  whistlers. 

4.  The  fourth  assumption  is  that  the  ionosphere  contains  plane- 
wave  energy,  with  all  wave  normals  contained  within  a  cone  whose  axis  is 
normal  to  the  boundary  and  whose  half  angle  is  equal  to  the  "critical" 
angle  for  total  internal  reflection.  Furthermore,  it  is  assumed  that 
the  transmitted  waves  may  be  considered  to  radiate  from  a  point  source. 

This  fourth  assumption  probably  does  not  completely  describe  the 

actual  conditions,  since  the  ionosphere  is  in  general  a  lossy  medium  and 

the  range  of  wave  normals  actually  available  to  excite  the  waveguide  may 

be  restricted  by  the  trapping  characteristics  of  a  duct.  When  a  plane 

wave  in  a  dielectric  medium  is  incident  on  a  lossy  medium,  the  wave  in 

18 

the  lossy  medium  is  what  Stratton  calls  an  "inhomogeneous  plane  wave." 

The  assumption  of  plane  waves  is  the  only  practical  assumption  in 

* 

the  present  state  of  knowledge.  The  Justification  must  again  be  the 
fact  that  the  results  are  reasonable  in  view  of  experimental  data. 

The  field  at  a  receiving  antenna  will  then  consist  of  the  sum  of  all 
the  rays  of  significant  amplitude  that  will  propagate  between  the  exit 
point  and  the  antenna.  The  assumption  that  all  necessary  wave  normals 
are  present  results  in  excitation  for  all  these  components. 

5.  It  is  assumed  that  all  the  energy  in  a  set  of  sky-wave  components 
is  of  the  same  frequency.  Actually,  the1  various  rays  contributing  to  a 
whistler  signal  at  a  given  time  are  of  slightly  different  frequency  be¬ 
cause  of  the  frequency  variation  of  the  whistler  and  the  difference  in 
propagation  times  of  the  rays. 

The  maximum  path  difference  between  two  adjacent -order  rays,  which 
occurs  for  vertical  incidence,  is  of  the  order  of  200  km  (one  round  trip 
between  earth  and  ionosphere).  A  typical  dispersion  of  100  sec^  yields 


* 

A  brief  discussion  of  critical-angle  phenomena  and  waves  in  lossy 
media  is  given  on  page  1 $6. 
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a  frequency  difference  of  about  13  cps  at  1C  ktf.  The  frequency  difference 
decreases  as  the  distance  is  increased  and  it  about  h  cps  at  a  distance 
of  300  km.  A  comparison  of  calculated  and  experimental  data  shews  that 
this  phenomenon  is  not  generally  important*  Further  discusalon  of  beating 
is  given  on  page  128. 

This  single -frequency  assumption  means  essentially  that  the  calculations 
may  be  performed  on  a  steady-state  basis,  and  it  is  not  in  general  neces¬ 
sary  to  solve  the  transient  problem.  The  transient  problem  is  extremely 
complicated,  and  attempts  to  solve  it  do  not  appear  to  be  warranted  until 
the  results  of  steady -state  solutions  have  been  more  completely  examined. 

C.  PATH  GEOMETRY 

It  is  now  possible  to  formulate  a  system  of  equations  for  calculating 

the  total  field  for  the  whistler-mode  signal  which  arrives  at  a  receiving 

20 

antenna.  The  equations  for  the  curved-earth  geometry  will  be  used  in 
order  to  determine  the  geometrical  characteristics  of  the  ray  path  it¬ 
self.  It  is  assumed  that  the  earth  is  a  perfect  sphere  of  radius 
6,367.39  km.10 

Refer  to  Figs.  2.1  and  2.2  for  an  example  of  the  ray  geometry  and 
for  the  notation.  The  angle  9^,  calculated  from  the  path  geometry,  is 
determined  by  the  virtual  height  of  the.  ionosphere  h',  the  ground  dis¬ 
tance  D  to  the  receiver  (mfcs  units  will  be  used,  except  for  h*  and  D  in 
km),  and  the  number  of  ionosphere  reflections  n. 

The  ionosphere  angles  of  refraction  0Q  and  0 ^  are  then  calculated 
from  Snell's  Law; 

u  sin  9  =  u  sin  0  »  sin  0 

e  “0  o  n 

where  |i  =  the  complex  refractive  index  of  the  ordinary  wave  (see  p.  16) 
o 

=  the  complex  refractive  index  of  the  extraordinary  wave. 
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D.  IONOSPHERE  MODEL  AND  ITS  REFRACTIVE  INDEX 


The  next  step  is  to  define  the  ionosphere  model  and  its  refractive 

index.  It  will  be  recalled  that  we  assumed  for  purposes  of  this  problem 

that  the  ionosphere  may  be  represented  by  a  homogeneous  Medium  with  a 

sharp  lower  boundary  and  a  vertical  magnetic  field. 

With  these  assumptions,  the  refractive  index  of  the  medium  has  been 
12 

expressed  by  Budden  as: 


where 


2  +JT 

U  *  1  -  j  (oy'oi)  £ 

=  1  +  CD CD  sin  T  -  jCD CD  COS  T 
tan  t  =  cd  /v 

±J 


CD 


CD  = 

r 


#7 


2 

“l 


x2  Ne 

o  /cm 
o 


(2.4) 
(2.4a) 

(2.5) 

(2.6) 


(2.7) 


H’He 

=  ST 


(2.8) 


N  =  number  of  electrons  per  meter 
e  =  charge  of  the  electron  in  coulombs 
m  =  mass  of  the  electron  in  kg 

K  «=  permittivity  of  free  space  in  farad/meter 

°  -1 
v  =  colli sional  frequency  in  sec 

2 

H  =  the  earth's  magnetic  field  in  amp/meter 
2 

This  expression  for  p.  is  double  valued,  as  is  shown  by  the  +  sign. 
The  whistler  mode  corresponds  to  the  +  sign  or  the  ordinary  (o)  mode. 


The  term  ordinary  is  used  here  to  conform  to  the  reference  from 
which  this  material  was  obtained.  Modern  practice, 36  however,  is  to 
label  the  whistler  mode  as  the  extraordinary  mode. 


(2.9) 


The  refrtdtlve  index  4  is  then  given  by 

0 

2 

=  1  +  (ay'aa)  (sin  t  -  j  cos  t) 

The  coordinate  system  and  field  components  are  shown  in  Fig.  2.3-  The 
polarization  for  the  case  of  a  downecmiog  wave  in  the  northern  hemisphere 
is,  then^ 


=  -  J 

^.0 


(2.10) 


where  the  II  and  i  indicate  that  the  electric  field  is  parallel  or  per¬ 
pendicular  to  the  plane  of  incidence.  Similarily,  for  the  extraordinary 
( e )  mode^ 

2 

H  =  1  -  (oj J<£>)  (sin  t  +  j  cos  t)  (2.11) 


Q 


e 


+  J 


The  refractive -index  expressions  are  thus  stated  in  a  convenient  form  as 
a  function  of  the  normalized  parameters  and  t. 


E.  TRANSMISSION  COEFFICIENTS 


The  expressions  for  the  transmitted  waves  may  now  he  written,  as¬ 
suming  that  the  incident  wave  is  plane  and  has  the  polarization  of  the 
"ordinary"  or  whistler  component.  Since  the  ql  approximation  is  valid 
this  wave  will  he  circularly  polarized.  The  polarization  at  the  boundary 
will  thus  he  such  that  it  will  excite  a  significant  reflected  wave  only 
in  the  same  mode,  since  the  e  mode  requires  the  oppositely  rotating 
circular  polarization.  The  e  mode  will  therefore  he  negligible  under 
ql  conditions.  The  wave  in  the  ionosphere  has  been  given  as: 


* 


A  discussion  of  the  point  is  given-  in  Appendix  E. 


*l|o  ”  *11  i 

Elo  -  Vh  •  \  1 

where  o  signifies  the  ordinary  Mode  of  propagation  and  the  It  and  i  indi¬ 
cate  whether  the  electric  field  is  parallel  to  or  perpendicular  to  the 
plane  of  incidence.  The  coordinate  system  is  defined  in  Fig.  2.3«  The 
x^--  x^  plane  is  the  lower  boundary  of  the  model  ionosphere,  and  the 
Xi~  plane  is  the  plane  of  incidence.  The  transmitted  components  be¬ 

low  the  ionosphere  may  then  be  expressed  by: 

V  - I.  (2-15) 


(2.13) 

(2.1*0 


Elt  =  Eli  Ti 


(2.16) 


where  T(|  and  T  are  the  transmission  coefficients  '.rich  are  to  be  deter¬ 
mined.  This  is  done  by  matching  the  four  field  components  parallel  to 
the' boundary .  These  components  are: 


C  E  -.  \  C  =  E  C 
i  111  Hr  r  It  t 


Eii+Elr  =Eit 


H..i  ci  -  ^lr  cr  *  V  Ct 


Hii  +  Hlr  =  Hit 


(2.17) 

(2.18) 

(2.19) 

(2.20) 


where  C  represents  cos  6,  and  the  subscripts  i,  r  and  t  refer  to  the 
incident,  reflected,  and  transmitted  waves,  respectively.  In  this  case 

“1  "  “r  ’  “o 

and 
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Tike  ratios  of 


where  the  o  refers  to  the  ordinary  node  ad T  propagation, 
the  components  of  E  to  H  may  now  be  expressed  by: 


E  u  Si 

- 1^7*  '  =  rj  /n 

TT  TT  W  V 


Hi  H„ 


where 


(2.21) 


H'  =  the  permeability  of  free  space, 

^  =  the  permittivity  of  the  medium. 

We  now  substitute  Eqs.  '2.21)  in  (2.19)  and  (2.20)  with  the  proper 
subscripts  and  divide  by  r}Q,  obtaining: 


Mi .  C .  E  . .  -  n  C  E.  *  C  E 
i  i  ii  r  rir  t  it 


(2.19a) 

+  “rV  *  El>t  (2>20a) 

Wow,  multiplying  Eq.  (2.17)  by  nr  and  Eq.  (2. 20)  by  Cr,  and  adding,  yields: 

“ill'rV'l'r'-AV'At 


T  .  *Jt  .  jVVl^fr 

“  E»i  “rCt  +  °r 


(2.22) 


Next,  multiply  Eq.  (2.18)  by  and  add  to  Eq.  (2.19a),  yielding: 


EiJ>rCr  +  -  Wr  +  °t> 


Tl* 


E  n  C  +  n  C 
it  r  r  ii 


(2.23) 


E 


'ii 


^rCr  +  Ct 


-19 


Now,  for  the  whistler  case  |i,  =  u  ■  u  and  C.  «  C  •  C  ,  by  Snell' •  ll». 

1  r  O  1  r  o 

Therefore,  T(|  and  may  be  expressed  as} 


C  (i  +  H  c 
o  o  o  o 

n  c.  +  c 

o  t  o 


2|i  C 

O  O 

n  c+  +  c 

o  t  o 


T 


2  p.  C 
o  0 

ji  C  +  C 
oo  t 


F.  RAY  ADDITIONS 


(2.25) 


The  resultant  transmitted  wave  expressed  by  the  preceding  equations 

is  then  reflected  n  times  from  the  earth  and  the  ionosphere  before 

arriving  at  the  receiver.  The  reflections  from  the  ionosphere  are  con- 
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veniently  described  by  the  reflection  coefficients  derived  by  Budden 
13 

and  used  by1  Wait.  There  are  four  of  these  coefficients,  designated 
II R  ||  ,  j^Rj^,  ||  Rj.,  j_  R|| .  The  first  subscript  indicates  whether  the  incident 
electric  field  is  parallel  or  perpendicular  to  the  plane  of  incidence, 
and  the  second  similarly  indicates  the  position  of  the  reflected  electric 
field.  The  existence  of  the  cross -coupling  terms  jjRjj  and  jR||  is  due  to 
the  effect  of  the  earth's  magnetic  field,  which  causes  the  characteristic 
waves  in  the  ionosphere  at  vlf  to  be  elliptically  polarized  in  general. 
Thus,  if  the  incident  wave  is  linearly  polarized,  the  reflected  wave 
must  contain  some  of  the  cross  polarized  component  in  order  to  match  the 
boundary  conditions.  Some  of  the  energy  in  a  given  linear  polarization 
is  therefore  converted  into  a  reflected  component  which  has  a  polarization 
perpendicular  to  that  of  the  incident  wave. 

The  reflections  from  the  earth  are  expressed  in  terms  of  the  Fresnel 
reflection  coefficients  R(|  and  R^  .  These  are  functions  of  the  dielectric 
constant  and  conductivity  of  the  ground,  the  wave  frequency,  and  of  the 
takeoff  and/or  arrival  angle  0^,  (hereafter  referred  to  as  the  takeoff - 
arrival  angle)  measured  from  the  horizontal.  The  subscripts  on  the 
coefficients  and  indicate  the  polarization  of  both  the  incident 
and  the  reflected  eleptric- field  at  the  ground. 
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The  expression  for  the  total  field  for  each  ray  may  nov  ha.  given, 
using  these  six  reflection  coefficients  and  the  transmission  coefficients 
derived  in  Section  II -E.  Equations  (2.15)  and  (2.1$)  express  the  trans¬ 
mitted  fields  as: 


E»t  ■  EmT« 

(2.15) 

Eit  ■  Vi 

(2.16) 

The  addition  of  components  may  now  be  expressed  as: 

kE.l  " 

k-1  EllRlillRll  +  k-lEiRllRl 

(2.26) 

ii 

k-iEiiRiinRi  +  k-iBA  A 

(2.27) 

Where 

oeh  *  Eiit 

(2.28) 

0Ei  -  Elt 

(2.29) 

The  subscript  k  is  a  dummy  index  which  indicates  that  the  calculation  is 
performed  repetitively  n  times  (i.e.,  until  k  *=  n).  The  quantities  kEM 
and  kEi  for  k  =  n,  represent  the  fields  of  the  n-hop  ray  uncorrected 
for  path  length.  The  fields  are  then  corrected  for  path  length  by: 


,.(n) 


E 


E  —  cos 
k  .1  Rn 


(n) 


A 

Ei  r  cos  A +  J  sin  A 


when  k  =  n 


when  k  =  n 


(2.30) 

(2.31) 


where 

*ir>  - 

n) 


a  = 
n 

R  = 

&■ 

R 


the  parallel  polarized  component  of  the  nth -order  ray 

the  perpendicular  polarized  component  of  the  nth-order  ray 

the  convergence  coefficient  of  the  nth-order  ray  • 

the  propagation  distance  of  the  nth-order  ray 

the  phase  of  the  nth-order  ray 

and  are  defined  in  Appendix  A. 
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This  method  of  expressing  the  general  n-hop  component  is  analogous  t# 
finding  the  output  of  a  multi-stage  network  by  splitting  it  into  a  number 
of  transfer  functions,  calculating  the  output  of  each  section  in  turn, 
and  using  the  output  of  that  section  as  the  input  to  the  next. 

The  convergence  coefficient  an  is  a  factor  derived  to  account  for 
the  effect  of  the  reflections  from  the  concave  ionosphere  surface  and 
the  convex  earth  surface  In  the  ground-to-ground  case,  the  coefficient 
expresses  the  net  convergence  due  to  n  reflections  from  the  ionosphere 
and  n  -  1  from  the  earth.  The  net  convergence  is  generally  very  small 
in  this  case,  except  near  the  critical  distance  where  the  expressions 
become  invalid  and  a  higher-order  approximation  must  be  used.  In  the 
present  case, .the  number  of  reflections  from  the  earth  and  the  ionosphere 
is  equal  and  the  two  effects  would  tend  to  cancel  each  other.  It  is 
therefore  reasonable  to  assume  a  value  of  unity  for  this  coefficient. 

This  will  lead  to  some  error  in  the  calculation  of  a  ray  as  its  critical 
distance  is  approached.  The  error  due  to  this  factor  is  opposite  in 
effect  to  that  of  neglecting  the  diffraction  factor  near  the  critical 
distance,  and  the  resulting  error  is  small. 

G.  CONVERGENCE  OF  THE  SERIES  OF  RAYS 

We  have  shown  how  the  field  is  calculated  as  the  sum  of  n  sky-wave 
components.  Each  succeeding  component  is  of  lower  amplitude  than  the 
previous  one,  due  to  the  l/R  amplitude  factor  and  the  multiplication  of 
further  reflection  coefficients,  which  are  always  less  than  unity  in 
magnitude.  Thus,  a  series  is  formed  with  terms  of  successively  decreas¬ 
ing  amplitude.  It  is  now  necessary  to  determine  the  point  at  which  the 
series  should  be  terminated.  The  method  used  is  as  follows:  The  six 
components  for  n  =  0  through  5  are  always  calculated.  The  field  from 
each  component  is  resolved  into  its  three  axial  components  (see  p.  24) 
and  each  axial  component  is  added  to  the  sum  of  that  axial  component 
built  up  by  previous  -rays.  The  amplitude  of  the  last  ray  is  then  com¬ 
pared  with  the  amplitude  of  the  sum  for  each  of  the  three  axial  components 
in  turn.  The  calculation  is  terminated  when  the  amplitude  of  the. last 
ray  is  less  than  a  specified  percentage  of  the  smallest  of  the  three  sums. 
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this  percentage  was  experimentally  studied  tad  the  value  of  1  per  a  ant 
was  adopted.  The  error  In  neglecting  further  components  Is  less  than 
1  per  cent  in  amplitude  and  about  0*1  degree  In  the  angles  calculated 

for  all  cases  studied. 

H.  CRITICAL  DISTANCE  EFFECTS 

The  critical  distance  is  that  which  yields  a  value  of  zero  for  6  , 

n 

the  takeoff -arrival  angle  at  the  earth's  surface.  When  the  critical 
distance  for  a  given  ray  is  exceeded,  a  diffracted  field  exists  at  the 
antenna.  The  expressions  for  the  diffracted  field  are  very  complicated. 
Recently,  Wait  *  has  published  approximate  solutions  to  this  problem. 
In  order  to  preserve  a  reasonable  simplicity v in  the  over-all  calculation, 
it  was  specified  that  the  field  for  a  given  ray  is  not  calculated  if  f*n 
is  less  than  zero.  The  effect  of  the  refracted  field  in  the  "shadow 
zone"  is  thus  ignored.  This  limits  the  accuracy  of  the  calculations  for 
distances  just  beyond  the  critical  distance.  There  is  an  additional 
effect  on  waves  which  have  a  low  takeoff-arrival  angle  which  is  caused 
by  the  presence  of  the  earth's  surface.  This  effect  was  also  ignored  in 
order  to  maintain  simplicity  of  the  calculations.  This  effect  is  similar 
to  that  dealt  with  by  Wait  in  his  consideration  of  the  "cutback"  factor, 
which  expresses  the  modification  of  the  pattern  of  a  radiating  antenna 
at- angles  close  to  the  earth's  surface.  Approximate  solutions  to  similar 
problems  also  exist,  but  these  were  also  discarded  in  order  to  preserve 
the  essential  simplicity  of  the  calculations.  The  error  due  to  these 
factors  is  believed  to  be  small.  No  noticeable  discontinuities  appear  at 
or  near  the  critical  distance,  so  that  it  appears  .reasonable  that  the 
results  are  useful  near  these  critical  distances. 

I.  TOTAL  FIELD  AND  ANTENNA  RESPONSE 

1.  GENERAL  CONSIDERATIONS.  We  have  expressed  the  electric  field 
due  to  all  significant  sky  wave  components.  Each  of  these  component* 
arrives  at  an  angle  g  ,  determined  by  its  path  geometry.  The  problem 
now  is  to  determine  the  total  field  and  the  response  of  the  antenna 
system  to  it. 
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Equations  will  be  derived  for  tUm  total  electric  field  at  the  antenna, 
which  is  assumed  to  he  small  compared  with  the  wavelength.  It  i s  im¬ 
portant  to  realize  at  this  poiat  that  the  field  which  is  effective  in 
inducing  a  voltage  in  an  antenna  will  not  necessarily  he  the  same  as  the 
total  electric  field  present  at  the  antenna.  The  type  of  antenna  used 
determines  the  exact  form  of  the  field  which  will  affect  it.  These 
differences  in  response  for  different  types  of  antennas  must  he  carefully 
considered  when  the  form  of  the  field  to  he  calculated  is  determined. 

2.  TOTAL  FIELD.  First,  we  consider  the  problem  of  the  behavior 
of  the  total  electric  field  as  a  function  of  time  in  an  area  small  with 
respect  to  a  wavelength.  It  is  commonly  expected  that  the  total  field 
vector  will  describe  a  complicated,  non-plane  figure  in  space,  even  in 
the  single -frequency  case.  This  is  not  true,  however.  It  can  be  shown 
that  the  figure  is  in  general  an  ellipse  contained  in  a  plane  whose 
orientation  is  readily  described  in  terms  of  the  component  rays.  A 
series  of  rays  may  thus  be  described  by  a  single  elliptically  polarized 
field.  We  will  therefore  consider  this  problem  in  detail. 

It  is  easily  seen  that  the  series  of  rays  arriving  at  different 
angles  and  each  containing  a  parallel  and  a  perpendicular  component 
of  electric  field,  will  have  components  in  all  three  coordinate  directions. 
Therefore,  we  first  resolve  the  waves  into  these  components: 

=  -  Z  E|,n^  sin  Pn  (2.32) 

n 

Ex2  =  2  ^  (2‘33) 

n 

E.  =  -  Z  E||n‘>  cos.Pn  (2.3*0 

J  n 

The  coordinates  are  the  same  as  in  Fig.  2.3.  The  three  resultant  com¬ 
ponents  are  independent  of  one  another  in  phase  and  amplitude. 

Now,  the  position  of  the  total  field  vector  as  a  function  of  time 
must  be  determined.  It  has  been  previously  assumed  that  the  waves  are 
all  sinusoidal  and  of  the  same  frequency.  Let  the  origin  of  the  coordin¬ 
ates  be  the'  base  of  the  vector.  The  tip  of  the  vector  will  then  describe 
a  pattern  in  three  dimensions,  with  its  displacement  in  each  of  the 


2k  -  ’ 


coordinate  directions  given  by  Eqs.  (2.32)  to  (2.3*0.  Sine#  the  Wre® 
are  sinusoidal,  each  of  these  displacements  must  be  symmetrical  about 
the  origin,  which  means  that  the  figure  passes  through  the  origin.  Also, 
each  component  has  the  same  number  of  maxima  and  minima  per  qyele,  since 
the  frequency  is  the  same.  These  restrictions  determine  that  the  figure 
is  a  plane,  whose  orientation  is  determined  only  by  the  relative  phases 
and  amplitudes  of  the  components.  (This  same  conclusion  is  reached 
mathematically  on  page  26.) 

Since  the  three  axial  components  of  field  strength  have  independent 

phases  and  amplitudes,  the  figure  which  the  field  vector  describes  will, 

in  the  general  case,  be  an  ellipse.  The  simplest  way  to  describe  the 

resultant  is  to  determine  the  orientation  of  the  plane  of  this  ellipse, 

and  then  rotate  coordinates  so  that  the  ellipse. can  be  described  by  two 

coordinates  and  two  angles  of  rotation. 

Consider  the  plane  of  the  ellipse.  This  plane  is  most  easily 

described  by  the  direction  cosines  of  its  normal.  Let  these  direction 

cosines  be  i,  m,  and  n,  in  the  x. ,  x  ,  and  x  directions,  respectively. 

1  23^ 

The  general  equation  of  a  plane  is  then: 

£x1  +  mx2  +  nx^  =  p  (2.35) 

where  p  is  a  constant  equal  to  the  length  of  the  normal. 

The  coordinate  values  then  may  be  replaced  by  the  actual  electric 
field  strength  values,  which  are  sinusoidal  functions  of  time.  Let: 

(2.36) 

(2.37) . 

(2.38) 

(2.39) 


x,  =  E  sin  cut  =  E  sin  cut 
1  xi  xi 


x  =  E  sin  (cot  +  8)  =  E  sin  cot 
2  x^  v  y '  x- 


x  =  E  sin  (cut  +  0  )  =  E  sin  cut 
3  x  ■  '  z'  x3 

where  E  ,  E  ,  and  E  are  complex  numbers. 

X1  X2  x3 

Then  (2.35)  becomes: 

sin  cut  (IE  +  mE  +  nE  )  *  P 
•  •  ■  X1  X2  x3  ' 
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This  equation  can  be  satisfied  only  by  p  -  0  for  all  value*  of  t.  This 
means  that  the  plane  must  pass  through  the  origin,  as  ve  found  previously. 
Nov  let: 


Ry  ■  —  ■  Ry  (cos  0y  +  j  sin  0  ) 
*1 


R  =  -*  =  R  (cos  0„  +  J  sin  0  ) 
z  z  z  z 

X1 


Now,  from  (2.39)  for  a  solution  to  exist: 


(.2.  ho) 

(2A1) 


and 


E  (i  +  mR  cos  8  +  nR  cos  8  )  -  0 

V  '  T  t  mr  mm  •  ' 


E  (mR  sin  6  +  nR  sin  9  )  =  0 

x.  y  .  y  z  z 


2  2  2 

1  +  m  +  n  =1 


(2.42) 

(2.43) 


(2.44) 

from  the  properties  of  direction  cosines.  Now,  multiply  Eq.  (2.42)  by 
sin  6  and  Bj,  (2.43)  by  cos  6  ,  yielding 

y  y  ^ 


y 


+  mR 

sin 

0 

cos  0 

+  nR 

sin 

0 

cos  0  =  0 

y 

y 

y 

z 

y 

z 

mR 

sin 

0 

cos  0 

+  nR 

cos 

0 

sin  0  =  0 

y 

y 

y 

z 

y 

z 

(2.42a) 

(2.43a) 


Then,  subtract  Eq.  (2.43a)  from  Eq.  (2.42a),  which  results  in: 


1  sin  9  +  nR  (sin  0  cos  0  -  cos  0  sin  0  )  =  0  (2.45) 

y  2  y  2  y  Z  ' 

This  gives  the  relation  between  n  and  1.  Next,  multiply  Eq.  (2.42)  by 

sin  0  and  Eq.  (2.43)  by  cos  0  ,  yielding: 
z  z 


+  mR 

cos  0 

sin  0  +  nR  sin 

0 

cos  0  *  0 

(2.42b) 

y 

y 

z  z 

z 

z 

mR 

sin  0 

cos  0  +  nR  sin 

0 

cos  0=0 

(2.43b) 

y 

y 

z  z 

z 

z 

Next,  subtract  Eq.  (2.43b)  from  Eq.  (2.42b),  which  yields: 
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I  «in  €  *  mfi  (eofl  6  sin  9  -  sin  9  soft  9  )  m  0  (2 .46) 

•» .  z  y  y  z  y  z 

Jfcft  next  step  is  to  substitute  Eq.  (2.45)  and  Eq.  (2.46)  into  Eq.  (2.44), 

2 

which  yields  an  equation  for  l  alone: 


.  J2.  J  2  „  _2  .  2  . 

~  1  +  R  sin  6  +  R  sin  6 

i2  - i - 1 - 1 - 1 - 

R  R  (sin  Q  cos  0  -  cos  0  sin  0 

y  z  y  z  y  z 

This  equation,  when  rearranged,  yields: 


(2.47) 


op  2 

R  R  (sin  0  cos  0  -  cos  0  sin  0  ) 

.  y  z _ _ z _ jr_  z 


R2  R2  (sin  0  cos  0  -  cos  0  sin  0  )2  +  R2  sin2  0  +  R2  sin2  0  . 

yz  y  z  y  z  y  yz  z 


Then  m  and  n  are  found  from  Eq.  (2.45)  and  Eq.  (2.46),  which,  when 
rearranged,  yield: 


(2.48) 


t  sin  0 


m  = 


R  (sin  0  cos  0.  -  cos  0  sin  0  ) 

y  y  z  y  z 


(2.49) 


and 


i  sin  0 


n  = 


R  (sin  0  cos  0  -  cos  0  sin  0  ) 

z  y  z  y  z 


(2.50) 


Now,  the  elevation  and  azimuth  angles  of  the  normal,  B  and  0  must 
be  found.  Refer  to  Fig.  2.4  for  an  illustration  of  the  nomenclature. 
The, direction  of  propagation  of  the  ray  is  the  positive  x^"  direction. 
The  line  a  is  the  projection  of  the  direction  of  propagation  on  the 

x,  --  x  .  plane.  The  angle  0  is  the  angle  of  rotation  in  that  plane. 

12.  e 

The  angle  -  is  the  corresponding  depression  angle.  It  is  readily 
apparent  that : 
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tan  0 

e 


ela  0 

e  _  ft 

*  a 

cos  0 

e 


a  _  o 
7*1 


a  =  cos  Pg  =  Vl  - 

1  -  cos^  P  =  sin^  p  =  1  -  (l  -  n^) 
e  e 


and 


sin  p  =  +  n 
Ke  - 


+  n 

tan  P  =  - 

e 


(2.51) 

(2.52) 


The  next  step  is  to  rotate  the  coordinates,  first  about  x^  through 

an  angle  0  ,  which  results  in  axes  x  ' ,  Then  we  rotate  the  coordinates 
e  n 

about  x'  through  an  angle  p  ,  which  results  in  coordinates  x".  These 
2  e  n 

operations  may  be  performed  as  follows:  first,  for  each  x  substitute 

13  n 

its  proper  replacement  in  the  primed  system. 


x|  cos  0e-  Xg  sin  0g=  Ex 

(2.53) 

xj,  sin  0e+  ^  cos  0e=  iT 

(2.5^) 

(2.55) 

which  expresses  the  first  rotation  of  0g  degrees  about  x^. 

Then  replace  each  of  the  primed  coordinates  by  its  proper 

in  the  double -primed  system: 

substitute 

(x£  cos  Pe-  x^  sinPe)  cos  0g-  x^  sin  0g=  IT 

(2.56) 

(x^  cos  Pe-  x^  sin  Pg)  sin  0Q+  x^  ,cos  0g=  IT 

(2.57) 

x”  sin  pe+  x”  cos  pg=  E^ 

(2.58) 
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Vfclch  expresses  the  rotation  of  Pe  degrees  about  x^.  Now,  rearranging 
and  collecting  tej»e  yields; 


x''  cos  P  cos  0  -  x"  sin  0  -  x"  sin  p  cos  0  =  E~  (2.59) 

1  e  e  2  e  3  e  e  x,  v 

x"  cos  P  sin  0  +  x"  cos  0  -  x"  sin  p  sin  0  =  E~  (2.60) 

1  e  e  2  e  3  e  e  x2 

x"  sin  p  +  x"  cos  p  =  E~  (2.6l) 

1  e  3  e  x^ 

These  equations  are  now  in  a  form  suitable  for  the  solution  of  the 
system  by  determinants. 

First,  the  denominator  of  the  system  is  given  by 


D  = 


cos  P  cos  0  -  sin  0 

6  e  6 


cos  p  sin  0 
e  e 

sin  p 


cos  0 


-  sin  p  cos  0 

e  e 

-  sin  p  sin  0 

e  e 

COS  P 


(2.62) 


2  2  2 
=  sin  0  (cos  P  sin  0  +  sin  B  sin  0  )  +  cos  0  (cos  p  cos  0 
e  e  6  6  e  e  0 .  e 


2  2  2 

+  sin  P  cos  0  )  =  sin  0  +  cos  0  =1 

e  e  e  .  e 


(2.63) 


Then  xl^  is  given  by 


x"  = 
X1 


E 


-  sin  0 


cos  0 


-  sin  p  cos  0 
e  e 


sin  p  sin  0 
Ke  e 


cos  P 


When  normalized,  this  yields: 
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E 


-  sin  0 

e 

cos  0 


-  sla  0  cos  0 

e  e 

-  sin  0  sin  0 

e 

cos  0 


(2.6k) 


x?  =  E  [sin  0  (R  cos  0  +  R  sin  0  sin  0  )  +  cos  0  (cos  0 
1  x^  e  y  e  z  Ke  e  e  e 

+  R  sin  0  cos  0  )]  =  E  (R  sin  0  cos  0  +  R~  sin  0  sin2  0 

z  e  e  x^  y  'e  e  z  e  'e 

_  p 

+  cos  0  cos  0  +  R  sin  0  cos  0  ) 

6  6  Z  6  © 


=  E  [R  sin  0  +  cos  0  (R  sin  0  +  cos  0  )  . 

1  x^^  z  *e  e  '  y 


(2.65) 


Similarly : 


cos  0 

cos  0 

E 

e 

e 

*1 

cos  0 

r  p 

sin  0 

'  p 

E~ 

C 

C 

X2 

sin  0 

E~ 

v 

V 

3 

cos  0 

cos  0 

1 

e 

e 

cos  0 

sin  0 

R~ 

e 

e 

y 

sin  0 

e 

R" 

z 

sin  0  cos  0 
e  'e 


sin  0  sin  0 
e  e 


cos  0 


sin  0  cos  0 
e  e 


sin  0  sin  0 
e  e 


cos  0 


(2.66) 


=  E  [-1  (cos  0  sin  0  +  sin  0  sin  0  )  +  R-  (cos2  0  cos  0 

c.  ©  ©  ©  0  y  ,  g  0 


.  +  sin  0  cos  0  )  -  R  (-  sin  0  cos  0  sin  0  cos  0 
©  e  z  ©  6  6  © 


+  sin  0  cos  0  sin  0  cos  0  ) 1 
e  e  e  e 


(2.67) 
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x"  -JT  <see  0  -  din  0 
2  y  9  e 


x"  =  ET 


• 

cos  0  cos  0 
e  e 

-  sin  0 

e 

1 

cos  0  sin  0 
e  e 

cos  0 

e 

R“ 

y 

(2.68) 

sin  0 

e 

0 

IT 

z 

'=  E  [sin  0  (R  cos  0  sin  0  -  R  sin  0  ) 

.  v  '  a  rr  r  ~  _  if  a  ' 


e  z 


e  y 


+  cos  0  (R  cos  0  cos  0  -  sin  0  )] 

e  z  e  e  e 


=  E  [R  cos  0  -  R  sin  S  sin  0  -  sin  0  cos  0  ]  (2.69) 

Xj_  z  ey  6  e  e  e 

x"  =  E  [R  cos  0  -  sin  0  (R  sin.0  +  cos  0  )  (2.70) 

3x1z  e  e  y  e  e 


The  x^  axis  is  now  aligned  with  the  normal,  and  the  vector  isy 
ideally,  zero,  or  very  small  in  an  actual  computation. 

The  polarization  of  this  total  field  is  then: 


(2.71) 


These  four  quantities,  Q^,  Pg,  0g,  and  then  describe  completely  the 
total  incident  field. 

3.  ANTENNA  RESPONSE.. '  In  order  to  compare  the  calculations  with 
measurements,  we  must  be  able  to  express  them  in  a  manner  that  will  yield 
induced'- antenna  voltages.  We  will  now  consider  the  response  of  the  crossed 
(mutually  perpendicular)  vertical  loops  and  vertical  monopole  system  most 
commonly  used  in  whistler  and  vlf  work.  The  coordinate  system  and  field 
directions  were  shown  in  Fig.  2.3*  Figure  2.2  shows  the  arrangement  of 
the  E  and  H  components  of  the  incident  and  reflected  waves. 
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responds  of  n  loop  ontonaa  to  a  sinueoidal  wave  is  proportional 
to  the  time  rate  of  change  in  the  3  field  threading  the  loop.  We  will 
use  the  H  field  to  determine  the  antenna  response,  and  then  express  the 
result  in  terms  of  the  E  field.  The  response  of  a  vertical-loop  antenna, 
making  an  angle  0  with  the  plane  of  incidence  of  the  n-hop  rays  is  given 
by: 


Vi  -  K[(HLi  +  Hir)  cos  0  -  ( Hj, ±  sin  0^  -  H|(r  sin  pjsin  0](2.72) 

The  subscripts  i  and  r  refer  to  the  incident  and  reflected  waves,  res¬ 
pectively.  This  equation  may  then  be  expressed  in  terms  of  the  E  field 
by: 

v/  "  hi  i(l  +  R„)  cos  0  -  EAl(l  -  Ra)  sin  Pn  sin  0]  (2.73) 

where  R|(  is  the  Fresnel  ground  reflection  coefficient  for  parallel  (vertical) 
polarization 

Ra  is  the  Fresnel  ground  reflection  coefficient  for  perpendicular 
(horizontal)  polarization 

h^  is  the  effective  "height"  of  the  loop  antenna  in  meters 
The  effective  height  of  the  loop  contains  the  constant  factor  K 
from  the  previous  expression  Eq.  (2.72)  as  well  as  the  factor  expressing 
the  ratio  of  E  to  H. 

The  preceding  expression  may  then  be  generalized  for  a  pair  of  crossed 
vertical  loops  and  for  a  number  of  waves  arriving  at  various  elevation 
angles  Pn  to  yield: 

vns  =  h£  (x3a  COS  0  "  X2a  8ln  0)  (2'74) 

vgw  =  h^  (x3a  sin  0  +  x2a  cos  0)  (2.75) 

The  vertical  monopole  antenna  responds  only  to  the  parallel  (or 
vertical)  component  of  the  total  electric  field,  which  is  proportional 
to  the  sum  of  the  incident  and  reflected  waves  a$  may  be  seen  from 
Fig.  2.2.  The  induced  voltage  vy  in  a  vertical  antenna  is  thus  given  by: 
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where: 


*2a  -  2  Bin>  t1  -  'In'  Sl”8n 

n 

x3a  ‘  2  Ei|n>  C1  +  Rlln> 

J  n 

x3y  ■  2  El|n)  t1  +  *«„)  008  Pn 
J  n 


(2.76) 

(2.77) 

(2.78) * 

(2.79) 


k.  THE  ARRIVAL  ANGLES  AND  POLARIZATION.  We  have  now  expressed  the 
induced  voltages  in  the  three  antennas  as  a  function  of  the  sum  of  the 
arriving  waves.  In  order  to  simplify  the  interpretation  of  these  results, 

we  now  seek  an  equivalent  plane  wave  which  will  produce  the  same  response 

voltages  as  the  "complex"  wave.  Let  us  rewrite  Eq.  (2.7^),  (2.75)>  (2.76) 

in  terms  of  an  equivalent  plane  wave  at  the  antenna,  and  use  the  subscript 

a  to  indicate  this  wave.  We  shall  first  define  a  polarization  Q  ,  which 

SL 

may  be  expressed  as: 


JY 


(2.80) 


We  must  also  define  an  effective  vertical  angle  of  arrival  B  and  an 
effective  azimuth  angle  of  arrival  0  .  Now,  we  write  Eqs.  (2.7k),  (2.75)> 

€t 

(2.76)  in  terms  of  these  equivalent  antenna  quantities,  and  combine  them, 
which  yields: 


ai  +  Jbi 


a2  *  Jb2 


“R  +  JbF 


(cos  0  -  Q  sin  B  sin  0  ) 

'  ~  a  d  n 


cos  pa 

( sin  0  +  Q  sin  B  cos  0  ) 
' _  a  a _ a _  a 

cos  B 

(cos  0  -  Q  sin  B  sin  0  ) 
_ &  8  ,  &  & _ 

(sin  0&  +  Qa  sin  Bfi  cos  0&) 


(2.81) 


(2.82) 


(2.83) 
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Sotloe  that  the  right  hand  fid®  of  these  equation*  Is  written  entirely 

in  terns  of  the  effective  antenna  quantities,  while  the  left  hand  side 
is  written  entirely  in  terns  of  the  actual  fields  present  at  the  antenna. 
We  will  nov  separate  the  real  and  imaginary  parts,  obtaining 


al  = 


cos  0  -  X  sin  B  sin  0 

a  a _ _a _ _a 

cos  0 
'  « 


(2.84) 


-Y  sin  p  sin  0 

b  =  — - - - - -  =  -  Y  tan  0.  sin  0 

1  „  a  a  a 


cos  P 

a 


(2.85) 


S2  = 


sin  0  +  X  sin  0  cos  0 
a  a  a  a 

cos  3 


(2.86) 


+  Y  sin  P  cos  0 


=  +  Y  tan  P  cos  0 


cos  p 


(2.87) 


From  Eq.  (2.85)  and  Eq.  (2.87)  it  is  easily  seen  that: 

\ 

tan  0  *  -  — 


2  ■ 


Equations  (2.84)  aid  (2.86)  yield: 


(2.88) 


cos  0  -  a,  cos  p  a_  cos  p  -  sin  0 

a  1  a  2  a  a 


sin  0 


cos  0 


a 


(2.89) 


2  2 
cos  0  a,  cos  P  cos  0  =  a_  cos  P  .sin  0  -  sin  0 

■a  1  a  .  a  2  a  a  a 

‘  2  2 

cos  p  (a.  cos  p  +  a  sin  0  )  =  sin  0  +  cos  0 
a  x  a  cl  a  **  a 


cos  P  =  — - 

a  a  sin  0  +  a  cos  0 


(2.90) 


a  1 


a 
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Equation  (2.85)  aad  (2.8?)  yield; 


tan  6  sin  0  tan  0  cos  0 
a  a  a  a 


(2.91) 


and  Eq.  (2.84)  and  (2.86)  yield: 


cos  0  -  a,  cos  0  a  cos  0  -  sin  0 

x  - a — i - a  =  J. - 2 - a  (2.92) 

a 


sin  6  sin  0 
Ka  a 


sin  6  cos  0 
Ka  a 


These  sets  of  equations  therefore  enable  the  four  unknown  quantities 

0  ,'0  ,  X  ,  and  Y  to  be  determined  for  a  plane  wave  when  voltages  in- 
a  a  a  a 

duced  in  the  antennas  are  known.  There  are  conditions  for  which  the 

values  of  the  fields  x^  ,  x„  .  and  x_  are  not  compatible  with  the 

2a'  3a'  3v 

assumption  that  the  field  present  at  the  antennas  is  a  plane  wave.  This 

phenomenon  is  discussed  on  p.  138. 

5-  GONIOMETER  RESPONSE  AND  ANTENNA  VOLTAGE  RATIOS.  The  response"  . 

of  a  goniometer  system  utilizing  crossed  vertical  loops  is  derived  on 

p.101.  0  is  the  goniometer  bearing  angle  for  maximum  signal  measured 

S 

clockwise  from  north.  This  corresponds  to  the  plane  of  the  electrically 

rotated  loop  lying  in  the  plane  of  incidence  and  receiving  maximum 

signal  from  a  vertically  polarized  wave  arriving  from  that  direction. 

The  calculations  are  arranged  to  give  the  antenna  responses  v  and  v 

ns  gw 

in  10°  steps  for  0  from  0  to  l80°.  This  is  equivalent  to  changing  the 
®  o 

direction  of  arrival  in  10  steps  in  the  opposite  rotation.  In  addition, 

the  ratios  v  ,  v  ,  and  v  ,  which  are  defined  by  Eq.  (2.81)  to  (2.83),  are 
n  .L  c. 

calculated  for  the  specified  values  of  0  .  They  correspond  to  the 

g 

quantities  measured  by  Lissajous  patterns  when  comparing  the  outputs  of 
the  ns  and  ew  loop  antennas  and  the  vertical  antenna  on  an  oscillograph, 
v^  corresponds  to  the  common  d-f  trace  often  used  for  sferics  d-f  systems 
at  vlf.  In  all  following  calculations,  it  will  be  assumed  that  the  North 
South  (ns)  loop  lies  in  the  plane  of  incidence,  and  the  East-West  (ew) 
loop  perpendicular  to  it.  These  are  sometimes  referred  to#  as  the  "normal 
and  "abnormal"  loops,  respectively. 
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III.  «*K1SULJ6  OP  TH3  CALCtJLATIOire 


®ie  majority  of  the  curves  shoving  the  results  of  the  calculations 
are  presented  in  this  chapter.  Similar  curves  are  grouped  together  to 
facilitate  comparison.  The  discussion  of  these  results  is  given  in 
Chapter  V. 

A.  THE  CLASSIFICATION  OF  THE  CALCULATIONS 

•There  are  seven  independent  variables  involved  in  these  calculations 
which  may  be  organized  into  three  groups:  (l)  The  normalized  factors 
t,  and  h',  which  describe  the  ionosphere  model.  (2)  The  ground  con¬ 
stants  k  and  a.  (3)  The  path  and  wave  constants  f  and  D.  The  calcul¬ 
ations,  have  been  organized  as  follows:  The  ionosphere  model  has  been 
chosen  from  available  literature  (Wait,  ■  Budden,  etc.)  giving  values 

for  typical  conditions.  The  majority  of  the  calculations  were  made 

5  0 

using  the  "summer  night"  constants.  These  are  =  5  x  10  ,  t  =  60  , 

and  h'  =  90  km.  These  values  correspond  to  effective  values  of  f„,  f 

n  O 

and  v  for  a  homogeneous 
specific1,  one  must  be 
is  fjj,  since,  its  value 
1.3^  Me  is  assumed,  the  values  correspond  to  fQ  =  352; 0  kc,  and 
v  =  4.855  x  10  .  These  are  physically  reasonable  values  for  the  night¬ 
time  layers.  A  "spring"  night  model  might  have  been  a  little  more 
typical  of  the  whistler  characteristics  discussed  in  this  .report  but 
these  established  values  appeared  to  be  the  most  reasonable  starting 
point.  Most  calculations  were  based  on  this  model.. 

In  additions,  some  calculations  were  made  for  the  summer-day  model. 
This  model  was  made  up  from  Wait' s’^'^  value  of  <jQr  =  1.5  x  10^,  t  =  15°, 

and  h'  =  70  km.  The  last  two  values  come  from  the  available  data  on  the 

2k  25 

D-layer  characteristics.  ’  This  corresponds  to  a  value  of  f„  =  1.35 

7  Zk  25 

Me,  f  =  350  kc,  and  v  =  3 .165  x  10  .  The  available  information  * 

on  D-layer  characteristics  contains  a  fairly  wide  range  of  choice  for 
these  factors.  The  values  were  chosen  to  represent  an  equi valent  homo¬ 
geneous  model,  and  probably  most  resemble  the  characteristics  of  the 


model.  Since  only  ratios  of  f„,  f  and  v  are 

rl  O 


assumed  to  find  the  others, 
is  essentially  constant. 


The  best  refi 


If  an  effective  f  of 
n 


Actual  layer  some  distance  above  the  as turned  virtual  height  of  70  km. 

Studies  of  the  results  for  variations  in  ai  f  t,  and  h*  about  the  summer- 

r 

night-model  values  were  also  made.  This  will  be  discussed  further  on 
p.  128. 

The  ground  constants  used  were  picked  for  values  typical  of  three 

general  cases  for  most  of  the  calculations.  These  were  (l)  poor  ground, 

-4  -2 

k  *  3.0,  0  *  1  x  10  mho/m,  (2)  average  ground,  k  =  20,  0  =  3x10  mho/m 

and  (3)  sea  water  k  =  80,  a  =  4.5  mho/m.  These  values  were  chosen  from 

32  26  '46 

the  constants  given  by  Terman,  Wait  and  others.  The  majority  of 
the  calculations  used  the  extreme  values  of  poor  ground  and  sea  water, 
which  seemed  to  be  of  the  most  interest  and  the  most  economical  for 
calculation.  The  results  for  average  ground  would  lie  approximately  mid¬ 
way  between  these  two  extremes.  Some  cases  were  calculated  with  ground 
constants  k  =  0  and  very  large  values  of  a,  which  approaches  perfect 

conductivity.  It  was  soon  determined  that  the  computer  program  used 

14 

would  not  handle  values  of  a  greater  than  1  x  10  . 

The  electric  field  strength  of  the  circularly  polarized  energy  in 
the  ionosphere  Ej^,  was  specified  as  2  +  JO  for  the  majority  of  these 
calculations.  This  is  an  arbitrary  normalizing  factor,  and  the  amplitude 
and  phase  of  the  total  field  and  antenna  response  are  relative  to  this 
value.  The  phase,  of  course,  has  an  ambiguity  of  +  360  p  degrees,  where 
p  is  an  integer. 

The  values  of  f  were  chosen  to  cover  the  major  part  of  the  whistler 

range,  from  500  cps  to  20  kc.  The  values  of  D  were  chosen  to  cover  the 

region  of  interest,  which  is  zero  to  2000  km.  The  latter  value  is  the 

2  3 

maximum  distance  mentioned  by  Storey  *  for  most  whistler  sources.  Most 
of  the  calculations  were  made  as  a  function  of  frequency  for  a  given 
distance.  Values  of  D.  =  0.1,  500,  1000,  1500  and  2000  km  were  originally 
chosen  as  being  typical  in  whistler  work.  The  first  was  essentially 
vertical  incidence.  The  next  two  lie  well  within  the  critical  distance 
for  most  ionosphere  models,  and  the  last  lies  between,  but  well  away 
from,  the  critical  distances  for  the  n  =  0  and  n  =  1  components.  One 
set  of  calculations  for  a  fixed  frequency  was  performed,  covering  a  range 

of  D  from  0.1  to  at  least  12000  km.  These  results  will  be  discussed  on 

page  132. 
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The  final  set  of  ealculation*  vets  performed  for  a  distance  of  14-7.23  lot. 
This  distance  was  chosen  so  that  the  n  *  0  ray  is  parallel  to  the  magnetic 
field  direction  for  the  latitude  of  Stanford  (44°N  geomagnetic),  in  order 
to  determine  the  characteristics  of  a  signal  arriving  from  the  direction 
of  the  magnetic  field.  The  calculated  curves  are  presented  together  to 
facilitate  comparison. 

V 

B.  FIELD  STRENGTH 

i 

The  field  strength  data  are  presented  in  Figs.  3-1  to  3.21  and  in 
3.64.  The  plotted  magnitude  is,  of  course,  relative  to  the  magnitude 
of  the  electric  field  in  the  ionosphere.  The  absolute  field  strength 
cannot  be  determined  without  knowledge  of  the  duct  size.  The  calculations 
and  curves  therefore  are  suitable  only  for  the  study  of  variations  in 
the  field  as  a  function  of  frequency,  distance,  and  ionosphere  and  ground 
parameters.  The  same  amplitude,  frequency,  and  distance  scales  will  be 
used,  when  the  range  of  the  results  permits  this,  to  facilitate  comparison 
of  different  curves 

1.  DISTANCE  VARIATIONS.  Figures  3.1  through  3-7  show  the  magntiude 

of  the  field  seen  by  the  normal  loop  ( |vns|  for  0 ^  =  0)  and  the  vertical 

monopole  ( | v^ | )  as  a  function  of  frequency.  The  curves  are  calculated 

for  poor  ground  and  sea  water,  using  the  summer-night  model  ionosphere, 

for  distances  of. 0.1,  500,  1000,  and  2000  km.  Figures  3.8  through  3*12 

show  the  magnitudes  of  the  fields  seen  by  the  normal  and  abnormal  loops 

(  |v  I  and  |v  I  for  0  =  0)  and  by  the  vertical  monopole  as  a  function 

\  1  ns 1  1  ew 1  /g 

of  distance  for  5  4c  and  15.5  4c.  Figure  3*8  was  plotted  from  results 
for  poor-ground  conditions,  and  the  rest  for  sea  water.  Figures  3*8 
through  3*10  were  plotted  from  results  calculated  for  the .summer -night 
model  ionosphere,  and  3«H  and  3*12  for  the  summer-day  model.' 

2.  IONOSPHERE  HEIGHT  VARIATION.  Figures  3-13  and  3.14  show  the 
normal-loop  and  vertical -monopole  fields  as  a  function  of  frequency  for 
poor  ground  at  a  distance  of  500  km.,  The  summer-night  model  ionosphere 
was  used,  except  that  h1  was  changed  to  89  and  91  4m,  respectively,  in 
Figs.  3.13  and  3.14. 
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FIG.  3.2.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP 
AND  MONOPOLE  VS  FREQUENCY,  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE, 

SEA  WATER,  D  *  0.1  km,  AND  4g  *  0. 
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FIG.  3.3-.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  AND  MONOPOLE  VS  FREQUENCY 
FOR  SUMMER-NIGHT  MODEL  IONOSPHERE,  POOR  GROUND,  D  =  500  km,  AND  <£  =  0. 


FIG.  3.4.  MAGNITUDE  OF  E- 
SUMME  R-NIGHT  MODEL  IONOSPH 


f  (kc) 

FIG.  3.5.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  AND  MONOPOLE  VS  FREQUENCY 
SUMMER-NIGHT  MODEL  IONOSPHERE,  SEA  WATER,  D  -  1000  ka,  AND  <£  ■  0. 
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FIG.  3.6.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  AND  MONOPOLE  VS  FREQUENCY 
SUMMER-NIGHT  MODEL  IONOSPHERE,  POOR  GROUND,  D  -  2000  ka,  AND  <£  -  0. 


FIG.  3.7.  MAGNITUDE  OF  E- FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  AND  MONOPOLE  VS  FREQUENCY,  FOR 
SUMMER-NIGHT  MODEL  IONOSPHERE,  SEA  WATER,  D  -  2000  k»,  AND  <£  -  0. 


FIG.  3.12  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  AND  ABNORMAL 
VERTICAL  LOOPS  AND  MONOPOLE  VS  DISTANCE,  FOR  SUMMER-DAY  MODEL 
IONOSPHERE.  SEA  WATER,  f  ■  15.5  kc,  AND  <£  ■  0. 
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FIG.  3.13.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  AND  MONOPOLE  VS  FREQU8NCY  FOR 
SUMMER-NIGHT  MODEL  IONOSPHERE,  EXCEPT  h'  *  89  k»,  POOR  GROUND,  D  -  500  k«,  AND  <£  =  0. 


SUMMER-NIGHT  MODEL  IONOSPHERE  EXCEPT  h  =  91  km,  POOR  GROUND, 
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FIG.  3.16.  MAGNITUDE  OF  E- FI  ELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  AND  MONOPOLE  VS  FREQUENCY 
SUMMER-NIGHT  MODEL  IONOSPHERE,  EXCEPT  u  =  5 x  lO6,  POOR  GROUND,  D  =  500  km,  AND  =  0. 
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FIG.  '3.1'.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  VS  FREQUENCY  FOR  SUMMER* NIGHT 
MODEL  IONOSPHERE,  SEA  WATER,  D  -  S00  Vm,  AND  (£  -  30. 


FIG.  3.18.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  VS  FREi 
MODEL  IONOSPHERE.  SEA  WATER,  D  *  500  ki,  AND<£  *  60. 


MODEL  IONOSPHERE,  SEA  WATER,  D  *  500  lea,  AND  <£a  *  90. 


f  (kc) 

FIG.  3.20.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  VS  FREQUENCY  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE,  SEA  HATER,  D  *  500  km,  AND  d  *  120. 
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FIG.  3.21.  MAGNITUDE  OF  E-FIELD  RECEIVED  BY  NORMAL  VERTICAL  LOOP  VS  FREQUENCY  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE,  SEA  WATER,  D  -  500  lea,  AND  <f>  -  150. 


3.  IONIZATION  DENS  IT*  VARIATION.  Figures  3. 15  aad  3.16  are  similar 

to  Fig.  3*13  and  Fig,  3.14  except  that  the  values  of  u>  have  been  changed 

4  r  & 

instead  of  those  of  h'.  The  values  used  are  5  x  10  and  5  x  10  for 
Figs.  3.13  and  3*14,  respectively,  a  factor  of  10  lover  and  higher  than 
for  the  summer -night  model. 


4.  RELATIVE  BEARING  VARIATION.  Figures  3.17  through  3.21  show 

the  field  strength,  as  seen  by  the  normal  (ns)  loop  as  a  function  of 

frequency  for  the  summer-night  model  ionosphere  and  sea  water  at  500  km. 

The  value  of  0  was  increased,  in  30-degree  steps  in  each  case,  from  30 
8 

through  150  degrees.  These  curves  show  the  effect  of  different  bearings 
relative  to  a  fixed  loop. 


C.  CROSSED -LOOP  VOLTAGE  RATIO  VARIATION  WITH  FREQUENCY 


1.  DISTANCE  VARIATION.  Figures  3-22  through  3.3O  show  the  values 
of  the  real  and  imaginary  parts  of  v  (v  =  v  /v  )  asi  a  function. of 

n  K  ns  GW 

frequency  for  various  conditions.  Figures  3-22  through  3*26  show  the 
values  calculated  for  the  summer -night  model  ionosphere  with  sea  water 
and  poor  ground  conditions  for  500,  1000,  and  2000  km. 

2.  IONOSPHERE  HEIGHT  VARIATION.  Figures  3.27  and  3.28  show  the 
values  of  the  real  and  imaginary  parts,  of  v  for  poor  ground  at  500  km, 
using  the  summer -night  model  ionosphere  except  that  h1  was  changed  to 
89  and  91  km,  respectively. 

3.  IONIZATION  DENSITY  VARIATION.  Figures  3.29  and  3.30  show  the 

values  of  the  real  and  imaginary  parts  of  v_  for  the  same  conditions  as 

R  4 

2,  above,  except  that  the  values  of  u>  were  changed  to  5  x  10  and 

6  r 

5  x  10 ,  respectively,  while  h'  was  90  km. 


D.  THE  WAVE  POLARIZATION  Q 

fit 


The  values  of  the  real  and  imaginary  parts  of  Q  ,  the  polarization 

Q 

of  the  incoming  wave  as  seen  by  the  crossed  loop  and  vertical  system  are 
presented  as  a  function  of  frequency  in  Figs.  3*31  to  3*34.  It  will  be 
recalled  that  Q  is  not  defined  for  conditions  for  which  there  is  no 
equivalent  plane  wave.  No  values  are  defined  for  these  frequencies. 
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FIG.  3.22.  REAL  AND  IMAGINARY  PARTS  OF  CROSSED- LOOP  VOLTAGE  RATIO 
VS  FREQUENCY  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE,  POOR  GROUND, 

D  ■  500  km,  AND  <£g  ■  0. 
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FIG.  3.24.  REAL  AND  IMAGINARY  PARTS  OF  CROSSED-LOOP  VOLTAGE  RATIO 
VS  FREQUENCY  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE,  SEA  WATER,  D  ■  1000  ka, 
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FIG.  3.25  •  REAL  AND  IMAGINARY  PARTS  OF  CROSSED- LOOP  VOLTAGE  RATIO  VS 
FREQUENCY  FOR  SUMMER  NIGHT  MODEL  IONOSPHERE,  POOR  GROUND,  D  *  '2000  km, 
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FIG.  3.26.  REAL  AND  IMAGINARY  PARTS  OF  CRQSSED-LOOP  VOLTAGE  RATIO 
VS  FREQUENCY,  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE,  SEA  WATER,  D  ■  2000  km, 
AND  <p  •  0.  l- 
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FIG.  3.27  REAL  AND  IMAGINARY  PARTS  OF  CROSSED-LOOP  VOLTAGE  RATIO  VS 
FREQUENCY,  FOR  SUMMER  NIGHT  MODEL  IONOSPHERE  EXCEPT  h'  «  89  kn,  POOR 
GROUND,  D  -  500  km,  AND  4  -  0. 
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FIG.  3.28.  REAL  AND  IMAGINARY  PARTS  OF  CROSSED-LOOP  VOLTAGE  RATIO 
VS  FREQUENCY,  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE  EXCEPT  h'  ■  91  k», 
POOR  GROUND,  D  ■  S00  kn,  AND  <£  ■(). 
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FIG.  3.29.  REAL  AND  IMAGINARY  PARTS  OF  CROSSED-LOQP  VOLTAGE  RATIO  VS 
FREQUENCY,  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE  EXCEPT  U(  -  5  x  104,  POOR 
GROUND,  D  ■  500  km,  AND  <t.  •  0. 
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FIG.  3.30.  REAL  AND  IMAGINARY  PARTS  OF  CROSSED-LOOP  VOLTAGE  RATIO  VS 
FREQUENCY,  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE  EXCEPT  U>r  -  5  x  1 06 ,  POOR 
GROUND,  D  »  500  km,  AND  $  ■  0. 
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FIG  3.31.  REAL  AND  IMAGINARY  PARTS  OF  POLARIZATION  Q,  AT  ANTENNA 
VS  FREQUENCY,  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE,  SEA  WATER, 

D  *  500  km . 
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FIG.  3.33.  HEAL  AND  IMAGINARY  PARTS  OF  POLARIZATION  Q,  AT  ANTENNA 
VS  FREQUENCY.  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE  EXCEPT  *  5  x  104, 
POOR  GROUND,  D  -  530  In. 
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FIG.  3.34.  REAL  AND  IMAGINARY  PARTS  OF  POLARIZATION  Q,  AT  ANTENNA 
VS  FREQUENCY,  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE  EXCEPT  uf  -  5  x  106, 
POOR  GROUND,  D  *=  500  km. 
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Figures  3«31  and  3.32  show  the  real  and  imaginary  parts  of  Q  for 

fit 

the  summer-night  model  and  sea  water  conditions,  at  500  and  2000  km 
respectively.  Figures  3-33  and  3*34  show  these  values  for  poor  ground 
at  500  km.  The  summer-night  model  ionosphere  is  used,  except  that 
os  =  5  x  10  and  5  x  10  ,  respectively. 

E.  THE  APPARENT  BEARING  ANGLE  0 

a . 

Figures  3-35  through-  3-45  show -the  values  of  0  vs  frequency  for 

& 

poor  ground  and  sea  water  conditions. 

1.  DISTANCE  VARIATION.  Figures  3*35  through  3-41  show  the  values 

of  0  vs  frequency  for  the  summer-night  model  ionosphere  at  distances  of 
& 

0.1,  500,  1000,  and  2000  km. 

2.  IONOSPHERE  HEIGHT  VARIATIONS.  Figures  3.42  and  3.43  show  the 

calculated  values  of  0  vs  frequency  for  poor  ground  at  500  km.  The 

& 

summer-night  model  ionosphere  was  used,  except  that  h'  was  changed  to 
89  and  91  km,  respectively. 

3.  IONIZATION  DENSITY  VARIATIONS.  Figures  3.44  and  3.45  show  the 
values  of  0  vs  frequency  for  poor  ground  conditions  at  500  km.  Summer- 

£L 

night  model  ionosphere  conditions  are  used  except  that  cjd  is  changed  to 

4  6  r 

5  x  10  and  5  x  10  ,  respectively. 

F.  TRANSMISSION.  COEFFICIENTS 

The  calculated  magnitude  of  the  two  transmission  coefficients 
(  |T|j  |  and  |TA  | )  is  shown  for  frequencies  of  0.5,1,  5>  10,  15,  and  20  kc. 
The  curves  are  plotted  as  a  function  of  6  ,  which  is  equivalent  to  6 , 
the  angle  of  transmission  into  the  waveguide.  Figures  3*46  through 
3.48  show  |t(,  I  for  three  different  ionization  densities,  and  Figs.  3.49 
through  3.51  show  I T . I  for  these  same  conditions.  The  conditions  are 
05^  =  5  x  10  ,  5  x  10^  ( summer -night  model)  and  5  x  10  ,  with  t  *  60 
degrees. 
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FIG.  3.36.  LATERAL  DEVI  ATION  ANGLE  <£,  VS  FREQUENCY,  FOR  SUMHC  R- NIGHT 
MODEL  IONOSPHERE,  SEA  WATER,  D  ■  0.  1  k». 


FIG.  3.37.  LATEBAL  DEVIATION  ANGLE  <£a  VS  FREQUENCY,  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE,  POOR  GROUND,  D  «  500  km. 
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FIG.  3.38.  LATERAL  DEVIATION  ANGLE  VS  FREQUENCY,  FOR  SUMMER-NIGHT 

MODEL  IONOSPHERE,  SEA  WATER,  D  ■  500  km. 
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FIG.  3.39.  LATERAL  DEVIATION  ANGLES,  VS  FREQUENCY,  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE,  SEA  WATER,  D  -  1000  km. 
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FIG.  3.40.  LATERAL  DEVIATION  ANGLE  <£,  VS  FREQUENCY,  FOR  SUMNER-NIGHT 
MODEL  IONOSPHERE,  POOR  GROUND,  D  ■  2000  ka. 
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FIG.  3.41.  LATERAL  DEVIATION  ANGLE  VS  FREQUENCY  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE,  SEA  WATER,  D  -  2000  kn. 
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FIG.  3.42.  LATERAL  DEVIATION  ANGLE  VS  FREQUENCY  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE  EXCEPT  h'  »  89  kn,  POOR  GROUND,  D  ■  500  kn. 
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FIG.  3.43.  LATERAL  DEVIATION  ANGLE  <£,  VS  FREQUENCY  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE  EXCEPT  h'  -  91  ka,  POOR  GROUND,  D  -  500  ka. 


(deg) 


FIG.  3.44.  .  LATERAL.  DEVIATION  ANGLES,  VS  FREQUENCY,  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE  EXCEPT  U)r  -  5  x  104,  POOR  GROUND,  D  ■  500  Its. 


FIG.  3.45.  LATERAL  DEVIATION  ANGLES,  VS  FREQUENCY,  FOR  SUMMER-NIGHT 
MODEL  IONOSPHERE  EXCEPT  «r  •  !x  106,  POOR  GROUND,  D  -  500  k». 
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FIG.  3.47.  MAGNITUDE  OF  TRANSMISSION  COEFFICIENT  T(,  VS  0n  FOR 
r  ■  60°.  tdr  ■  5  x  105  (SUMMER-NIGHT  MODEL  IONOSPHERE). 


-  86  - 


0  10  20  30  40  50  60  70  80  90 

6n  ( deg) 

FIG.  3.48.  MAGNITUDE  OF  TRANSMISSION  COEFFICIENT 
T ||  VS  £?n  FOR  r  -  60°,  wr  ■  5  x  106. 
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FIG.  3.49.  MAGNITUDE  OF  TRANSMISSION  COEFFICIENT  ^  VS 
8n  FOR  r  «  60°,  Wr  ■  5  x  104. 

,  88  - 


<9„  (deg) 


FIG.  3.50.  MAGNITUDE  OF  TRANSMISSION  COEFFICIENT  T  VS  &n  F.OR 
r  •  60°,  ur  -  5  x  10®  (SUMMER- NIGHT' MODEL  IONOSPHERE). 
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FIG.  3.51.  MAGNITUDE  OF  TRANSMISSION  COEFFICIENT 
Tx  VS  8n  FOR  T  «  60°,  Wr  -  5  x  106. 
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G.  IONOSPHERE  REFLECTION  COEFFICIENTS 


The  calculated  magnitudes  of  the  ionosphere  reflection  coefficients 
are  shown  in  Figs.  3*52  through  3*59  fis  a  function  of  the  ionosphere 
angle  of  incidence  9 for  frequencies  of  0.5*  1*  5*  10*  15*  20  kc. 

They  are  grouped  by  ionosphere  conditions.  Figures  3*52  through  3-55 
show  the  values  for  to  =  5  x  10  ,  and  Figs.  3*56  through  3*59  show  those 

c  ** 

for  o>r  =  5  x  10^  (summer-night  model),  with  t  ■  60  degrees. 

H.  GROUND  REFLECTION  COEFFICIENTS 

The  calculated  magnitude  of  the  two  ground  reflection  coefficients 
is  shown  in  Figs.  3*60  through  3*63  as  a  function  of  p  ,  the  takeoff- 
arrival  angle  at  the  earth's  surface.  Frequencies  of  0.5*  1*  5*  10*  15* 
and  20  kc  are  generally  shown.  Figures  3-60  and  3-6l  show  the  values 
for  poor  ground,  and  Figs.  3*62  and  3-63  show  the  values  for  sea  water. 

I.  THE  CHARACTERISTICS  OF  A  WAVE  ARRIVING  FROM  THE  MAGNETIC  FIELD 
DIRECTION 

The  curves  in  Figs.  3.64  and  3.65  were  calculated  for  a  distance 

of  47.23  km.  At  this  distance  the  n  =  0  ray  is  parallel  to  the  earth's 

magnetic  field  at  the  geomagnetic  latitude  of  Stanford  (44°N).  Summer- 

night  model  ionosphere  and  sea  water  conditions  were  vised.  The  calculated 

magnitude  of  the  field  seen  by  the  antennas  is  shown  as  a  function  of 

frequency  in  Fig.  3.64.  Figure  3.65  shows  the  real  and  imaginary  parts 

of  v  for  the  same  conditions. 
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FIG.  3.52.  MAGNITUDE  OF  IONOSPHERE  REFLECTION  COEFFICIENT 
II  R||  VS  &n ,  FOR  T  *  60° ,  AND  Ur  *  5  x  104. 


FIG.  3.53.  MAGNITUDE  OF  IONOSPHEHE  REFLECTION  COEFFICIENT 
||  R  JL  F0R  T  "  60°,  AND  (i)r  ■  5x  104 
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FIG.  3.54.  MAGNITUDE  OF  IONOSPHERE  REFLECTION  COEFFICIENT 
iR)|  VS  0n,  -FOR  t  ■  60°.  AND  ilr  .  5,  104. 


FIG.  3.55.  MAGNITUDE  OF  IONOSPHERE  REFLECTION  COEFFICIENT 
iRl  VS  9n‘  FOR  T  ■  60°,  AND  U)r  ■  5  x  104. 
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FIG.  3.56.  MAGNITUDE  OF  IONOSPHERE  REFLECTION  COEFFICIENT 


||  R ||  VS  &n,  FOR  T 
MODEL  IONOSPHERE). 
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FIG.  3.57  MAGNITUDE  OF  IONOSPHERE  REFLECTION  COEFFICIENT 
,, R  i  VS  6n,  FOR  r  -  60°,  wt  -  5  *  105  (  SUMMER-NIGHT 
MODEL  IONOSPHERE). 


FIG.  3.59.  MAGNITUDE  OF  IONOSPHERE  REFLECTION  COEFFICIENT 
R.  VS  6n,  FOR  r  -  60°,  U)r  »  5  x  105  (SUMMER-NIGHT  MODEL 
ONOSPHERE). 
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FIG.  3.61.  MAGNITUDE  OF.  GROUND  REFLECTION  COEFFICIENT 
Rj.  vs  fin  ’  F0B  P00B  GROUND,  K  ■  3.0,  cr  ■  1  x  104  bIio/b. 
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FIG.  6.62.  MAGNITUDE  OF  GROUND  REFLECTION  COEFFICIENT  R  (|  VS  /?„  FOR 
SEA  WATER,  K  «  80.0,  a  ■  4.5  mho/m. 
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FIG.  3.  63.  MAGNITUDE  OF  GROUND  REFLECTION  COEFFICIENT  Rx  VS  /Sn  FOR 
SEA  WATER,  k  ■  80.0,  cr  ■  4.5  aiho/ra. 
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FIG.  3  64 
FREQUENCY 


FIG.  3.65.  BEAL 
SEA  WATER,  AND  S 


IV.  THE  d-f  EXPERIMENT 


A.  SOURCE  LOCATION  BY  d-f 

As  we  have  seen  in  Chapter  I,  it  should  be  possible  to  determine 
the  exit  point  at  which  whistler-mode  energy  is  injected  into  the  earth- 
ionosphere  waveguide.  There  are  several  common  methods  for  doing  this. 

The  two  most  common  methods  are  called  "inverse -Loran"  or  time  delay 
location  and  "d-f"  or  direction  finding. 

The  inverse  Loran  method  requires  the  measurement  of  differences  in 
phase  or  time  of  arrival  of  a  signal  at  spaced  stations.  Each  difference 
in  arrival  time  at  two  stations  will  place  the  source  on  a  hyperbola 
drawn  with  the  two  stations  as  foci.  If  more  sets  of  differences  are 
available,  the  location  can  be  determined  at  the  intersection  of  the 
proper  hyperbolas.  This  method  has  the  advantage  of  not  requiring  special 
goniometer  equipment  and  antennas. 

Large  spacings,  of  the  order  of  150  km  or  more,  are  necessary  for 

reasonable  phase  differences  at  vlf.  The  required  precision  is  such 

that  ordinary  tape  recording  of  signals  produces  too  much  error.  The 

only  practical  way  appears  to  be  to  relay  signals  to  a  central  point 

for  direct- recording.  The  problem  of  separation  and  identification  of 

components  in  a  whistler  further  complicates  the  measurement  of  relative 

phase.  In  addition,  the  phase  velocity  of  earth -ionosphere  waveguide 

34  - 

signals  is  not  precisely  known,  which  contributes  more  error.  These 
complications  led  to  the  conclusion  that  this  method  should  be  used  only 
if  simpler  methods  were  not  successful. 

The  d-f  method  relies  on  direction -of -arrival  measurements  at 
different  points  to  locate  the  exit  point  by  triangulation,  d-f  systems 
utilizing  loops  or  Adcock  antennas  are  common  at  the  higher  frequencies. 
These  are  generally  narrow-band  systems  and  use  some  type  of  goniometer. 
The  use  of  wideband  d-f  systems  at  vlf  for  the  location  of  atmospherics 
has  become  widespread.  Wideband  systems  have  the  advantage  of  providing 
better  separation,  in  time,  of  the  various  sky  wave  components  of  im¬ 
pulsive  signals,  except  for  sources  very  close  to  the  receiver.  A 
narrow-band  system  often  yields  confused  results  because  of  its  inherent 
ringing  and  slow  rise  time. 
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A  wideband  system  would  be  advantageous,  since  a  whistler  is  a 

wideband  phenomenon.  The  problem  is  complicated  by  the  existence  of 

different  components  in  the  whistler  which  have  traveled  different  paths 

and  thus  exhibit  different  arrival  angles. 

The  use  of  a  goniometer  system  for  rapidly  rotating  the  antenna 

.  * 

pattern  was  suggested  by  the  author  in  1956,  but  no  practical  method 
of  recording  and  analyzing  the  output  was  available.  It  remained  for 
J.  M.  Watts  to  suggest  that  the  wideband  output  of  the  goniometer  be 
recorded  on  magnetic  tape  in  the  same  manner  as  routine  whistler  record¬ 
ings,  and  that  a  spectrum  analyzer  be  used  to  separate  the  multiple  com¬ 
ponents  and  to  determine  the  location  of. the  nulls.  Watts  was  the  first 

4l 

to  perform  this  experiment. 

This  system  has  the  advantage  of  being  simple  enough  to  permit 
routine  recordings  of  whistler  phenomena.  The  goniometer  system  was 
thus  selected  for  an  experimental  test  program.  It  has  the  advantage 
that  it  can  often  utilize  an  antenna  system  which  is  normally  used  for 
routine  vlf  recordings.  The  present  system  utilizes  the  crossed  vertical 
loops  and  vertical  monopole  antenna  which  were  designed  for  routine 
recording  on  another  experiment. 

B.  THE  GONIOMETER  SYSTEM  RESPONSE 

In  order  to  discuss  the  goniometer  system  itself,  it  is  necessary 
to  examine  the  response  to  vlf  signals  of  the  two  types  of  antennas. 

The  loop  antenna  responds  to  the  magnetic  field  threading  it,  and  so  the 
induced  voltage  may  be  expressed  by: 


nd0 

dt 


-  nA 


dB 

dt 


Now  for  sinusoidal  variation,  —  =  -  J01®  and  B  =  p'H 

dt 


(4.1) 


The  goniometer  equations  stated  in  this  report  are  not  correct, 
•  however. 

— ** 
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Therefore : 


v  =  -  jonAp'H  (4.2) 

where : 

n  =  the  number  of  turns  of  the  loop 

2 

A  =  the  area  of  the  loop  in  m 
H'  =  the  permeability  of  free  space 
03  =  2«f 

f  --  the  wave  frequency 

H  =  the  magnetic  field  strength  of  the  wave 
T)  =  the  intrinsic  impedance  of  free  space 
Now,  substituting  r|  =  E/H  yields : 

vt  = .  (U.3) 

or 

vJ  =  h^  E  (4.4). 

.  .uxiAll  ' 

Where  E  is  the  electric  field  of  the  wave  in  v/m  and  h'  =  -  J  -—'■■■■  is  the 

*  n 

effective  height  of  the  loop  in  meters.  Now  the  induced  voltages  may  be 
expressed  in  terms  of  the  parallel  and  perpendicular  polarized  components 
of  a  wave  incident  at  azimuth  0  and  vertical  angle  e(d  =  90  -  3).  (See 
Fig.  2.2.)  For  the  ns  and  ew  loops,  oriented  at  an  angle  0  to  the  plane 
of  incidence  and  perpendicular  to  this  plane,  respectively,  the  induced 
voltages  are: 

vns  =  (E|jcos0  -  E^  cos©  sin0)  (4.5) 

vew  =  (E||  sin0  +  E^  cose  cos0)  ( 4.6) 


The  goniometer  system  utilizes  a  sine -cosine  capacitor  whose  output 
voltage  is  given  by: 


v  =  v  cos0  •+  v  sin0 
g  ns  g  ew  g 


(*.7) 


where  0  is  the  goniometer  coupling  single,  measured  clockwise  from  zero 
S 

(or  north).  Without  loss  of  generality  we  may  now  assume  that  h^  is 
equal  to  unity.  When  Eqs.  (4.5)  euid  (4.6)  are  substituted  in  Eq.  (4.7) 
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and  the  equation  simplified,  the  result  is: 

Vg  =  E||  cos0  cos0^  -  E^  cos0sin0^  +  E(|  sin0sin0g  +  cos8cos0sin0g 

=  E||  ( COS0CO30  +  sln08ln0  )-  E^  cos@(sin0cos0  -  cos0sin0  ) 

S  8  8  8 

=  E||-cos(0  -  0  )  -  E,  cos@sin(0  -  0  )  (4.8) 

8  6 


Therefore  we  see  that  the  goniometer  output  is  the  same  as  that  for  vng 

with  the  angle  0  in  the  original  expression  replaced  by  0  -  0  ,  and  that 

6 

rotating  the  goniometer  has  the  same  effect  as  rotating  the  pattern  of 
the  antenna  system  or  changing  the  direction  of  arrival  with  respect  to 
the  antenna. 

Now  let  us  express  the  output  voltage  in  terms  of  the  polarization 
Q  of  the  wave,  with  the  following  results: 


where 


v  =  E||  [cos(0-0  )-Qcos0sin(0-0  )  ]  (4.9) 

O  6  ’  D 

Q  =  X  +  JY 


From  the  form  of  this  equation  it  is  evident  that  the  value  of  the  com¬ 
plex  polarization  Q  will  have  a  great  effect  on  the  voltage  of  the 


goniometer  system.  Let  us  find  the  maximum  and  minimum  points  of  this 

voltage  with  respect  to  (0  -  0  ).  The  value  of  (0  -  0  )  corresponding 

8  8 

to  the  minimum  is  the  hearing  error  angle.  The  maximum  and  minimum 


points  are  most  easily  found  by  taking  the  square  of  the  magnitude  of 

Eq.  (4.9),  dividing  through  by  the  square  of  the  magnitude  of  E(|,  and 

then  differentiating  with  respect  to  (0  -  0  ).  This  process  yields  the 

6 

following  equations: 


|v|2=  -S  a [cos (0-0  ) -Xcos6sin(0-0  ) ]+Y^cos20sin2(0-0  ) 

tji  6  6  8 


(*.io) 


=  0  =  [cos(0-0  )-Xcos6sin(0-0  ) ] [-sin(0-0  ) -XcosS cos(0-0  ) ] 
8  8  8  6 


+  Y2cos20sin(0-0  )cos(0-0  ) 

6  8 
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Then 


2 

g-^-y  =  -  Xcose[cos2(0-0g)-8in2(j2(-0g]+8in(0.0g)co8(0.0g)[.l+(X2+Y2)eo820] 

s 

Now  the  substitution  of  the  trigonometric  identities  sinX  cosX  «  sin2X 

2  2  2 
and  cos  X  -  sin  X  =  cos2X  yields : 

0  =  -Xcos6cos2(0-0  )+  i[(X2+Y2)cos20-l]sin2(0-0  ) 

s  _  s 

This  equation  may  be  rearranged  to  yield: 

tan2(0-0  )=+2Xcose/(X2+Y2)cos2e  -1  (4.11) 

g 

This  shows  that  the  polarization  affects  the  bearing  error  angle. 

Note  especially  the  following  property  which  is  present  in  general  and 
is  important,  as  will  later  be  shown.  The  sign  of  Y,  the  imaginary  part 

of  the  polarization,  determines  the  direction  of  the  rotation  of  the 

/ 

phase  shift  Lissajous  pattern  when  the  two  antenna  voltages  are  compared. 
Note  that  Y  appears  only  to  an  even  power,  which  means  that  the  direction 
of  rotation  does  not  affect  the  sign  of  the  bearing  error.  The  sign  of 
X,  the  real  part  of  the  polarization,  does  affect  the  sign  of  the  .error 
since  it  appears  to  an  odd  power. 

Two  special  cases  are  of  interest.  First,  when  X  =  0(Q  =  +  j Y ) , 
and  6  does  not  equal  zero,  then  the  following  results  are  obtained: 


tan2  (0-0  )  =  0 

6 

(4.11a) 

(0“0g)l  =  n* 

(4.11b) 

(0-0g)2  =  (2n  -  1)  | 

(4.11c) 

where  n  =  0,  1,  2,  3> • • • 

Taking  the  first  values  of  each  of  the  above  solutions,  the 
voltages  obtained  from  Eq.  (4.9)  for  these  two  values  is: 

ratio  of 

=  Y  cosQ 

Vgl 

-  10  4  - 


Thus,  the  null  for  Y  cosS  near  unity  is  very  shallow.  The  bearing  error 
is  either  0°  or  90°  when  the  magnitude  of  Y  cos 6  is  greater  than,  or 
less  than,  unity,  respectively.  For  circular  polarization  (Y  ■  +  1), 
the  error  is  always  90°. 

When  the  polarization  is  linear,  Y  *  0.  For  vertical  polarization 
X  =  0  also,  and  Eqs.  (4. 11a)  to  (4.11c)  again  apply.  However,  in  this 
case,  if  the  first  value  of  each  solution  is  taken,  the  result  is: 

Isi 

Vg2 

Thus,  the  expected  results,  that  the  bearing  error  is  zero  and  the 
minimum  is  a  true  null,  are  shown  to  be  correct  for  the  case  of  parallel 
polarization.  Calculations  of  bearing  error  and  null  ratio  have  been 
made  for  a  wide  range  of  polarizations.  The  results  are  discussed  on 
page  110. 

C.  THE  GONIOMETER -VERTICAL  SYSTEM 

One  logical  improvement  to  the  goniometer  system  would  be  to  re¬ 
move  the  180°  ambiguity  by  some  method  of  "sense".  The  polarity  of  the 
output  voltage,  as  the  goniometer  rotates  through  half  a  revolution,  will 
be  opposite  to  that  in  the  other  half.  Therefore,  if  the  output  of  a 
vertical  antenna,  which  does  not  change  with  the  direction  of  arrival, 
were  added  to  the  goniometer  voltage,  one  of  the  two  minima  would  tend 
to  be  canceled  and  the  other  doubled  in  amplitude.  The  exact  amount  of 
cancellation  would  depend  upon  the  vertical  angle  of  arrival  and  the 
polarization,  since  the  induced  voltage  in  a  vertical  antenna  system  is 
given  by: 

v^  =  h^  E||  sin@  (4.12) 

We  can  assume  here  without  loss  of  generality  that  h^.  is  equal  to  unity. 
In  this  case  if  the  angle  9  is  9'"’  and  the  polarization  is  parallel, 
the  combined  output  voltage  will  have  one  null  which  is  wider  than  the 
original  goniometer  null  and  one  maximum  which  is  twice  as  large  as  the 


cos  0 

- >  00 

cos  n/2 
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original  goniometer  maadmitn.  The  null  ratio  will  therefore  be  twice 

that  with  the  goniometer  alone.  If  ©  is  lees  than  90° >  two  equally 

spaced  nulls  will  appear  on  either  side  of  the  true  direction  of  arrival. 

On  further  consideration  of  the  details  of  the  response,  another 

theoretical  advantage  appears.  The  best  nulls  would  tend  to  result  where 

the  goniometer  and  vertical  voltages  are  in  exactly  opposite  phase. 

Consideration  of  Eq.  (4.8)  shows  that  this  must  occur  for  (0  -  0  )«  n,3rt, 

8 

...etc.,  or  when  the  plane  of  the  loop  is  in  the  plane  of  incidence  of 
the  signal,  so  that  the  loop  has  no  response  to  the  perpendicular  com¬ 
ponent  of  electric  field  strength.  It  would  therefore  be  expected  that 
this  system,  under  the  condition  of  a  single  ray,  would  exhibit  less 
polarization  error  and  better  nulls  in  general  than  the  goniometer 
system  alone,  since  there  should  always  be  an  angle  where  the  voltages 
are  exactly  out  of  phase,  and  this  angle  is  the  true  null  bearing.  It 
is  likely  that  subsidiary  nulls  would  appear,  due  to  minima  caused  by 
amplitude  variations  of  the  out-of -phase  signals  at  other  angles  also. 
These  would  tend  to  be  of  lesser  amplitude  than  the  main  null  and  should 
be  relatively  unimportant. 

Let  us  determine  the  response  of  the  system  just  described  by  the 
method  used  to  determine  the  goniometer  response.  After  adding  Eqs. 

(4.9)  and  Eq.  (4.12)  the  voltage  of  this  system  becomes: 

v  =  E||  [sine  +  cos(0-0  )  -  Q  cos©sin(0-0  )]  (4.13) 

gv  g  g 

The  simplest  way  of  determining  the  nulls  is  again  to  take  the  square 
of  the  amplitude  divided  by  the  square  of  the  amplitude  of  E|(,  which 
yields: 

v  2 

-SL  =  [cos(0-0  )+sin©-Xcos©sin(0-0  )  ]^+Y^cos2©sin^(0-0  )  (4.14) 

r>  8  8  6 

k  „ 

Now  let  us  make  the  substitution  a  =  (0-0  ).  Then,  taking  the  derivative 

D 

with  respect  to  a  and  setting  the  result  equal  to  zero,  yields: 

( [ { X2+Y2 ) cos2© -1 ]sina-Xsin0  cos©  }  cosa=+Xcos© ( l-2sin2a) +sin©sinO! 


The  next  step  is  to  substitute  til  -  sin2a  »  cosa  and  square  toe  re stilt, 
which  yields: 

{  [ ( X2+Y2 ) cos20 -1 ]  2sin2cr-2Xsin0 cos0 [ (X2+Y2 ) cos20 -1 ]sina+X2sin29 cos 2 9  } 

.  (l-sin2a) 


.2  2  2  2  2  2 
=a  cos  0(l-2sin  a)  +2Xsin0cos0(l-2sini  a)sina+sin  Ssina 


When  the  multiplications  are  performed  and  the  terms  are  collected  the 

4  o  2 

A  sin  a+BsinJa+Csin  a+Dsina+E  =  0  (4.15). 

A  =  [(X2+Y2)cos20-l]+4x2cos20 
B  =  +  {2Xsin0cos0[(X2+Y2lcos20+l]} 

C  =  [(X2+Y2)cos20-l]2+X2cos20(4-sin20)-sin20 

D  =  -  2Xsin0cos30(X2+Y2) 

2  4 

E  =  -  X  cos  0 

This  resulting  equation,  a  quartic  in  sin  a,  is  obviously  more  com¬ 
plicated  than  the  result  for  the  goniometer  system,  Eq.  (4.11).  Methods 

are  available  for  solving  quartic  polynomials,  'and  the  method  given  by 
27  * 

Milne -Thomson  was  employed  here.  The  resulting  solution  yields  four 

values  of  sin  Q,  where  a  is  (0-0  ),  the  bearing  error  angle.  There  are 

6 

eight  possible  values  of  a  which  satisfy  this  equation,  of  which  four 
are  maxima  and  minima.  The  other  four  solutions  proved  to  be  trivial. 

The  suspected  subsidiary  "nulls"  are  therefore  shown  to  exist.  Extensive 
calculations  of  theoretical  response  to  different  polarizations  and  verti 
cal  angles  of  arrival  were  made.  A  typical  example  of  the  maximum  and 
minimum  positions  is  shown  in  Fig.  4.1  and  the  maximum -to -minimum  ratios 
in  Fig.  4.2.  The  positions  of  the  maxima  and  minima  for  different  values 
of  Y  vs  the  goniometer  angle  are  shown  in  Fig. .4.1  for  a  single  ray  with 

* 

The  calculations  were  programed  on  a  digital  computer  and  are 
described  in  Appendix  C. 


result  is: 


where : 
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(  deg) 


NULL  BEARING  ERROR  (deg) 


polarization  X  =  1.0  and  6  ■  70  degrees.  Figure  4.2  shows  the  magnitude  of 
the  goniometer  null  ratio  and  the  ratio  of  the  goniometer -vertical  maxima 
and  minima  to  the  main  maximum.  The  correct  bearing  angle  of  the  gonio¬ 
meter-vertical  null  would  be  180  degrees.  A  scale  is  shown  to  indicate 
the  actual  null  bearing  error.  The  position  of  the  goniometer  null  as 
a  function  of  Y  is  shown,  using  the  same  null  bearing  error  scale.  The 
positions  of  the  maxima  and  minima  of  the  goniometer-vertical  are  shown 

by  the  0  scale  as  a  function  of  Y. 

S 

The  main  null  for  the  goniometer  and  vertical  has  a  bearing  error 
about  half  as  great  as  the  goniometer  null,  and  the  maximum -to -minimum 
ratio  is  improved  by  roughly  a  factor  of  ten  over  the  goniometer. 

The  only  general  exception  to  this  advantage  is  for  the  cases 
wherein  two  identical  nulls  appear.  These  are  closely  spaced  and  are  of 
equal  amplitude.  They  would  be  recognizable  by  their  symmetry,  however. 

The  true  bearing  would  be  njidway  between  them.  The  main  minimum  is  much 
deeper  than  the  subsidiary  ones  for  the  great  majority  of  cases,  and  so 
the  subsidiary  minima  are  not  likely  to  be  significant.  A  further  dis¬ 
cussion  is  given  on  page  143 .  ' 

It  was  noted,  in  the  section  on  the  goniometer,  that  Y  appears  only 
in  even  powers  while  X  appears  in  both,  odd  and  even  powers.  The  odd 
powers  occur  only  in  coefficients  of  odd-power  terms  in  sina.  The  effect 
of  a  change  in  the  sign  of  X  will  again,  as  in  the  case  of  the  goniometer, 
reverse  the  pattern  of  maximum  and  minimum  about  the  true  bearing  or,  in 
other  words,  change  the  sign  of  the  error.  The.  significance  of  this 
property  is  discussed  later,  on  page  l4l. 

D.  THE  EXPERIMENTAL  PROGRAM 

Direction  finding  measurements,  using  a  goniometer  system  of  the 
type  previously  described,  were  made  at  Stanford,  California,  and  at  a 
cite  near  Seattle,  Washington.  The  Stanford  station  utilized  the  antenna 
and  associated  equipment  built  for  Lissajous -pattern  studies  of  whistlers 
and  vlf  signals  and  for  whistler  studies  in  the  30-60  kc  range.  The 
Seattle  station  utilized  a  standard  two-channel  IGY  system.  Regular 
operation  of  the  two  stations  began  in  April'  i960  after  a  trial,  period 
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I  of  several  Veeks,  and  ffcntinfced  until  the  end  of  January  1961,  A  similar 
experiment  was  started  several  months  earlier  by  J.  M.  Watts  of  NBS  Boulder 
Laboratories.  Simultaneous  observations  among  the  three  stations  were 
obtained  for  a  considerable  time  by  mutual  agreement.  Hie  goniometer- 
vertical  system  was  added  to  the  Stanford  system  at  a  later  time. 

1.  EQUIPMENT.  Figure  4.3  shows  a  block  diagram  of  the  Stanford 
equipment.  The  outputs  of  the  goniometer  and  the  goniometer -vertical 
were  simultaneously  recorded  on  a  dual  channel  tape  recorder.'  The. 
operation  was  automatically  programed  to  record,  on  the  standard  IGY 
whistler  recording  schedule  of  2  minutes  per  hour.  A  calibration  signal. 

was  regularly  inserted  into  the  antenna  circuit  by  the  constant  current 

5  •  '  ' 

method  as  a  test  of  the  proper  adjustment.  The  tape  recordings  were 

played  back  through  a  spectrum  analyzer.  A  set  of  spectrograms  is  shown 
in  Figs.  4.4  and  4.5.  The  spectrograms  show  voltage,  measured  by  the 
relative  darkness  of  the  display,  as  a  function  of  frequency  and  time. 

These  spectrograms  illustrate  some  of  the  best  nulls  obtained  on  the 
goniometer-vertical  system.  . 

The  three  spectrograms  of  Fig.  4.4  show  the  signals  obtained  simul¬ 
taneously  from  the  Stanford  goniometer  and  goniometer-vertical  and  the 
Seattle  goniometer  system.  The  vertical  bars  on  the  spectrograms  pre 
the  atmospherics  from  lightning  discharges,  which  contain  broadband 
energy.  The  broken,  horizontal  lines  in  the  15-20  kc  region  of  all  the 
spectrograms  are  the  Morse  code  signals  from  navy  vlf  transmitters.  The 
one  in  Fig.  4.5  is  station  NPG,  18.6  kc,  located  near  Seattle,  Washington. 
The  recording  level  is  varied  on  spectrograms  so  that  the  desired  signal 
(in  this  case,  the  whistler)  appears  at  the  proper  amplitude  (or  inten¬ 
sity).  The  vlf  signals  are  therefore  too  strong,  in  this  spectrogram, 
to  show  nulls. 

The  whistler  shown  in  Fig.  4.4  is  a  good  example  of  the  multipath 
type  of  whistler  which  shows  several  traces  with  different  delays  which 
result  from  different  propagation  paths.  The  nulls  for  the.  different 
traces  are  easily  seen  to  occur  at  different  bearing  angles. 

Figure  4.5  shows  a  whistler  with  a  single  trace  but  a  very  regular 
null  pattern  which  varies  slightly  with  frequency. 
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FIG.  4. S.  BLOCK  DIAGRAM  OF  STANFORD  GONIOMETER  AND  GONIOMETER 
VERTICAL  SYSTEM. 


t _  CLOCKWISE  ROTATION  klSEG-l 

NORTH  REFERENCE 


FIG.  4.4.  SPECTROGRAMS  OF  WHISTLER  SHOWING  DIFFERENCE  BETWEEN 
NULLS  ON  THE  GONIOMETER  AND  GONIOMETER- VERTICAL  SYSTEMS  AT 
STANFORD  AND  THE  GONIOMETER  SYSTEM  AT  SEATTLE. 
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FIG.  4.5.  SPECTROGRAMS  OF  WHISTLER  WITH  WELL  DEFINED  NULLS  ON 
THE  GONIOMETER-VERTICAL  AND  POOR  NULLS  ON  THE  GONIOMETER  SYSTEM 
AT  STANFORD. 


A  variable -frequency  signal  from  an  oscillator  is  inserted,  with 
equal  amplitude  and  phase,  into  each  antenna  circuit,  .typical,  calibration 
spectrograms  made  of  this  signal  are  shown  in  Fig.  4.6.  The  position 
of  the  nulls  marks  the  coupling  of  equal  and  opposite  signals  from  the 

two  channels  into  the  goniometer  output,  where  the  goniometer  angle  0 

0  ^ 
is  45  .  The  series  of  marks  near  the  bottom  edge  of  the  spectrogram 

(520  cps  and  harmonics)  is  the  train  of  index  marks.  The  goniometer  is 

rotated  at  10  rps  by  a  synchronous  motor.  A  pulse  is  generated  by  a 

magnet  mounted  on  a  disc  rotating  with  the  goniometer,  and  passing  by  a 

coil.  The  voltage  induced  in  the  coil  is  amplified,  shaped  and  used  to 

trigger  the  one-shot  multivibrator  shown  in  the  block  diagram.  The 

resulting  pulse  turns  on  the  " gated -Hartley "  oscillator  and  produces  an 

index  mark  once  in  each  revolution.  Bearings  are  scaled  by  measuring 

the  position  of  a  null,  relative  to  the  index  marks,  and  converting  this 

to  angular  position. 

The  calibration  is  used  to  determine  the  correction  to  the  position 
of  the  index  mark  relative  to  the  goniometer  position  and  to  check  the 
equialization  of  the  channels.  Scaling  of  the  calibration  nulls  at 
different  frequencies  shows  that  the  shift  of  position  with  frequency  is 
not  significant  compared  to  scaling  error. 

Recorded  whistlers  were  analyzed  on  a  regular  basis,  and  a  systematic 
search  made  for  those  with  "nulls"  at  both  stations  so  that  bearings 
could  be  computed  and  the  position  of  the  exit-point  triangulated  from 
both  stations.  The  error  in  the  index -mark  position,  determined  from 
the  calibration  spectrogram,  was  applied  as  a  bearing  correction  factor 
during  scaling.  It  is  estimated  that  a  high-quality  spectrogram  can  be 
scaled  to  a  maximum  accuracy  of  +  1  degree  in  bearing;  the  average 
scaling  accuracy  is  estimated  to  be  of  the  order  of  +  2-4  degrees.  The 
resulting  error  in  position  depends  entirely  upon  the  position  of  the 
"target"  relative  to  the  station.  The  worst  cases  occur  when  the  "target" 
is  close  to  one  of  the  stations  or  the  line  Joining  them. 

Several  nulls  are  usually  obtained  on  one  whistler  trace.  Each  of 
these  gives  a  bearing  for  the  corresponding  frequency.  The  bearings 
obtained  for  the  whistlers  showing  nulls  at  both  stations  are  listed  in 
Table  I . 
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TABLE  I  MEASURED  NULL  BEARINGS  OF  LOCATED  WHISUErR 


Seattle 

Stan  ford 

Universal 
date  and 

Sonagran 

Frequency 

Searing 

Sonagran 

Frequency 

Bearing 

time 

nunSer 

nuBDer 

25  April  I960 

1A 

8.8 

1A 

9.0 

60 

1135/84 

7.0 

7.  4 

55 

5.0 

BiH 

5.  5 

62 

4.  2 

113 

4.  5 

54 

3.  5 

116 

4.  1 

62 

3.0 

140 

3.8 

56 

2.8 

122 

2.  5 

124 

ID 

4.0 

60 

_ 

3.8 

51 

Average 

Average 

58 

Trace  A(l) 

Trace  A(  1) 

3  July  1960 

1A 

6.8 

1A 

9.0 

144 

1435/89 

6.0 

8.0 

114 

5.0 

6.2 

140 

4.0 

5.  2 

127 

Trace  D(4) 

4.8 

104 

7.3 

4.0 

163 

6.2 

4 

Trace  D(4) 

4.8 

9 

6.0 

153 

IB 

Trace  D(4) 

5.  2 

147 

6,0 

20 

4.  5 

150 

5.0 

■  ■ 

4.  0 

162 

4.  5 

mm 

Average  A 

96 

Average  A 

110 

Average  D 

9 

Average  D 

151 

12  Aug  1960 

1A 

6.  2 

154 

1A 

6.9 

82 

1135/43 

5.  2 

134 

5.  5 

101 

4.2 

168 

4.  5 

83 

3.7 

157 

4.  0 

57 

3.2 

97 

2.  8 

112 

2.  5 

93 

. 

2.0 

76 

Average 

153 

Average 

88 

26  September 

1A 

5.8 

147 

IB 

8.  5 

124 

1960 

159 

6.  5 

105 

0435/43 

4.0 

139 

5.0 

147 

3  5 

144 

3.7 

139 

152 

3.  0 

125 

£-8 

142 

1.6 

131 

z  5 

117 

Average 

U3 

Average 

129 

0735/24 

IB 

7.  0 

177 

1C 

5  4 

154 

6.1 

174 

4.  2 

127 

5  5 

145 

3.  3 

118 

4.6 

117 

2.  5 

140 

4.  2 

133 

3.  5 

144 

Average 

mm 

Average 

135 

0835/29 

IB 

5.4 

■SB 

8.  4 

202 

4.5 

6.  5 

148 

4.0 

■fMf 

4.9 

148 

2.  5 

2.8 

202 

1.  5 

■vr®r 

Average 

132 

Average 

175 
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TABLE  I.  MEASURED  NULL  BEARINGS  OF  LOCATED  WHISTLERS  (Cont'd) 


Seattle 

Stanford 

Universal 
date  and 

Sonagran 

number 

Frequency 

Bearing 

Sonagram 

number 

Frequency 

Bearing 

tine 

0835/65 

IB 

5.8 

149 

10.  3 

5.0 

146 

7.7 

4.4 

142 

5.7 

3.  9 

96.9 

4.  3 

3.4 

155 

3.2 

2.9 

144 

2.7 

BBtf  £ 

2.6 

150 

2.0 

2.  3 

131 

2.  1 

127 

1.8 

123 

1.6 

144 

1.4 

117 

1  1 

Average 

135 

Average 

169 

0835/109 

IB 

8.  5 

175 

8.8 

154 

7.7 

108 

6.  2 

197 

5.5 

130 

4.6 

189 

3.  5 

153 

3.  5 

192 

3  2 

128 

2.8 

183 

2.8 

153 

2.  5 

134 

2.  2 

108 

2.0 

129 

1.7 

109 

Average 

133 

Average 

183 

1335/103 

1A 

4.5 

175- 

2B 

4.8 

207 

4.0 

182 

2.0 

224 

3.4 

201 

1.7 

189 

3  0 

185 

1.3 

198 

Average 

186 

Average 

204 

In  a  few  cases,  bearings  were  taken  from  two  different  spectrograms 
made  at  different  level  settings.  This  sometimes  produces  very  different 
bearings  at  nearly  the  same  frequencies.  This  is  associated  with  poor 
signal-to-noise  ratio  and  with  errors  such  as  flutter  in  the  tape  record¬ 
ing  and  reproducing  process.  The  variation  of  the  individual  bearings 
about  the  average  is  very  large,  as  is  expected  from  the  results  of  the 
calculations.  This  is  discussed  on  page  l4l. 

The  average  bearing  for  each  trace  is  used.  The  locations  are  then 
determined  by  solving  the  spherical  trigonometry  on  a  digital  computer. 


E.  EXPERIMENTAL  RESULTS 


During  the  measurements  program,  which  lasted  approximately  10 
months,  approximately  170  whistlers  were  recorded  at  Stanford;  130  of 
these  showed  scalable  nulls.  A  total  of  only  10  whistlers  yielded  nulls 
on  the  same  trace  that  were  scalable  at  both  Stanford  and  Seattle,  as 
Table  II  shows.  No  whistlers  were  obtained  which  showed  scalable  nulls 
at  Stanford  or  Seattle  and  also  at  the  Boulder  Station  run  by  Watts. 

The  positions  of  the  exit-point  locations  were  plotted  on  a  map  from 
which  the  magnetic  latitude  of  the  position  could  be  determined.  The 
map  is  shown  in  Fig.  4.7. 

The  magnetic  latitude  of  the  path. of  propagation  was  also  calculated 

by  analysis  of  the  spectrograms,  using  a  method  developed  by  Carpenter 

4  33 

and  Smith.  1  In  their  method  .the  nose  frequency  is  first  calculated 

and  then  used  to  estimate  the  latitude  of  the  path,  assuming  a  dipole 

4 

field  and  a  particular  form  for  the  electron  distribution.  Each  location 
shown  in  Fig.  4.7  is  connected  by  a  line  with  the  nearest  point  which  has 
the  magnetic  latitude  so  determined.  These  two  values  are  also  compared 
in  Table  II. 


TABLE  II.  COMPARISON  OF  MAGNETIC  LATITUDE  OF  WHISTLER  PATHS 
DETERMINED  BY  d-f  LOCATION  AND  NOSE  WHISTLER  THEORY. 

Ident. 

D-f 

Nose 

D-  f-nose 

map 

Date 

UT 

Trace 

e0 

e 

Mo 

no. 

(deg) 

(deg) 

.  (deg) 

1 

25  Apr  60 

1135/84 

1 

49.0 

45.5 

+3.5 

2 

3  Jul  60 

1436/29 

1 

37.0 

47.5 

-10.5 

3 

4 

46.0 

54.  5 

-8.5 

4 

12  Aug  60 

1135/43 

1 

44.0 

46.0 

-2.0 

5 

.  26  Sep  60 

0435/43 

■ 

30.0 

47.0 

-17.0 

6 

07  35/24 

■ 

29.0 

47.0,  46.5 

-17.8 

7 

0835/29 

K 

53.0 

46.3,  46.5 

+6.5 

8 

08  35/65 

’  ■ 

54.0 

46.  0 

+8.0 

0835/109 

52.0 

45.5 

+6 . 5 

1335/103 

■ 

40.0 

46.0 

-6.0 

Sub 

Average 
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The  differences  in  latitude  range  from  2.1®  to  17.8°.  The  average 
is  -3.7°  and  the  standard  deviation  is  9.2°.  .The  number  of  cases  is  too 
small  for  the  average  to  be  statistically  significant.  In  a  large  sample, 
it  could  be  significant  (see  page  149). 

A  further  result  of  the  experiment  is  shown  in  Fig.  4.8.  This  curve 
indicates  the  range  of  23  Stanford  ..goniometer  bearings' for  vlf  station  NSS 
(15.5  he)  at  Annapolis.  These  data  were  taken  on  4  different  days  spaced 
over  a  3-month  period.  The  scalings  were  made  without  prior  knowledge 

of  the  true  bearing.  The  average  bearing  for  this  set  is  73*4°,  and  the 

o  o 

standard  deviation  is  4.3  .  The  true  great  circle  bearing  is  72.99  • 

The  measurement  of  vlf  station  bearings  is  generally  more  accurate  than 
whistler  bearing  measurements,  since  the  field  strength. is  generally  ■ 
greater  and  always  more  constant.  The  results  illustrate  the  improvement 
which  may  result  from  averaging  bearings.  The  expected  measuring  error 
is  +  1  to  2  degrees.  The  remaining  error  can  be  attributed  to  polarization 
and  lateral  deviation  error.  Compared  with  whistlers,  the  error  would 
be  expected  to  be  smaller  since  vlf  stations  generate  vertically  polar¬ 
ized  waves. 


>•  • 
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v.  Diflctaaioi 


A.  THE  EFFECT  OF  DIFFERENT  ANTENNAS  ON  FIELD  STRENGTH 


The  difference  in  the  response  of  different  antennas  is  one  of  the 
most  striking  features  of  the  curves  of  field  strength.  It  is  evident 
that  the  loop  and  the  vertical  monopole  measure  different  fields,  even 
when  the  plane  of  the  loop  lies  in  the  plane  of  propagation.  An  examin¬ 
ation  of  the  results  shows  that  the  field  measured  by  the  vertical  antenna 
is  less  than  that  measured  by  the  loop  antenna  by  a  factor  of  approximately 
1  db  for  the  majority  of  cases;  Figure  3- ^  shows  one  point  where  the 
difference  is  roughly  6  db  in  the  opposite  direction.  Experimental 
measurements  at  Stanford  have  shown  a  discrepancy  of  this  magnitude  for 
vlf  station  NPG  on  18.6  kc. 

This  difference  in  the  fields  is  predicted  by  the  equations  for 
antenna  response.  The  difference  exists  even  when  the  incident  field  is 
vertically  polarized.  Consider  the  two  equations  for  the  response  to 
the  parallel  component  of  electric  field: 


3a 


?(n) 

'll 


=  2  E (1  +  R(n) 


(3.7) 


x_  =2  En  (.1  +  H,  )  cos  P 
3v  Hn  I'  n  n 

n 


(3.8) 


It  might  seem  at  first  that  the  magnitude  of  x^v  would  always  be  less 


than  the  magnitude'  of  x^a>  since  cos  P^ 


1.  Figure  5*1  shows  the  vector 
diagram  of  the  two  voltages  for  a  typical  case,  in  terms  of  the  incre¬ 
ments.  For  clarity,  only  the  first  5  components  are  shown.  For  each  n, 
the  component  of  both  fields  will  have  the  same  phase  angle  (which  measures 
time  phase),  since  x^  is  merely  x^a  multiplied  by  the  scalar  function 

cos • p  .  The  two  increments  for  each  n  will  thus  be  parallel,  with  x„ 
n  7  3a 

always  longer  than  x^.  The  magnitude  of  each  voltage  is  the  distance 
from  the  origin.  The  example  clearly  shows  that  either  voltage  may  have 
the  greater  magnitude. 
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The  difference  between  the  measured  field  strengths  of  vlf  signals 
on  vertical  and  loop  antennas  has  been  generally  explained  in  terms  of 
errors  in  the  effective  height  values,  site  errors  (ground  irregularities, 
etc.)/  or  polarization  errors.  However  the  difference  in  the  response 
caused  by  this  new  factor  should  be  considered. 

The  results  o?  the ' calculations  make  possible  an  evaluation  of 
different  antennas  for  particular  applications.  One  very  important 
question  that  has  often  arisen  in  whistler  work  is  the  relative  sensi¬ 
tivity,  to  whistler  signals,  of  a  vertical  loop  and  a  vertical  monopole. 

This  is  best  determined  from  Figs.  3.8  through  3.12,  where  field  strength 
is  plotted  vs  distance  for  5  kc  and  15.5kc.  It  is  evident  that,  the 
response  of  the  vertical  antenna  will  be  not  more  than  6  db  below  the 
normal  loop  response  for  D  greater  than  about  6o  km,  for  both  5  kc  and 
15.5  kc.  The  vertical  is  equally  good  for  the  great  majority  of- whistler 
signals,  therefore,  and  has  the  advantage  that  its  response  is  independent 
of  the  azimuth  of  arrival.  The  choice  can  then  be  b^sed  on  consider¬ 
ations  of  sensitivity  to  hum  and  other  interference.  Since  the  vertical 
cannot  be  nulled  for  him  or  atmospheric  noise,  the  interference  may  well 
be  the  deciding  factor. 

A  second  important  question  involves  the  orientation  of  a  loop 

antenna.  Figures  3-17  through  3*21  show  the  response  of  a  loop  vs  frequency 

for  different  values  of  0  .  It  is  evident  that  while  the  loop  orientation 

6 

affects  the  details  of  the  response  curves  as  a  function  of  frequency 
or  distance,,  the  general  behavior  is  the  same  for  all  orientations.  An 
examination  of  Figs.  3.9  and  3. 10  shows  that  the  amplitude  of  the  signal 
in  the  abnormal  (EW)  loop  varies  considerably,  relative  to  that  in  the 
normal  (NS)  loop.  It  may  be  seen  from  Fig.  3.9/  the.  curve  for  5  kc, . 
that  the  ratio  of  the  NS  to  EW  loops  is  not  greater  than  6  db  for  any 
distance.  However,  Fig.  3. 10  for  15.5  kc  shows  that  the  ratio  is  gen-  . 
erally  greater  than  6  db  for  distances  greater  than  600  km.  This  indi¬ 
cates  that  the  loop  orientation  is  important  in  the  upper  part  of  the 
vlf  spectrum.  There  is  some  experimental  evidence  to  support  this  theory. 

At  Byrd  Station,  Antarctica,  no  "whistler  mode"  signals  were  received 
from  station  IBS  at  15-5  kc. until  the  loop  antenna  was  oriented  so  that 
it  pointed  at  the  area  where  the  ducts  could  be  expected  to  terminate. 
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Vhistler  and  vlf  activity  vLi eh  aectfi’S  at  the  lew  frequencies  was  normal 

for  this  location. 

B.  FIELD  STRENGTH  VARIATIONS 

Rapid  variations  of  field  strength  with  frequency  are  characteristic 
of  all  but  the  0.1  km  calculations.  The  shape  of  the  curves  for  D  *  0.1 
km,  shown  in  Figs.  3«1  and  3*2,  is  of  special  interest.  The  conditions 
are,  of  course,  very  close  to  vertical  incidence.  For  this  case,  the 
earth-ionosphere  waveguide  appears  to  the  transmitted  wave  as  a  trans¬ 
mission  line  with  the  earth  as  a  termination  at  the  end.  If  the  ground 
is  relatively  good,  as  in  sea  water,  the  physical  situation  closely 
resembles  the  shorted  transmission  line,  and  a  loop  antenna  on  the  ground 
would  be  measuring  the  equivalent  of  the  magnetic  field  at  the  short. 

This  analogy  is  demonstrated  by  expressing  the  magnetic  field  strength 
at  the  earth  in  terms  of  ordinary  transmission  line  equations. 

The  transmission  line  equations  represent  the  effect  of  the  multiple 
reflections,  present  in  the  line  as  a  change  in  the  impedance  at  the 
input  to  the  transmission  line.  By  contrast,  the  ray -theory  method 
requires  the  determination  of  the  effect  of  each  ray.  For  this  reason, 
the  ray -theory  calculation  is  very  sensitive  to  small  changes  in  the 
model.  There  are  differences  in  the  two  situations,  due  to  the  assumed 
l/R. attenuation  of  the  field  and  the  imperfect  ground  in  the  ray -theory 
case'.,  However,  t'he  situations  are  still  very  similar.  A  typical  curve 
(based  on  ray  theory)  for  the  summer-night  model  and  sea  water  ground 
is  shown  in  Fig.  3*2.  Figure  5.2  shows  a  portion  of  the  shorted  trans¬ 
mission  line  calculation  for  a  perfect  conductor  at  the  short,  with  the 
summer-night  model  ionosphere  constants.  The  similarity  in  the  curves 
is  obvious,  and  shows  the  analogy  between  the  two  viewpoints  in  this 
case. 

Variations  in  the  field  strength  similar  to  those  shown  in  Figs. 

3.2  and  5*2  have  been  seen  in  at  least  one  whistler.  This  whistler 
shows  several  fairly  sharp  nulls  in  its  pattern  at  anhannonic  frequencies 
which  are  spaced  by  0.8  -  1.2  kc  in  the  frequency  range  below  5  he. 
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This  could,  be  explained  by  the  shape  of  the  curves  of  field  strength 
for  vertical  incidence.  .The  minima  in  the  field  strength  curves  tend 
to  be  rather  broad,  however,  while  the  nulls  in  the  whistler  are  fairly 
sharp.  The  nulls  could  also  be  produced  by  beating  of  adjacent  rays  of 
different  frequencies.  It  was  shown  in  Section  II -B  that  the  difference 

at  vertical  incidence  for  a  frequency  of  10  kc  and  a  typical  dispersion 

•  l/2 

of  106  (sec)  '  is  about  13  cps.  This  difference  is  reduced  to  about 

2  ops  at  3  kc.  The  nuil  separation  is  approximately  correct  for  the 
•variation  due  to  the  transmission  line  effect  but  the  period  is  probably 
too  short  for  beating  effects  at  this  frequency. 

The  beating  effect  alone  would  not  produce  very  large  variations, 
since  the  ray  for  n  =  1  is  weaker  than  that  for  n  =  0  by  at  least  a 
factor  of  three  due  to  the  increased  propagation  distance.  It  is  there¬ 
fore  more  likely  that  the  cause  is  primarily  the  transmission -line 
effect.  The  appearance  of  these  nulls  is  an  easy  way  to  identify  short- 
distance  propagation  from  routine  spectrograms. 

The  behavior  of  the  polarization  and  the  crossed-loop  Lissajous 
patterns  is  Also  unique  for  near -vertical  incidence,  and  will  be  dis¬ 
cussed  on  page  133. 

The  curves  for  greater  distances  show  the  characteristic  rapid 
variations  with  frequency  due  to  interference  of  different  rays.  A  com¬ 
parison  of  the  curves  in  Figs.  3-^  and  3*7  shows  that  this  effect  is 
not  noticeably  greater  at  2000  km  than  at  500  km.  The  most  evident  effect 

is  that  of  changes  in  ionization  density.  The  curve  in  Fig.  3«15  is  for 
'  k 

cd  '=  5  x  10  .,  a  factor  of  10  lower  than  the  summer-night  model  for  the 
same  distance  (500  km)  as  in  Fig,  3*3*  This  corresponds  to  a  reduction 
in  f  of  Vio)  or  of  10  in  the  ion  density.  It  is  immediately  evident 
that  the  rapid  variations  are  much  smaller  and  the  field  about  6  db 
weaker  than  for  the  summer-night  model.  The  slope  and  position  of  the 
maxima  and  minima  are  also  different.  These  changes  are  reasonable 
since  the  transmission  coefficients  decrease  slightly  with  decreasing 
oor,  while  the  reflection  coefficients  (|R||-  and  ^R^  decrease  fairly 
rapidly  with  decreasing  cnr  as.  shown  in  Figs.  3*38  through  3 - 55-  The 
magnitude  of  the  cross -coupling  reflection  coefficients  u  and  ^  R  || 
change  very  slowly  with  j  and  their  magnitude  is  only  of  the  order  of 
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one -fourth  that  of  the  other  two  coefficients .  thus,  tht  amplitude  of 
the  rays  decreases  much  more  rapidly  with  number  of  reflections,  and  a 
smaller  number  of  rays  contributes  to  the  signal.  This  makes  the  inter¬ 
ference  patterns  less,  complex. 

The  corresponding  curve  for  an  increased  density,  a>r  b  5  x  10  ,  is 
shown  in  Fig.  3*16*  As.  would  be  expected,  the  pattern  is  more  complicated 
and  the  variations  dre  more  rapid.  • 

Rapid-  variations  in  field  strength  are  seen  in ' experimental  data. 
Figure- 5*3  shows  the  signal  voltage  vs  time  for  the  parts  of  a  whistler 
of  4  March  1956  at  O4o4  UT.  An  examination  of  a  typical  section  shows 
that  the  field  strength  variation  takes  place  in  a  minimum  period  of 
about  3  ms.  The  amplitude  fades  from  the  peak  values  down  to  noise 
level  during  these  variati :ns,  covering  a  range  of  about  20  db.  Fluctu¬ 
ations  of  this  magnitude  are  very  rare  in  the  calculated  results,  and 
do  not  occur  with  such  short  periods. 

It  appears  that  these  variations  are  due  to  some  factor  other  than 
the  interference  between  rays.  The  period  of  3  ms  corresponds  to  a 
frequency  difference  of  about  300  cps  if  this  were  caused  by  "beating" 
of  signals  of  different  frequency.  The  whistler  spectrogram.  Fig.  5*4, 
shows  that  the  whistler  is  about  1  kc  wide  at  this  frequency.  It  Is 
quite  possible  that  this  variation  is  due  to  the  interference  of  differ¬ 
ent  frequency  components.  Such  rapid  variation  cannot  be  seen  on  the 
observed  crossed-loop  Lissajous  patterns  probably  because  of  the  low 
sampling  rate. 

1.  THE'  EFFECT  OF  GROUND.  The  electrical  characteristics  of  the 
ground  have  little  effect  on  the  detailed  variations  in  the  calculated 
data.  The  main  effect  is,  of  course,  the  relative  attenuation.  Sea 
water  adds  about  0.2  db  to  the  signal  at  500  cycles,  and  1  db  at  20  kc, 
over  poor  ground  at  a  distance  of  0.1  km.  At  a  distance  of  2000  km, 
the  increase  is  that  0.3  db  at  500  cycles  and  15  db  at  20  kc.  These 
values  are  typical  in  that  the  relative  attenuation  is  low  at  low  fre¬ 
quencies.  One  feature  is  immediately  evident:  that  the  curves  for  poor 
ground  show  a  definite  slope  toward  greater  loss  as  the  frequency  is 
increased.  ThiB  does  not  appear  on  the  curves  for  sea  water.  The  effect 
is  especially  evident  on  the  curves  for  long  distances,  such  as  Fig.  3*6 
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FIG.  S. 3,  WAVEFORMS  OF  STANFORD  WHISTLER  OF  4  MARCH  1956,  0404  UT 


SPECTROGRAM  OF  STANFORD  WHISTLER  OF  4  MARCH  1956,  0404  UT 


for  2000  km.  This  indicates  that  there  is  considerable  advantage  in 
locating  middle-  or  low -latitude  stations  near  a  sea  coast,  especially 
when  high  frequency  data,  such  as  nose  whistlers  or  vlf  station  echoes 
are  desired.  This  finding  agrees  in  general  with  the  experimental  data, 
indicating  that  the  majority  of  whistler  sources  are  generally  over  the 
ocean,  to  the  west  of  Stanford. 

2.  FREQUENCY  EFFECTS.  The  curves  of  field  strength  vs  frequency 
(Figs.  3-1  through  3*21)  do  not  show  any  tendency  to  exhibit  a  peak  in 
any  part  of  the  frequency  range.  Figures  3.38  through  3.59  show  that 
no  peak  exists  in  the  transmission  or  reflection  coefficients.  The 
rapidly  varying  nature  of  the  curves  of  antenna  voltage  vs  frequency 
make  it  difficult  to  define  a  figure  for  the  loss  at  a  particular  fre¬ 
quency.  It  is  evident  that  the  loss  always  goes  down  near  the  low  end 
of  the  frequency  range,  around  0.5  kc.  A  fairly  deep  null  appears  above 
this  frequency  in  the  curve,  and  the  position  and  depth  depend  upon  the 
distance,  the  ionosphere  parameters,  and  the  antenna  orientation  (for 

a  loop).  It  is  therefore  very  difficult  to  make  general  statements  about 
transmission  as  a  function  of  frequency.  The  experimental  data  on 
whistlers  generally  show  a  drop  off  in  intensity  at  the  low  end  of  the 
spectrum  near  the  "tweek"  cutoff,  where  the  earth-ionosphere  waveguide 
becomes  a  half  wave  length  in  height.  The  interpretation  of  this  fact 
is  complicated  by  the  energy  spectrum  of  the  source,  which  is  known  to 
drop  off  at  low  frequencies.  This  experimental  effect  could  be  produced 
by  the  null  in  the  curves,  together  with  the  reduction- in  the  excitation 
at  low  frequencies.  The  rise  at  lower  frequencies  could  be  masked  by  the 
equipment,  which  generally  begins  to  roll  off  Just  above  0.5  kc,  and  by 
60  cycle  harmonic  components. 

3.  EFFECTIVE  AREA.  The  "effective  area"  of  a  whistler  is  defined 

as  that  area  over  which  a  whistler  may  be  detected.  This  area  is  assumed 

to  be  circular,  and  the  value  of  1000  km  has  often  been  quoted  for  the 

12  3 

radius  as  determined  from  experimental  data.  }  *  If  a  reasonable  value 
of  20  db  is  chosen  for  the  range  of  detectable  whistler  amplitudes,  the 
curve's  of  Figs.,  3.8  through  3. 10  show' that  this  change  in  amplitude  occurs, 
in  approximately  1000  km. 
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C,  SBff  CROSSED -LOOP  USSAJQUS  HlTTEfcBS 


One  of  the  results  of  the  experimental  studies  was  a  number  of 
recordings  of  the  crossed-loop  Lissajous  patterns  as  a  function  of  fre¬ 
quency.  A  particularly  good  example  of  this  is  the  whistler  of  4  March 
1956,  shown  in  Fig.  5  >5.  This  shows  the  waveforms  of  the  signal  on 
two  beams,  and  the  resulting  Lissajous  pattern  on  the  third  beam  of  a 
multibeam  oscilloscope.  The  beams  are  gated  on  only  during  the  200 -n 
sec  sweep.  The  sweep  repetition  rate  is  25  per  second,  which  corresponds 
to  0.04-second  intervals  between  sweeps.  The  whistlers  were  first  tape- 
recorded  on  a  carefully  equalized  tape  recorder..  They  were  then  playedi 
back  and  filmed. 

The  upper  strip  shows  a  10  kc  calibration  signal  introduced  into 
the  system.  The  effects  of  tape  flutter  are  evident  in  the  changes  in 
the  axis  ratio  of  the  ellipses  and  in  the  amplitude  of  the  waveforms. 

It  is  also  evident  that  the  flutter  effects  are  small  in  comparison  to 
the  total  changes.  The  striking  thing  about  these  records  is  the  fairly 
regular  change  in  the  axis  ratio  and  angle  of  the  ellipse.  Figure  5*6 
shows  the  characteristics  of  the  ellipse  for  different  sweeps.  The  axis 
ratio  varies  from  almost  zero  to  about  0.8,  and  the  angle  varies  from 
zero  to  almost  180  degrees.  Periods  in  the  range  of  0.04  -  0.13  seconds 
are  evident  when  this  and  other  examples  are  examined.  The  Isolated 
cases  of  opposite  rotation  are  also  of  interest,  as  discussed  on  page  137. 

The  calculated  real  and  imaginary  parts  of  vR,  the  voltage  ratio 
of  the  MS  and  EW  loops,  are  shown  in  Figs.  3.22  through  3*30  for  various 
distances,  ground  parameters  and  ionosphere  conditions.  It  is  evident 
that  the  behavior  is  very  complicated,  except  for  very  short  distances. 

At  a  distance  of  0.1  km,  the  polarization  is  circular  and  the  loop 
voltages  are  essentially  equal  and  90  degrees  out  of  phase.  The  imagin¬ 
ary  part  of  v_  is  then  always  equal  to  unity,  and  the  real  part  is  very 
small.  When  the  distance  is  increased,  the  ratio  departs  from  this  con¬ 
stant  value  by  only  about  %  at  20  km.  The  difference  increases  to 
about  5C$  at  100  km.  The  behavior  Is  characterized  by  rapid  changes, 
as  a  function  of  frequency,  for  greater  distances. 
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FIG.  5.5.  WAVEFORMS  AND  LISSAJOUS  PATTERNS  OF  STANFORD  WHISTLER  OF 
4  MARCH  1956  0404  UT.  ~ 


FIG.  5.6.  CHARACTERISTIC  OF  LISSAJOUS  PATTERN  ELLIPSE  OF  STANFORD  WHISTLER 
OF  4  MARCH  1956  0404  UT. 


Ltt  ve  consider  the  behavior  of  v  la  terms  of  the  Lissejons  pattern. 

A 

Each  time  the  real  part  changes  sign  from  +  to  -  it  indicates'  a  shift 
of  the  ellipse  from  the  first  and  third  quadrants  to  the  second  and 
fourth.  Each  change  in  sign  of  the  imaginary  part  indicates  a  change 
in  the  direction  of  rotation,  and  would  not  be  visible  in  the  photo¬ 
graphs,  but  can  be  determined  from  the  waveforms.  Appendix  D  gives 
details  and  curves  for  converting  the  real  and  imaginary  (XandY)  no¬ 
tation  to  the  axis  ratio,  e,  and  tilt  angle,  0,  of  the  major  azis  for 
any  ellipse.  Essentially  linear  polarization  occurs  where  Y>  is  very 
small  compared  to  X.  Circular  polarization  occurs  where  X  is  very  small 
and  Y  is.  near  unity.  Both  these  cases  occur  in  the  calculated  curves 
and  in  the  experimental  data. 

A  linear  pattern  occurs  in  Fig.  3-23  for  2.2  kc,  and  an  essentially 
circular  one  for  about  2.35  he.  The  only  features  of  the  Lissajous 
pattern  that  can  be  called  characteristic  for  any  parameter  are  those 
for  short  distances  aiid  for  low  values  of  such  as  in  Fig.  3.29.  The 
latter  occurs  for  values  of  ^  too  low  for  actual  ionosphere  conditions. 

The  curves  for  D  =  2000  km  are  considerably  smoother  than  the  500- 
km  and  1000-km  curves,  especially  for  the  higher  frequencies.  The  values 
sometimes  fluctuate  considerably,  and  sometimes  form  a  fairly  smooth 
curve  over  part  of  the  frequency  range.  Consider  Fig.  3«25*  The  curves 
are  fairly  smooth  in  the  range  around  10  kc.  The  pattern  at  10.22  kc 
corresponds  to  a  roughly  linear  polarization  with  0  =  -  58°  and  €  =  0.12. 
The  pattern  then  changes  smoothly,  with  frequency,  with  X  passing  through 
0  at  f  =  10.5  he.  This  pattern  would  be  an  ellipse  with  the  axis  horiz¬ 
ontal.  The  pattern  remains  approximately  the  same  in  a  range  around 
12  kc,  where  .0  =  -  60°  and  e  =  0.40.  It  then  changes  smoothly  again  until 
f  =  14  kc,  where  the  pattern  is  circular.  It  then  passes  through  a  very 
small  linear  pattern  at  f  =  14.5  he,  and  then  starts  more  rapid  changes. 

At  15.4  kc  the  pattern  Is  tilted  to  0  =  85°  and  e  =  0.3*  By  15.5  kc, 
it  has  changed  to  0  =  +  79°  and  €  =  0.025.  This  yields  an  ellipse  with 
its  axis  almost  vertical  changing  into  a  similarly  oriented,  almost 
linear  pattern. 
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It  is  evident  from  the  above  description,  that  the  variations  In 
the  calculated  patterns  very  much  resemble  those  seen  in  the  experimental 
data,  as  in  Fig.  5-4,  which  shows  the  Lissajous  patterns  of  the  4  March 
1956  whistler.  The  variations  are  sometimes  fairly  periodic,  and  a 
period  which  occurs  in  a  frequency  interval  1  to  1.5  he  can  be  identified. 
Let  us  determine  the  time  interval  involved. 

The  dispersion  equation  for  a  whistler  is:^" 

D  =  t  VF 

where  t  =  elapsed  time  from  the  occurrence  of  the  initiating  pulse  in  sec. 
f  =  the  instantaneous  frequency  in  cps 

1/2 

D  =  the  dispersion  of  the  whistler  in  sec  ' 

The  variations  of  the  4  March  1956  whistler,  shown  in  Fig.  5*4  and 

occurring  at  about  1.7  sec,  correspond  to  the  second  trace  of  the  whistler 

1/  2 

with  a  dispersion  =  175  sec  '  and  a  frequency  of  about  10  kc.  When 
values  of  the  1  and  1.5  frequency  change  are  substituted  into  Eq.  (5.1), 
the  values  of  At  are  approximately  0.10  and  0.15  second.  This  is  in  the 
upper  range  of  the  experimental  period  of  0.04  -  0.13  second  determined 
for  the  whistlers. 

The  isolated  points  in  Fig*  5.6,  which  show  a  different  direction 
of  rotation  of  the  ellipse  from  those  on  either  side^  are  also  of  inter¬ 
est,  as  the  same  sort  of  behavior  is  evident  at  times  in  the  calculated 
curves.  The  reversal  of  rotation  occurs  each  time  the  imaginary  part 
of  v-j  changes  sign.  Examples  of  the  isolated  points  occur  at  11.6  and 
15.4  kc  in  Fig.  3.23.  Other  cases  of  reversal  were  evident  on  curves 
not  reproduced  here.  The  only  feature  of  the  experimental  data  on  the 
Lissajous  patterns  that  is  not  generally  reproduced  is  that  which 
occurred  in  the  13  September  1955  whistler.  No  reversal  in  the  direction 
of  rotation  is  apparent  for  that  one  whistler.  Some  reversal  is  generally 
expected.  ,  The  only  exception  which  would  be  expected  in  general  is  for  . 
whistlers  close  to  vertical  incidence,  for  which  the  pattern  is  essen¬ 
tially  constant.  This  may  mean  that  the  distance  was  in  the  intermediate 
range,  where  not  enough  variation  in  the  imaginary  part  of  v^  is  present,  so 
that  its  sign  does  not  change. 


Another  feature  of  Interact  is  the  algebraic  signs  of  the  reel  and 

imaginary  parts  of  v_.  It  might  seem  reasonable  at  first  that  the 

JR 

imaginary  part  should  have  the  sign  corresponding  to  the  direction  of 
rotation  of  the  signal  in  the  ionosphere.  An  examination  of  the  curves 
shows  that  the  calculated  direction  of -rotation  of  the  crossed-loop 
Lissajous  pattern  is  as  often  of  one  sign  as  the  other  when  the  entire 
frequency  range  is  considered. 

Let  us  examine  the  meaning  of  the  sign  of  the  real  part  of  vR.  It 
indicates  which  set  of  quadrants  the  ellipse  occupies.  Changing  0  has 

o 

the  effect  of  rotating  the  ellipse  relative  to  the  coordinates.  There¬ 
fore,  the  sign  and  range  of  the  real  part  depends  upon  the  direction  of 

arrival  relative  to  0  .  For  a  particular  direction,  then,  it  is  possible 

S 

that  the  value  of  the  real  part  may  have  an  average  that  is  positive, 
negative,  or  zero. 

The  behavior  of  vR  is  of  interest  also  because  of  its  effect  on  the 
error  in  goniometer  systems .  It  was  noted  on  page  10 ^  and  page  HO  that 
the  sign  of  the  theoretical  error  in  the  goniometer  and  goniometer- 
vertical  system  was  determined  by  the  sign  of  X,  the  real  part  of  the 
complex  polarization.  In  the  actual  situation,  the  polarization  is  not 
defined  for  those  conditions  wherein  there  is  no  equivalent  plane  wave. 

If  the  meaning  of  the  goniometer  maximum/minimum  ratio  is  considered,  it 
is  readily  seen  that  the  null  occurs  at  the  angular  position  of  the  minor 
axis  of  the  crossed-loop  Lissajous  pattern  ellipse.  The  Lissajous  pattern 
ellipse  is  actually  a  representation  of  the  goniometer  voltage  as  a 
function  of  goniometer  angle.  There  should  be  a  relation  between  vR  and 
the  bearing  error.  This  relation  will  be  examined  in  the  next  section. 


GONIOMETER  AND  GONIOMETER -VERTICAL  NULL  CHARACTERISTICS 


Now  let  us  consider  the  goniometer  null  bearings  as  calculated. 
Figures  5-7  and  5.8  show  typical  goniometer  voltage  curves.  Figure  3.22 
shows  the  teal  suid  imaginary  parts  of  vR  vs  frequency  for  the  same  con¬ 
ditions  as  in  Figs.  5.7  and  5.8.  A  comparison  of  Figs.  5.7  and  5.8  with 
Fig.  3-22  shows  that  the  bearing  error  changes  sign  when  the  reed,  part 

of  v_  changes  sign.  This  is  to  be  expected  because  of  the  an £dogy 
R 
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FIG.  5.8.  CALCULATED  MAGNITUDE  OF  GONIOMETER  VOLTAGE  VS  4 FOR 
f  -  3.0,  4.0,  S.O  kc,  SUMMER-NIGHT  MODEL  IONOSPHERE,  POOR  GROUND, 
D  -  500  k m . 
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between  and  the  Q  uaad  In  the  eingle-rny  goniometer  aquation* .  The 
sign  of  the  bearing  error  changes  whenever  the  sign  of  X  In  the  equation* 
for  bearing  error  is  changed.  The  behavior  of  v^  is  therefore  significant, 
since  the  sign  of  the  real  part  of  v^  determines  the  sign  of  the  bearing 
error.  It  was  previously  noted  that  the  average  value  of  the  real  part 
of  Vp  could. tend  toward  any  vaiue.  The  average  depends  upon  the  distance, 
frequency,  and  all  the  waveguide  parameters.  The  average  error  could 
therefore  be  improved  if  bearings  of  a  single  exit  point  are  obtained 
over  a  wide  range  of  conditions  and  frequencies. 

Sometimes  several  whistlers  are  received  over  a  period  of  several 
minutes  to  several  hours.  Each  has  the  same  characteristic  frequency - 
vs-time  variation.  There  is  reason  to  believe  that  these  arise  from  a 
single  propagation  path  and  from  a  single  exit  point.  The  goniometer 
nulls  would  occur  at  different  frequencies  on  each  whistler,  since  the 
rotation  is  random  with  respect  to  the  time  of  occurrence.  The  averaging 
of  bearings  from  these  whistlers  should  reduce  the  error.  This  is  not 
true  in  the  case  of  the  series  of  whistlers  received  on  26  September 
i960.  The  exit-point  is  apparently  near  the  line  Joining  the  two  d-f 
stations  and  the  d-f  error  is  very  large. 

Figure.  5*9  shows  the  deviation  of  the  bearing  of  a  whistler 
(Stanford  10  April  i960,  0941/ 46UT),  obtained  at  different  frequencies, 
from  the  average  for  that  whistler.  The  bearing  errors  computed  for 
conditions  of  Figs.  5*7  and  5*8  are  also  shown  in  Fig.  5*9*  A  wide  vari¬ 
ation  in  apparent  direction  occurs  in  both  the  measured  and  the  calcu¬ 
lated  points.  The  standard,  deviation  of  the  measured  points  is  approxi¬ 
mately  20  degrees,  while  that  of  the  calculated  points  is- approximately 
10  degrees.  This  difference  is  too  great  to  be  explained  by  measuring 
error  which  was  estimated  on  page  115  to  be  +  2  -  4  degrees.  The  bearings 
on  another  whistler  recorded  an  hour  later  had  a  standard  deviation  of 

15* 

The  calculated  results  show  that  the  bearing  errors  are  a  very  com-  . 
plicated  function  of  all  the  variables.  The  magnitude  of  the  error 
varies  widely  under  all  conditions.  The  large  variations  of  the  measured 
bearings  is  what  would  be  expected  for  any  conditions  other  than  nearly 
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FIG.  5.9.  CALCULATED  BEARING  ERRORS  FOR  SUMMER-NIGHT  MODEL  IONOSPHERE,  POOR  GROUND, 
D  -  500  km,  AND  MEASURED  BEARING  ERRORS  FOR  WHISTLER  OF  10  APRIL  1960  AT  0941:46  UT, 
VS  FREQUENCY. 


vertical  incidence.  Uhls  follows  from  the  fact  that  the  behavior  is 
related  to  ve  and  0  ,  which  vary  widely  with  distance,  ionosphere 
parameters,  etc. 

The  next  step  is  to  compare  some  of  the  calculated  bearing  errors 
and  null  ratios  for  the  goniometer  system  with  those  calculated  for  the 
goniometer -vertical  system.  Some  typical  curves  for  the  latter  are 
plotted  in  Fig.  5*10  for  the  same  conditions  as  in  Figs.  5.7  and  5.8  for 
the  goniometer.  These  values  are  shown  in  Fig.  5*10  and  compared  in 
Table  III..  It  is  evident  from  Fig.  5*9  and  Table  III  that  a  result 
not  forseen  from  the  simple  plane -wave  theory  of  the  goniometer -vertical 
system  has  appeared  here.  The  calculated  errors  for  the  goniometer- 
vertical  are  generally  considerably  larger  than  those  for  the  goniometer 
alone.  When  the  presence  of .  multiple  rays  is  considered,  this  result  is 
reasonable.  The  main  null  of  the  goniometer -vertical  depends  upon  the 
difference  between  the  normal  loop  and  the  vertical  monopole  antenna 
voltages.  When  multiple  rays  are  present,  these  voltages  do  not  have 
a  fixed  relationship,  in  contrast  to  the  single -ray  situation.  Since 
the  null  depends  upon  the  difference  between  these  voltages,  the  error 
would  be  expected  to  vary  considerably.  Note  that  the  computed  average 
bearing  error  is  still  of  the  order  of  the  average  scaling  error  (.2  -  4°). 

It  will  be  noted,  however,  that  the  nulls  are  generally  much  deeper 
for  the  goniometer-vertical  than  for  the  goniometer  alone.  The  nulls 
for  the  goniometer  are  always  symmetrical  and  narrow.  The  nulls  for  the 
goniometer-vertical  are  generally  broad  and  asymmetrical  near  the  bottom, 
due  to  the  multiple  nature  of  the  nulls. 

The  records  must  be  analyzed  on  a  spectrum  analyzer.  In  normal 
operation  it  is  the  practice  to  reduce  the  signal  to  the  level  at  which 
only  the  peaks  of  the  desired  signal  are  recorded.  Interfering  signals 
are  reduced  as  much  as  possible  in  this  manner.  Spectrum  analyzers 
always  have  a  limited  dynamic  range.  A  value  of  6  db  is  typical.  There¬ 
fore,  only  the  upper  part  of  the  pattern  is  generally  evident  on  the 
spectrogram.  If  the  measurement  is  made  on  the  upper  part,  say  the 
upper  6-db  range,  the  asymmetry  is  almost  removed.  The  error  measured 
in  this  way  is  generally  of  the  same  order  as  the  goniometer  null  error. 
The  fourth  column  in  Table  III  indicates  this  fact. 
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FIG.  S.XO.  CALCULATED  MAGNITUDE  OF  GONIOMETER- VERTICAL  VOLTAGE  VS  <£  FOR 
f  *  1.0,  l.S,  2.0,  2.5,  3^0,  4.0  kc  AND  5.0  kc,  SUMMER-NIGHT  MODEL  IONOSPHERE. 
POOR  GROUND,  D«  500  k«.  5 


TABLE  III.  COMPARISON  OF  BKARING  ERROR  CALCULATED  FOB  GONIOMETER 
AND  GON IOMETER -VESTI CAL  SYSTEMS  FROM  SUMMER-NIGHT  MODEL  IONOSPHERE, 
POOR  GROUND,  D  -  S00  km. 


Bearing 

error 

Max/nin  ratio 

f  (kc) 

Goniometer 

Goniometer- 

vertical 

-6  db  goniometer- 
vertical 

Goniometer 

Goniometer- 

vertical 

1.0 

-3±1 

-9±2 

-  3.  5±2 

16 

154 

1.5 

-3±1 

-8±2 

-3.5±2. 

19.5 

500 

2.0 

-18±1 

-38±2 

- 18. 5±2  • 

9.5 

15 

2.5 

+  30±2 

+28±3 

+  1.513 

3 

74 

3.0 

+13±2 

+40±2 

+  13,  5±2 

11 

55.3 

4.0 

-7±2 

-20*2 

8. 0±2 

.7 

67.7 

5.0 

-7±2 

-15±3 

+  5. 0±3 

3 

23.7 

Average 

+0.7 

wm 

-1.9 

E.  THE  ANGLE  0  AND  LATERAL  DEVIATION 
a 

The  angle  0  is  the  hearing  angle  of  the  equivalent  plane  wave  as 
defined  on  page  34.  It  depends  on  the  relative  amplitudes  and  phases 
of  the  voltages  from  the  three  antennas.  In  general,  it  differs  from 
the  true  hearing  of  the  individual  rays  which  are  assumed  to  arrive  from 
the  great  circle  hearing.  This  phenomenon,  which  we  shall  call  apparent 
lateral  deviation,  is  separate  from  the  polarization  error,  which  is 
accounted  for  in  Eqs.  (2.81)  through  (2.88). 

The  phenomenon  of  lateral  deviation  has  been  observed  over  a  wide 

07  oft 

range  of  frequencies  for  several  years.  ’  It  is  characterized  as 
the  arrival  of  a  signal  from  a  bearing  apparently  different  from  the 
great  circle  bearing.  At  vlf,  this  phenomenon  is  associated  with  the  . 
difference  in  the  sums  of  the  multiple  rays  that  make  up  the  total  field 
at  different  antennas.  The  non -plane -wave  character  of  the  received 
field  for  some  conditions  is  a  related  result.  Apparent  lateral  devia¬ 
tion  is  the  result  of  the  presence  of  multiple  rays.  Polarization  error 
can  occur  from  a  single  ray. 
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The  calculated  angle  0.  Is  platted  as  a  function  of  frequency  for 
different  distances,  Ionosphere  conditions,  and  ground  constants  In 

Figs.  3.35  through  3. 45.  It  is  readily  apparent  that  the  value  of  0. 
is  highly  variable.  It  is  also  apparent  that  its  average  value  tends 
to  zero  (i.e.,  zero  average  error)  for  a  large  range  of  frequencies, 
distances  and  waveguide  parameters.  The  significance’ of  this  is  again 
that  an  average  over  a. large  number  of  measurements  under  varied  con¬ 
ditions  will  approach  the  true  bearing,  even  though  individual  errors 
may  be  large.  This  behavior  is  similar  to  the  effect  noted  for  gonio¬ 
meter  bearing  error.  The  errors  due  to  apparent  lateral  deviation  are 
sometimes  difficult  to  separate  from  polarization  error.  The  best 
example  of  this  separation  is  in  the  case  of  D  *  0.1  km.  The  variation 
in  angle  0  ,  shown  in  Fig.  3-36  is  similar  to  the  field  strength  vari- 

cl 

ation,  while  the  polarization  is  essentially  constant. 

F.  THE  CHARACTERISTICS  OF  A  WAVE  ARRIVING  FROM  THE  MAGNETIC  FIELD 
DIRECTION 


A  set  of  calculations  was  carried  out  for  a  distance  of  47.23  km  for 
sea  water  and  summer-night  model  ionosphere.  At  this  distance  the 
n  =  0  ray  is  parallel  to  the  earth's  magnetic  field  at  the  geomagnetic 
latitude  of  Stanford  (44°N).  This  was  done  for  the  particular  purpose 
of  determining  the  characteristics  of  a  wave  arriving  from  the  magnetic 
field  direction.  Such  a  signal  would  actually  be  made. up  of  the  n  =  0 
ray  plus  other  higher -order  rays  of  significant  amplitude. 

The  fields  are  shown  in  Fig.  3.64  as  a  function  of,  frequency.  The 
patterns  for  the  NS  (normal)  and  EW  (abnormal)  loops  are  similar,  with 
the  EW  loop  voltage  about  '20$  less  than  the  NS  loop.  The  pattern  of  the 
vertical-incidence  field  variation  is  quite  well  retained,  with  the 
peaking  near  the  frequencies  for'  which  h'  is  an  odd  multiple  of  one- 
half  wavelength.  The  apparent  polarization  is  almost  constant,  as 
shown  by  Fig.  3*65.  The  Lissajous -pattern  ellipse  has  an  axis  ratio  of 
about  O.85  -  0.9  and  its  major  axis  is  oriented  within  +  15°  from  the 
vertical  for  0  '  equal  to  zero. 


The  behavior  of  these  curves  It  typical  of  a  vide  range  of  ionosphere 
and  ground  conditions.  It  Is  thus,  apparent  that  If  a  whistler  came  from 
the  field  direction,  It  would  have  the  characteristics  of  short -distance 
propagation,  with  the  distinctive  apparent  polarization  pattern  and 
variations  in  field  strength.  The  latter  should  be  apparent  on  spectro¬ 
grams.  In  addition,  the  vertical -antenna  voltage  would  be  only  about 
half  the  loop-antenna  voltage.  These  facts  make  it  quite  apparent  that 
only  a  very  few  whistlers  have  the  field  characteristics  compatible  with 
arrival  from  the  magnetic  field  direction,  and  none  has  been  observed 
with  the  proper  apparent  polarization.  It  is  concluded  therefore  that 
whistlers  do  not  arrive  from  the  direction  of  the  magnetic  field. 

G.  WHISTLERS  WITH  FAIRLY  CONSTANT  APPARENT  POLARIZATION 

A  few  whistlers  have  been  observed  with  fairly  constant  apparent 
polarization.  The-  crossed-loop  Lissajous  patterns  were  observed  visually 
by  the  author  at  Stanford  on  27  May  1959,  about  0730Z.  All  showed 
elliptical  patterns  on  the  crossed  loops.  The  axis  ratio  varied  from 
about  0.1  to  0.9.  The  major  axis  tended  to  be  oriented  at  about  U50 
clockwise  from  north.  There  were  some  rapid  variations  in  the  pattern, 
but  the  visual  integration  effect  made  it  difficult  to  notice  details. 

The  exit -point  of  these  whistlers  was  probably  nearby,  to  the  northwest 
or  southeast.  Their  pattern  orientations  were  not  compatible  with  the 
possibility  of  arrival  from  the  magnetic  field  direction.  One  almost 
linear  pattern  was  observed  (axis,  ratio  about  0.1)  which  was  oriented 
almost  due  east-west.  An  explanation  for  this  might  be  a  long  distance 
with  low  ionosphere  reflection  coefficients  and  poor  ground  conductivity. 
None  of  the  available  calculations  explain  this  observation. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCUSSIONS 

The  ray -theory  method  used  in  this  report  utilizes  the  following 
simplifying  assumptions:  (1)  each  ray  is  a  plane  wave  in  the  antenna 
area,  (2)  the  ql  approximation  is  valid  (3)  the  ionosphere  may  be  repre¬ 
sented  by  an  equivalent  homogeneous  model  with  a  sharp  lower  boundary, 

(.4)  the  necessary  cone  of  plane  wave  energy  is  available  for  excitation, 
(5)  all  the  rays  are  of  a  single  frequency.  This  simplified  method  of 
adding  multiple  ray.s  can  still  account  for  most  of  the  previously  un¬ 
explained  characteristics  found  in  a  variety  of  experimental  data.  These 
include,  the  variations  with  f  requency  in  the  shape  and  position  of  the 
.crossed-loop  Lissajous  patterns,  the  changes  in  the  direction  of  rotation 
of  these  patterns  and  the  wide  variation  with  frequency  of  the  goniometer 
null  bearings  and  ratios.  This  agreement  with  experimental  data  is  not 
on  a  point -for -point  basis.  The  behavior  under  specified  conditions 
is  an  extremely  complicated  function  of  all  the  variables.  It  is  possible 
that  a  correlation  on  a  point -for -point  basis  might  be  found  between  a 
particular  whistler  and  a  particular  set  of  conditions.  The  investi¬ 
gation  of  this  possibility  would  be  far  too  costly,  and  the  required 
conditions  might  not  be  unique.  The  agreement  is  in  the  types  of  be- 
.  havior,  particularly  that  of  the  crossed-loop  Lissajous  patterns. 

This  type  of  behavior  was  not  explained  prior  to  this  work. 

The  results  show  that  whistlers  do  not  generally  arrive  from  the 
direction  of  the  magnetic  field,  contrary  to  the  conclusion  which  has 
sometimes  been  reached  from  earlier  theories.  Any  whistler  which  does 
arrive  from  the  magnetic  field  direction  can  be  identified  by  its 
characteristic  variations  in  field  strength  and  crossed-loop  Lissajous 
pattern  with  frequency.  The  results  are  also  consistent  with  the 
hypothesis  that  whistlers  emerge  from  a  relatively  small  area  in  the 
lower  edge  of  the  ionosphere,  and  that  this  exit  point  may  be  at 'a  con¬ 
siderable  distance  from  the  receiver.  The  average  difference  between 
the  measured  d-f  latitudes  and  those  calculated  independently  from  the 
dispersion  characteristics  is  roughly  what  would  be  expected  from  the 


calculated  bearing  errors.  Therefore,  the  exit  point  location,  end 
presumably  the  path  location,  la  veil  described  by  this  calculated 
magnetic  latitude.  The  occasional  large  discrepancies  between  the  A -t 
and  dispersion  results  can  be  accounted  for  by  d-f  error. 

The  variation  in  the  observed  bearings  with  frequency  does  not 
indicate  a  variation  in  the  exit  point  location,  since  the  variation  in 
bearing  is  consistent  with  that  calculated  for  a  fixed  location-. 

Consideration  of  the  series  of  rays  which  determines  the  field  com¬ 
ponents  for  the  vertical  loops  and  monopole  yields  an  explanation  for 
the  phenomenon  of  apparent  lateral  deviation  at  vlf .  No  actual  lateral 
deviation  need  be  present.  The  same  effect  produces  the  difference  in 
field  strength  measured  from  a  normal  loop  and  from  a  vertical  monopole, 
which  is  a  previously  unknown  factor  in  the  error  in  field-strength 
measurement  at  vlf.  A  consideration  of  the  total  wave  field  present 
should  be  made  before  interpreting  measurements  made  with  antenna 
systems  measuring  more  than  one  component  of  the  total  field.  The 
coupling  of  cross -polarized  components  due  to  the  earth's  magnetic  field 
is  important  in  this  case.  Care  must  be  exercised,  when  using  the  results 
of  measurements  on  fields  generated  from  a  single  polarization  component, 
to  predict  the  characteristics  of  a  wave  in  which  both  components  are 
originally  excited. 

The  results  of  the  calculations  show  that  a  vertical  monopole  is 
as  effective  as  a  vertical  loop  for  whistler  signals  coming  from  all 
distances  beyond  about  60  km,  which  includes  almost  all  whistlers. 

The  orientation  of  a  loop  antenna  is  not  important  at  the  low  frequencies 
(around  5  kc)  but  can  reduce  the  received  signal  by  at  least  6  db  at  the 
higher  frequencies  (around  15*5  kc).  The  radius  of  the  effective-area 
circle  for  a  whistler,  for  20  db  reduction  in  signal  from  the  center  to 
the  edge,  is  shown  by  the  calculations  to  be  approximately  1000  km. 

This  value  is  consistent  with  experimental  data.  The  theoretical  re-  . 
suits  for  the  goniometer  and  goniometer-vertical  systems  show  .that  the 
null  bearing  error  in  the  case  of  a  single  ray  is  much  smdller  for  the 
goniometer -vertical  system  than  for  the  goniometer.  The  calculations 
show  that  the  actual  null  bearing  error  of  the  goniometer-vertical  in 
the  case  of  a  series  of  incident  rays  is  considerably  greater  than  that 
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of  the  goniometer.  This  results  freo  the  fact  that  the  parallel  component 
of  the  field  seen  by  the  loop  is  different  in  both  amplitude  and  phase 
from  that  seen  by  the  vertical  monopole.  The  calculated  value*  of 
bearing  error  show  that  the  two  systems  have  nearly  the  fame  average 

i 

error  when  the  nulls  are  measured  properly.  The  theoretioal  increase  in 
null  ratio  of  the  goniometer-vertical  system  over  the  goniometer  for  a 
single  ray  is  also  realized  for  multiple  rays.  The  calculated  results 
show  that  the  signal  strength  tends  to  be  appreciably  greater  over  a 
sea-water  path  compared  with  a  path  over  poor  ground  at  frequencies  in 
the  15  -  20  kc  range.  This  increase  may  amount  to  as  much  as  1$  db  for 
long  distances  (2000  km).  At  the  lower  frequencies,  the  difference  is 
relatively  small,  being  about  10  db  at  5  kc. 

The  results  show  that  the  best  d-f  system  would  be  one  that  utilized 
identical  antennas.  In  this  case  the  amplitude  of  the  signal  would  be 
independent  of  the  bearing  of  the  signal.  The  bearing  would  therefore 
depend  only  upon  the  relative  phases  of  the  antenna  voltages.  A  system 
that  fits  these  requirements  is  an  Adcock  array,  made  up  of  four  vertical 
monopoles  arranged  in  a  square.  The  difference  between  the  signals  from 
antennas  in  opposite  corners  is  taken,  and  the  two  resulting  differences 
applied  to  a  goniometer  to  determine  the  bearings. 

The  errors  in  a  properly  adjusted  Adcock  array  are  more  subtle  than 
those  in  a  loop  system.  The  polarization  and  lateral -deviation  errors 
are  eliminated.  The  antennas  respond  only  to  the  parallel  component  of 
the  field  seen  by  a  vertical  monopole.  The  difference  in  the  signals 
obtained  from  a  pair  of  antennas  in  the  plane  of  an  arriving  plane  wave 
is  thepretically  zero.  Any  null  bearing  error,  therefore,  must  arise 
either  from  a  non-plane  wave  front  in  the  antenna  area,  or  from  site 
errors. 

There  are  numerous  practical  problems  involved  in  the  Adcock  system. 
The  antenna  spacing  must  be  large  enough  to  make  voltages  large  compared 
to  system  noise.  Small  antenna  spacings  increase  the  relative  errors 
due  to  unbalance  in  the  system. 
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3.  RECOMMENDATIONS  FOR  FUTURE  WORK 


1.  Further  study  of  ray  theory  la  needed  In  order  to  establish  in 

a  more  quantitative  way  the  regions  of  validity*  The  possible  advantages 

of  mode  theory  for  the  study  of  whistler  propagation  should  be  studied. 

•  • 

There  is  great  need  for  a  definitive  work  on  this  subject  which  states 
the  restrictions  on  each  method.  More  detailed  physical  interpretations 
of  the  ray  and  mode  theories  Beem  to  be  needed. 

2.  A  more  accurate  model  of  the  lower  region  of  the  ionosphere 
should  be  developed.  The  use' of  a  homogeneous  model  with  a  sharp  lower 
boundary  at  vlf  is  a  rather  rough  approximation  and  probably  results  in 
significant  errors  for  some  conditions  of  interest,  particularly  at  low 
frequencies.  Consideration  should  be  given  to  the  use  of  a  two  or  three 
slab  model  to.  account  for  the  increased  penetration  at  low  frequencies 
and  for  absorbing  regions  below  the  reflection  point. 

3.  The  characteristics  of  propagation  for  conditions  under  which 
the  ql  approximation  is  not  valid  should  be  studied  further.  This  should 
probably  be  dene  in  combination  with  the  multi-slab  model,  since  both 
approximations  are  not  necessarily  applicable  under  some  conditions  of 
interest. 

4.  In  problems  for  which  the  sharp -boundary  homogeneous  model  and 
ql  approximation  are  valid,  better  values  of  the  parameters  csr  and  t 
are  needed.  It  is  suggested  that  total  field  information  be  obtained 
from  ground-to-ground  measurements  to  provide  an  experimental  determin¬ 
ation  of  T. 

5.  The  same  ray -theory  methods  may  also  be  used  to  calculate  the 
excitation  of  a  whistler  duct  by  a  series  of  rays  transmitted  into  the 
ionosphere  from  the  waveguide.  The  results  would  be  valuable  for  the 
determination  of  transmission  losses  and  duct  sizes  in  whistler -mode 
propagation. 

6.  '  The  effect  of  the  difference  in  the  frequency  of  the  rays  for 
.a 'whistler  at  near-vertical  incidence  deserves  more  investigation.  The 
steady-state  calculations  show  regular  sharp  maxima  and  rather  broad 
minima  at  harmonic  frequencies.  The  introduction  of  different  frequencies 
would  undoubtedly  change  the  shape  of  and  displace  the  positions  of  the 
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maxima  and  minima.  It  would  also  produce  an  amplitude . fluctuation  of  a 
few  db  as  a  result  of  beating  of  different  frequencies.  It  could  be 
studied  by  the  present  method  if  the  exact  frequency  were  used  for  each 
ray  end  a  very  large  number  of  calculations  were  made  at  very  slightly 
different  frequencies.  The  cost  of  the  calculation  would  be  prohibitive 
at  present.  It  is  probable  that  useful  information  could  be  obtained 
by  the  use  of  a  simple  analog  device,  consisting  of  a  shorted  trans¬ 
mission  line. 

7.  Further  studies  of  whistlers  with  the  present  antenna  system 
are  not  warranted.  The  accuracy  is  too  limited  to  be  of  further  interest. 
The  development  of  an  Adcock  system,  as  previously  described,  would  be 
of  great  interest.  The  increased  accuracy  of  each  bearing  and  the  in¬ 
crease  in  the  number  of  whistlers  showing  nulls  should  make  the  project 
well  worth  the  development  effort.  The  system,  could  also  be  applied  to 
the  study  of  whistler-mode  echoes  from  fixed-frequency  vlf  transmitters. 
Presentation  of  the  Lissajous  patterns  from  the  two  pa^rs  of  antennas 
would  provide  a  presentation  like  that  of  the  standard  d-f  trace  which 
would  be  free  of  polarization  error . 

A  recording  technique  using  an  intensity -modulated  sweep  triggered 
by  the  goniometer  indexing  pulse  would  probably  enable  direction  measure¬ 
ments  to  be  made  on  vlf  transmitters  when  ordinary  cw  traffic  is  being 
transmitted.  This  would  permit  continuous  recording  instead  of  requiring 
special  pulse  programs  as  at  present. 
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APPENDIX  A;  PATH  GEOMETRY  EQUATIONS 


The  path  geometry  equations  used  in  this  report  are  slightly  modified 
versions  of  those  used  for  the  usual  case  of  ground-to-ground  trans¬ 
mission.  The  expressions  for  a  single  trip  between  earth  and  ionosphere 

20 

are  given  by  Helliwell  in  this  form: 

P2  =  2R(h'  +  R)[l  -  cos(l80D/jtR)]+h'2  (A.l) 

9n  =  tan’^sinUSOD/jtR)/!  +  h'/R  -  cos  (l80D/«R)]  (A. 2) 

Pn  =  90  '  9n  "  (l80D/ItR)  (A.  3) 

where 

P  =  ray  path  distance  for  a  single  trip  from  ionosphere  to  ground, 
n 

9  =  angle  of  incidence  at  the  ionosphere  in  degrees, 

n 

Pn  =  angle  of  take-off  or  arrival  at  the  earth,  in  degrees, 
h'  =  virtual  height  of  the  ionosphere  in  km. 

D  =  total  ground  distance  in  km. 

R  =  mean  earth  radius,  in  km,  equal  to  6367*39  km. 

These  equations  must  now  be  generalized  for  the  case  of  n  reflections  from 
the  earth  and  ionosphere,  with  transmitter  located  at  the  ionosphere, 
and  receiver  at  the  ground.  The  situation  is  illustrated  in  Fig.  2.1, 
where  the  case  of  n  =  2  is  illustrated.  It  is  readily  seen  that -the 
equations  may  be  generalized  to  this  case  by  merely  dividing  the  total 
distance  by  the  factor  (2n  +  1).  The  equations  then  become: 

p2  =  2R(h'  +  R)  {l  -  cos[l80D/(2n+l)nR]  }  +  h'2  (A.4) 

0  =  tan-1  {sin[l80D/(2n+l)nR]/l+h'/R-cos[180D/2n+l)jtR]}  .  (A.5) 

=  90  -  0n  -  l80D/(2n+l)nR  (A.6) 
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Then  the  total  path  distance  R 

n 


and  the  phase 


A 


are  given  by 


R  =  (2n+l)P 
n  *  n 

y/T^  =  -  Rn  ^  degrees 


(A.7) 


These  equations  yield  -11  the  necessary  information  about  the  path  geo¬ 
metry  for  solution  of  the  problem. 


APPENDIX  B.  THE  REFLECTION  COEFFICIENTS 

The  ionospheric  reflection  coefficients  uead  in  this  report  are 

12  13 

those  derived  by  Budden  and  are  listed  in  the  fora  given  by  Wait: 


II  3,1  =  [(4o+Ue)(C2-C0Ce)+(noHe-l)(Co+Ce)C]/D  (B.l) 

•  „RA  -  2JC(noCo-ueCe)/D  (B.2) 

lR..=  2jC(%Ce  <B*3) 

iRi  -  [(^o%)(C2-C0Ce)-(no4e-l)(Co+Ce)C]/D  .  ....  (B.L) 


where 

D  =  ( u  +u  ) ( C2+C  C  ) +( u  H  +1 ) ( C  +C  ) C 
06  oe  06  0  6 

C  =  cos  9 

n 

C  =  cos  9 
o  o 

C  =  cos  9 
e  e 

The  subscripts  o  and  e  refer  to  the  ordinary  and  extraordinary 
waves,  respectively.  As  in  Chapter  II,  the  whistler  mode  is  called  the 
ordinary  wave,  in  keeping  with  the  notation  and  derivation  used  in  the 
references. 

The  ground  reflections  are  expressed  by  the  Fresnel  coefficients 
42  43 

as  given  by  Wait  and  Jordan,  with  the  take  off  or  arrival  angle  Pn 

substituted  in  the  proper  manner  for  the  angle  of  incidence.  This  substi 
tution  yields : 

=  sinPn+(N2-cos2pn)1^2/sinPn-(N2-cos2pn)1//2 

'  R  „  =  N2 sinf3  -(N^cosS  J^/l^Binp  +(N2-cos2p  )1/2 
11  n  n  '  a  n 

where 


(B.5) 

(B.6) 
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^  a  ^conductivity  a t  the  ground  la  mho/m. 

"  k/k  =' dielectric  constant  of  the  ground. 

'  o 

kq  =  permitivity  of  free  space  in  fd/m. 
f  =  wave  frequency' in  cps. 

The  other  conventions  such  as  the  field  component  indicated  are . the  same 
as  those  in  previous  sections. 


APPENDIX  C.  THE  COMPUTER  PROGRAM 


The  saleulatioos  of  ‘the  total  field  and  the  antenna  response  were 
carried  out  at  first  on  the  ISM  6^0  digital  computer  at  the  Stanford 
Computation  Center*  Tie  Bell  system  with  a  complex  arithmetic  subroutine 
was  used.  The  running  times  were  found  to  average  at  least  2  minutes 
per  n,  where  n  is  the  number  of  reflections  for  each  reflection  from 
the  earth  and  ionosphere.  The  total  time  per  calculation  was  therefore 
very  long,  the  minimum  being  about  20  minutes.  A  hand  calculation  was 
performed  as  a  check  on  the  programing  and  operation.  A  standard  test 
program  was  run  regularly  as  a  check. 

During  the  calculations  a  Burroughs  220  computer  became  available.  • 
The  program  was  then  re-written  for  the  220,  using  the  "Bleap"  assembler 
system.  A  complex  arithmetic  subroutine  was  written  for  Bleap  for  this 
project.  Intermediate  data,  including  all  the  important  path  parameters, 
the  results  for  each  value  of  n,  and  the  final  results,  were  recorded  on 
magnetic  tape  by  the  program.  These  results  could  later  be  taken  off 
the  tape,  the  decimal  point  fixed,  and  cards  punched  for  use  on  an  X-Y 
plotter^  Considerable  time  was  saved  by  reproducing  most  of  the  curves 
directly  on  Multilith  mats  with  this  machine. 

The  new  program  gave  an  approximate  speeding-up  factor  of  thirty 
to  one  in  calculating  time.  The  running  times  range  from  about  fifty 
seconds  for  n  =  4-  to  about  400  seconds  for  n  =  50.  The  program  can 
calculate  a  maximum  n  equal  to  97. 

The  programing  was  organized  in  a  system  of  subroutines  to  enable 
the  use  of  these  subroutines  for  other  programs,  if  desired. 

The  calculations  of  the  curves  for  the  phase  shift  ellipse  conver¬ 
sion  were  performed  on  the  IBM  650,  again  using  the  Bell  system  routine. 
The  calculations  of  theoretical  null  ratio  and  bearing  error  were  also 
done  in  the  same  system  on  the  650. 

The  d-f  location  program  was  aided  by  another  pair  of  routines, 
also  written  in  the  Bell  system  for  the  65O.  These  routines  solved  the 
spherical -trigonomeJixy . equations  to  determine:  (l)  the  bearing  and 
distance  between  any  pair  of  points  on  the  earth's  surface,  when  the 
latitude  and  longitude  are  given;  and  (2)  the  latitude  and  longitude 


of  the  "target"  when  the  latitude  and  longitude  of,  and  the  hearing 
from,  up  to  five  stations  are  given.  The  equations  are  solved  for  each 
pair  of  stations  giving  bearings j  for  example:  one  location  with  two 
stations*  three  with  three  stations;  six  with  four  stations,  etc.  This 
■olutioa  was  always  accurate  to  at  least  the  nearest  one-tenth  degree, 
so  that  the  location  accuracy  was  never  limited  by  these  calculations. 
The  running  times  were  5-10  seconds  and  about  20  seconds  for  the  two 
respective  routines. 
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APPENDIX  D.  THE  PHASE  SHIFT  ELLIPSE  CHARACTERISTICS 


A  common  problem  in  the  use  of  Lissajous  pattern  ellipses  Is  that 
of  converting  the  pattern  specified  by  an  axis  ratio  6  and  an  orientation 
of  the  major  axis  into  the  equivalent  complex  form:  ^  ■  X  +  JY.  In 
order  to  facilitate  this  conversion,  a  pair  of  conversion  formulas  were 
derived  in  the  following  manner: 

Let  the  vertical  . nd  horizontal  deflections  be  given  by: 

y  =  A  sin  (D.l) 

x  =  A  sin  (oj  +  0)  (D.2) 

Equation  (D.2)  can  be  expanded  to  yield: 


x  =  A  (sin  cu  cos  0  +  cos  cu  sin  0)  (D.3) 

X  Z  X 

Now,  substituting  (D.l)  in  (D.3)  and  substituting  R  =  A J k,  rearranging 
and  squaring,  yields: 


^y  cos  0  +V  1  -  ^L.  '  sin 


A  =  A  I  y  cos  6 
x  x  v 


A2x  -  AA  y  cos  0  =  A  -  y2>  sin  0 


A2(x  -  R  cos  0)  =  AR  VA2  -  y2  sin  0 

y 

p  p  p  p  p  p 

A  (x  -  R  cos  e)  =  R  sin  0(A  -  y  ) 

Substituting  Q  =  X  +  JY  =  R(cos  0  +  J  sin  0)  yields: 


Y2(A2  -  y2)  =  (x  -  Xy)2 


(D.5) 


It  can  be  shown  that  Eq.  (D.5)  describes  an  ellipse;  so  the  next 
step  is  to  rotate  the  coordinate  axes  so  that  the  major  axis  is  aligned 
with  the  x  axis.  A  general  rotation  of  0  degrees  is  expressed  by: 

x  =  x'cos  0  +  y'  sin  0 

y  «  x'  sin  0  +  y'  cos  0  • 
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After  this  substitution  is  made*  the  terms  x  and  y  are  understood  to 
mean  the  new  coordinate  system,  which  yields: 

x2  ( -Y2 sin20  -  cos20  -  2X  BlafS  coe0  -  X2 sin20) 

+  xy  ( +2Y2sin0cos0  -  2sin0coa0  +  2Xa©s20  -  2Xsi a20  +  2X2sin0cos0) 

+  y2  ( -Y2cos20-sin20  +  2Xsin0cos0  -  X2cos20)  -i-  A2Y2  0  (D.6) 

Now,  to  determine  the  angle  of  rotation,  the  xy  term  is  equated  to  zero, 
which  puts  the  equation  almost  in  standard  form.  This  step  yields: 

Y2sin20  -  sin20  +  2Xcos20  +  X2sin20=  0 

sin20(X2  +  Y2-l)  =  -2Xcos20 

oy 

tan20  =  — -£--p  (D.7) 

1-X  -x 

The  tilt  angle  0  of  the  major  axis  is  determined  from  Eq.  (D.7).  The 
next  step  is  to  find  the  functions  of  20  in  terms  of  X  and  Y,  which 
yields : 

'  2  1  1  (l-X2 -Y2)2  (l-X2^2)2 

cos  20  =  - ^ -  =  — - - - (D.8) 

sec  20  tan  20+1  4X+(l-0T-Y^r  D 

2  1  1  4X2  4X2 

sin  20  =  ■  2  =  2  =  2  ~2  2~ 2  ~  (®«9) 

csd  20  ctn  20+1  4X  +(1-)T-Y  )  D 


where  D  =  4X2+(l-X2-Y2)2 

V/hen  Eq.  (D.8)  and  Eq.  (D.9)  are  substituted  into  Eq.  (D.6)  the  result 
is: 


x2[-Y2(^-|cos20)-(|+icos20)-Xsin20  -  X2(|4cos20)  ] 

+  y2[-Y2(  |+icos20)  -  ($-$cos20)  +  Xsin20  -  X2(£+|cos20)  J+A^  =  0 
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* 


2 

x 


2  yr  J 

P  P  P  P 

Yd  r(i-3r.r) 
2  #  2  fp 


(l-y?-Y2)  2X2  X2  X^l-X2-^) 

2  f¥*  YF  2  VT 


(l-^-Y2) 

z'fiT 


2X2  X2  X^d-X2-^) 

flT  2  2  ITd1 


+  A2^  «  0 

2  2 

The  results  must  now  he  manipulated  into  the  standard  form  of  x  /a 

2  2 

+  y  /b  =1  for  an  ellipse  where  a  and  b  are  the  axes.  This  operation 
yields: 


2  2 
AY 


(1+X^-fY2)  y  4X2+(1-X2-Y2)2  +'  [l+X^l-X2.*2)2] 


2  2 

aV 


(1+X2+Y2)  y  4x2+(1-X2-Y2)2-  f 4X2+( 1 -X2 -Y2 ) 2 ] 


The  axis  ratio  e  may  then  be  determined  from: 


2 

€ 


a2  (1+X2+Y2)  -  V^d-X^Y2)2  ' 
b  (1+X2+Y2)  +  y4X2+(l-X2-Y2)2' 


(D.10) 


These  two  equations  now  make  simple  the  conversion  from  the  complex  form 
to  the  ellipse  constants.  The  reverse  operation  is  still  difficult.  The 
best  solution  is  to  plot  the  curves.  The  result  is  Fig.  D.l,  which  shows 
curves  of  constant  X  and  Y  as  a  function  of  €  and  0.  Note  that  the  curves 
of  constant  X  and  Y  are  nearly  orthogonal  for  most  conditions  and  one 
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FIG.  D. 1.  PHASE  SHIFT  ELLIPSE  CONVERSION  CHART. 


sheet  can  be  used  for  both  directions  of  conversion.  3he  result  is 
plotted  for  positive  X  and  Y.  It  is  evident  from  Eq.  (D.7)  that  a 
change  in  the  sign  of  Y  has  no  effect  on  $  or  €,  and  a  change  in  the 
sign  of  X  merely  changes  the  sign  of  0.  This  curve  is  therefore  usable 
for  all  cases. 


APPENDIX  E.  THE  TRANSMISSION  COEFFICIENT  APPROXIMATION 


i 


The  derivation  of  the  transmission  coefficients  for  the  V  field  of 

vhiatler-ooda  waves  traveling  from  the  Ionosphere  Into  free  space  is 

given  la  Section  II-E,  page  16.  She  derivation  uses  two  facts  which 

follow  from  the  ql  approximation.  The  practical  necessity  for  using  the 

ql  approximation  arises  from  the  fact  that  the  refractive  index  in  the 

general  case  is  a  function  of  the  wave  normal  angle.  This  oase  would 

17 

require  solution  of  a  complicated  equation  for  each  transmission  angle, 
which  would  complicate  the  solution,  for  several  rays,  to  an  impractical 
point.  Hie  ql  approximation  was  shown  to  be  a  valid  approximation  for 
most  of  the  whistler  range. 

When  the  ql  assumption  is  applied,  the  refractive  index  becomes 
independent  of  the  wave  normal.  As  a  consequence,  the  two  characteristic 
waves  are  circularly  polarized  with  opposite  rotation.  It  is  clear  that 
for  vertical  incidence,  the  ordinary -wave  whistler  energy  would  excite 
no  reflection  of  the  other  characteristic  wave.  For  angles  other  than 
the  vertical,  the  amount  of  the  other  wave  excited  should  still  be  small. 
It  was  thus  assumed  that  the  field  of  the  other  characteristic  reflected 
wave  would  be  negligible.  A  measure  of  the  quality  of  this  approximation 
may  be  determined  by  finding  the  corresponding  reflection  coefficients 
for  the  assumed  ordinary  wave.  Since  the  ••  and  i  field  components  in 
the  incident  wave  are  equal  in  magnitude,  the  actual  polarization  of  the 
reflected  ordinary  wave  calculated  under  these  assumptions  will  be 
determined  by  the  ratio  of  the  reflection  coefficients,  which  we  shall 
label  Rj_  and  R(|  .  Rj^  may  be  easily  found  from  Eqs.  (2.17)  and  (2.20a), 
and  Rj^  from  Eq.  (2.18)  and  Eq.  (2.19a). 


Ri  - 


Ei  [i  C  ■  C , 
ir  _  00  '  t 

E  -  ,  |i  C  +  G . 
li  00.  t 


(E.l) 


!lk  .  co  ~  Moct 

*1li  Co  +  ^oCt 


(E.2) 
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The  ratio  is  then: 


*1  <“oCo  -  ct>  <Co  *  ^ 

R»  ’  *  <V  (°o  •  “0ct> 


|1J  C  0,  -  tl  cf  +  |1  C8  -  0  c 
0  O  t  *0  t  O  0 


O  t 


-  n2c  ct  ■  (i  cf  +  n  cs  +  o  a 

0  O  t  O  t  0  o  o  t 


c  c.  (ii‘  -l)  +  n  (cr-cf) 

o  t'  o  '  *o  o  %' 


C  C. (n2-l)  -  H  (c2 
O  t  o  f  O 


c?> 


(*.3) 


The  magnitude  of  this  quantity  will  be  very  close  to  unity  for  large 
(iQ  and  angles  of  incidence  not  too  large.  The  range  of  validity  there¬ 
fore  is  approximately  that  of  the  ql  approximation. 


t 


APPENDIX  F.  THE  CRITICAL  ANGLE  EFFECT  IN  THE  IONOSPHERE 

i 


The  ionosphere  is  assumed  in  this  report  to  be  homogeneous  with 
a  vertical  magnetic  field,  with  losses  represented  by  v,  the  effective 
collisional  frequency.  Not  very  much  is  known  about  refraction  and 
tranamission  at  vlf  in  such  a  medium.  The  assumption  of  infinite  plane 
waves  is  not  strictly  correct,  because  of  the  loss  factor. 

The  refraction  properties  of  such  a  medium  are  treated  by  Stratton. 
He  considers  the  case, of  transmission  from  a  dielectric  medium  (2)  into 
a  medium  possessing  conductivity  (l).  The  wave  which  is  transmitted 
into  the  refracting  medium  is  what  he  calls  an  "inhomogeneous  plane 
wave",  in  which  the  planes  of  constant  amplitude  are  parallel  to  the 
interface.  The  normals  to  the  planes  of  constant  phase  which  define  the 
direction  of  propagation,  are  inclined  at  an  angle  f  to  the  normal  to  the 
interface.  The  angle  f  is  defined  by: 


cos  -ft  = 


sin 


3 

o 


t 


4 


sin  1 1/ 


sin  0 

_ o _ 

+  OU  sin  0 
d  O 


where  0  = 


o 


constant  in 


cos©^  = 
P  = 

<1  = 


the  ang|Le  of  incidence  in  medium  (2) 

the  angle  of  refraction  in  medium  (l)  as  given  by  Snell's 
law 

2  2 

a)  +  icuo^i  =  (o^  +  JP^)  ,  the  square  of  the  propagation 
medium  (l) 

2 

a)  the  square  of  the  propagation  constant  in  medium  (2). 

pe1*  .  Vl  -  sinV 
p(P1  cosy  +  o^siny) 
p(a1  cosy  “  f^sinY) 


y 


-  1 66  - 


r 


* 


The  above  equation  from  Stratton  may  be  simplified  by  further  sub¬ 
stitution.  The  first  step  is  to  express  tan  f  as  follows. 


,  OLsin  0 
,  sin  $  2  o 

tan  ib  =  - 7  =  - 

cos  if 


(F.l) 


The  quantity  q  is  easily  seen  to  -be  the  real  pn^t  of  k^  cos  0^. 
Therefore : 


tan  1 1/  = 


0!  sin  0 

_2 _ o 

Re^cose^^ 


c*  sin  9 
2  o 


Rek,  i/l  -  sin 


a.  sin  9 
2  o 


✓  2  2  2 
k^  -  k^  sin  0^ 


(F.2) 


Now,  by  the  use  of  Snell's  law,  this  yields: 


tan  if/ 


OL  sin  0 
2 _ o 


(F.3) 


Now,  substitute  for  the  propagation  constants  k  =  —  ,  where  u  is  the 

c 

refractive  index  of  the  medium.  If  medium  (2)  is  free  space  and  we  assume 
that  medium  (1)  is  the  ionosphere,  then  for  the  whistler  mode: 


-  -  sin  0 

tan  -  - - - — — — - 

Re  A/  \J?  -  sin  0 

This  expression  is  readily  evaluated  by  using  the 
Eq.  (2.9)  for  uq. 


(F.lf) 


ql  approximation  and 
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FIG.  F.l.  CRITICAL  ANGLE  VS  FREQUENCY  FOR  VARIOUS  IONOSPHERE  MODELS 
AND  ANGLE  OF  EARTH’S  MAGNETIC  FIELD  FROM  VERTICAL  VS  GEOMAGNETIC 
LATITUDE. 
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The  critical  angle  is  given  when  9q  =  90°,  and  is  plotted  in 

Fig.  (F.l)  for  the  summer-night  model  ionosphere,  the  summer-day  model, 

and  two  variations  of  o>  from  the  summer-night  model.  In  addition,  the 

r 

angle  of  the  earth's  magnetic  field  measured  from  the  vertical,  is 

shown  as  a  function  of  geomagnetic  latitude  A. 

This  figure  shows  several  interesting  features.  The  summer-night 

model  curve  intersects  the  magnetic  field  curve  at  f  =  10.6  kc,  and 

A  =  530  for  i/  =20.8.  Let  us  consider  the  effect  on  whistler  excitation, 
^cr 

Suppose  that  a  field-aligned  duct  is  present  which  may  trap  waves  with 

normals  reasonably  close  to  the  field  angle,  say  5  degrees.  Trapping 

would  then  be  possible  over  a  range  of  latitudes  above  29  degrees  and  for 

frequencies  above  5-3  kc.  If  the  trapping  angle  were  very  limited,  the 

trapping  could  be  confined  to  frequencies  and  latitudes  above  10.6  kc  and 
o 

53  ,  respectively.  In  the  opposite  case,  of  transmission  into  the  wave¬ 
guide,  this  effect  could  also  limit  the  number  of  rays  which  contribute 
to  the  field  at  an  antenna,  and  would  modify  Ibe  assumption  about  tte  exciting 
energy;  it  could  cause  complete  cutoff  in  some  cases.  The  trapping 
area  or  the  exciting  area  would  correspond  to  the  intersection  of  two 
cones  of  angles,  one  around  the  field  direction  corresponding  to  the 
trapping  angle,  and  the  other  around  the  vertical  corresponding  to  the 

angles  of  refraction.  Partial  intersections  would  cause  marked  directional 
* 

.  effects. 

Next  let  us  consider  the  effect  of  changes  in  the  ion  density.  Note 

that  a  decrease  in  the  ion  density  shifts  the  curve  upward,  as  shown  by 

4  o 

the  curve  for  =  5  x  10  ,  t  =  60  .  This  would  have  the  effect  of 
lowering  the  minimum  latitude  for  trapping.  A  decrease  in  the  ionization 
density  should  cause  a  lowering  of  the  minimum  latitude  and  of  the  average 

latitude  of  received  whistler  traces.  Such  a  decrease  in  density  may 

39 

occur  after  magnetic  storms,  and  the  number  of  low-latitude  whistlers 

1+0 

received  at  Byrd  Station,  Antarctica,  does  indeed  increase.  Other 
factors  could  cause  this  same  effect,  however. 


•ft 

Some  details  of  the  trapping  cones  and  the  directional  effect  were 
given  by  R.  A.  Helliwell  in  the  paper  "Synoptic  observations  of  whistler 
mode  signals  from  vlf  transmitters."  URSI-IRE  meeting,  Washington  D.C., 
May  1961. 
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